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FOREWORD 


Los  Angeles  Aircraft  Division  of  Rockwell  International  personnel  developed  the 
Mission  Damage  Effectiveness  Model  computer  simulation  under  Contract 
No.  1265R175900  issued  by  Armament  Systems  Incorporated,  Anaheim,  California.  The 
v work  was  performed  between  April  1975  and  August  1976.  Marvin  Gove,  Analysis  Branch, 
NWC  (Naval  Weapons  Center),  acted  as  contract  administrator  for  the  documentation  of  the 
program  under  NWC  Contract  N00123-76-C-0159. 

The  work  was  sponsored  by  the  JTCG/AS  as  part  of  a 3-year  TEAS  (Test  and 
Evaluation  Aircraft  Survivability)  program.  The  TEAS  program  was  funded  by  DDR&E/ 
ODDT&E.  The  effort  was  conducted  under  the  direction  of  the  JTCG/AS  Survivability 
Assessment  Subgroup,  as  part  of  JTCG/AS  Project  SA-6-02,  Trade  Studies. 

The  purpose  of  the  user  and  analyst  manuals  is  to  provide  a current  documentation  of 
the  methodology  and  easy  update  maintenance  for  future  program  applications  on  a 
page-by-page  basis. 

G.  L.  Gallien  and  S.  C.  Silver,  with  the  program  managership  of  R.  L.  Moonan,  were 
the  key  Rockwell  International  personnel  responsible  for  the  computer  program  develop- 
ment and  documentation. 
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NOTE 

This  technical  report  was  prepared  by  the  Survivability  Assessment  Subgroup  of 
the  Joint  Technical  Coordinating  Group  on  Aircraft  Survivability  in  the  Joint 
Logistics  Commanders’  organization.  Because  the  Services’  aircraft  survivability 
development  programs  are  dynamic  and  changing,  this  report  represents  the  best 
data  available  to  the  subgroup  at  this  time.  It  has  been  coordinated  and  approved 
at  the  JTCG  subgroup  level.  The  purpose  of  the  report  is  to  exchange  data  on  all 
aircraft  survivability  programs,  thereby  promoting  interservicc  awareness  of  the 
DOD  aircraft  survivability  program  under  the  cognizance  of  the  Joint  Logistics 
Commanders.  By  careful  analysis  of  the  data  in  this  report,  personnel  with 
expertise  in  the  aircraft  survivability  area  should  be  better  able  to  determine 
technical  voids  and  areas  of  potential  duplication  or  proliferation. 
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INTRODUCTION 


The  attrition  of  large  numbers  of  fixed  and  rotary  wing  aircraft  during  the  Southeast 
Asia  conflict  has  focused  attention  on  the  requirement  to  provide  design  features  that  will 
enhance  the  survivability  of  these  systems,  curtailing  the  losses  in  hostile  engagements. 
Implementing  the  design  features,  whether  in  initial  design  or  the  more  costly  retrofit 
programs,  usually  impacts  aircraft  weight  and  cost.  Because  of  these  impacts,  trade  studies 
are  required  to  develop  a prioritized  listing  of  survival  enhancement  features  that  increase 
mission  effectiveness.  The  priority  list  will  identify  those  features  that  provide  the  greatest 
increment  to  improved  effectiveness,  and  provide  the  data  base  from  which  a bound  can  be 
established  on  the  amount  of  hardening  that  is  effective  from  a mission-cost  standpoint.  The 
data  base  will  also  permit  a determination  of  the  best  mix  of  hardening  features  when  a 
budget  constraint  is  imposed. 

The  MISDKM  (Mission/Damage  Effectiveness  Model)  is  a simulation  of  an  aircraft  (and 
its  subsystems)  experiencing  a time-series  of  events.  The  events  may  include  targets  to  be 
attacked,  threat  weapons  to  be  encountered,  refueling,  recovery  at  an  air  base,  or  certain 
events  selected  by  the  user. 

MISDEM  evaluates  the  capability  and  effectiveness  of  an  aircraft  system  throughout  a 
mission  scenario  provided  by  the  user.  A schematic  illustrating  where  MISDEM  fits  into  the 
total  mission  effectiveness  analysis  procedure  is  displayed  in  Figure  1.  MISDEM  performs  a 
statistical  bookkeeping  function,  aggregating  the  detailed  results  of  systems  and  operations 
analyses  carried  out  for  various  threat  elements  and  target  combinations.  The  generation  of 
the  inputs  may  require  the  user  to  exercise  several  other  simulations.  For  example,  various 
options  of  aircraft  survivability  features  would  be  subjected  to  trade-off  studies  using 
models  such  as  GVAP,  GPBP,  SHOTC.EN.  MAGIC,  and  COVART  (JTCG/AS  approved 
programs).  An  end-game  simulation  such  as  ATTACK  would  then  be  utilized  to  produce 
threat  effectiveness  numbers  to  be  input  to  MISDEM.  The  output  of  the  model  is  several 
scalar  effectiveness  parameters  which  the  user/analyst  must  assemble  and  combine  to  make  a 
determination  of  mission  effectiveness  for  his  particular  problem. 

The  model  was  developed  initially  to  analyze  the  impact  on  system  survivability  and 
mission  effectiveness  of  hardening  various  aircraft  subsystems  to  the  weapon  effects 
produced  by  nuclear  weapons.  The  model,  as  presented  here,  has  been  extended  to  include 
nonnuclear  weapons  effect  capability.  The  approach  is  an  extension  of  concepts  developed 
by  the  WSEIAC  (Weapon  System  Effectiveness  Industry  Advisory  Committee).'  The  basic 
theory  of  MISDEM  has  been  described  by  Rockwell  International.2,3,4 


' Air  Force  Systems  Command  Weapon  Syatem  Effectiveness  Industry  Adviaory  Committee  (WSEIAC).  Final  Report 
of  Task  Group  IL  Andrews  AFB,  DC  January  1*65,  <AFSC-TR-«5-2  (Volume  II)). 

^Rockwell  International  (Los  Angeles  Aircraft  Division)  Description  of  an  improved  Effectiveness  Model, 
November  1973,  (TFD-74-62). 

^Rockwell  International  (Los  Angeles  Aircraft  Divirion)  Mission / Damage  Effectiveness  Model,  1974,  (NA-74-62). 
4RockweB  International  (Los  Angeles  Aircraft  Division).  Mission /Damage  Effectiveness  Model  (Sample  Case),  for 
Joint  Technical  Coordinating  Group/ Aircraft  Survivability.  Survivability  Assessment  Subgroup.  May  1974,  (NA-74-358). 
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Figure  1.  Schem Stic  of  Mission  Effectiveness  Analysis. 


Figure  2 is  a brief  summary  of  the  basic  inputs  and  outputs  to  the  MISDEM  simulation. 
The  mission  scenario  consists  of  a time/event  series,  in  which  the  events  are  either  offensive 
(aircraft  system  against  the  enemy)  or  defensive  (enemy  defensive  system  against  the  aircraft 
system).  The  aircraft  system  must  be  defined  in  terms  of  identifying  electronic  "black  box” 
functions  and  mechanical  functions  that  must  be  performed  at  specified  times  during  the 
mission  (e.g..  terrain  following  radar,  1R  target  acquisition,  engine  must  operate,  and 
weapon  must  launch).  Each  combination  of  electronic  and  mechanical  functions  define 
system  modes  of  operation.  The  optimum  mode  would  be  to  have  all  functions  available  and 
operating.  However,  due  to  enemy  damage  inflicted  on  the  aircraft  or  system  breakdown, 
some  subsystems  may  not  function  and  the  aircraft  must  operate  in  a degraded  mode.  An 
example  of  a degraded  mode  of  operation  would  be  the  pilot’s  use  of  a visual  navigation  fix 
plus  dead  reckoning  instead  of  utilizing  inertial  navigation  with  a radar  update.  Thus, 
M1SDEM  is  particularly  adapted  to  analyzing  aircraft  systems  that  are  multiply  vulnerable, 
although  singly  vulnerable  systems  are  readily  handled  by  the  probabilistic  nature  of  the 
model.  In  the  situation  where  several  functions  are  lost  during  the  mission,  the  user  may 
elect  to  abort  the  mission  and  determine  the  probability  that  the  aircraft  can  be  recovered.5 


5Jotat  Technical  Coordinating  Group/ Aircraft  Survivability.  MISDEM  Computtr  Smularton,  Kot  /.  Uter't  Manual. 
by  G.  L Gaittrn  and  S.  C Si  vet,  Rockwall  International,  Los  Angeles  Aircraft  DhrMon  Washington,  DC,  JTOG/AS,  (in 
process).  (JTCG/AS-76-S-003.  publication  UNCLASSIFIED.) 
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INPUTS.  SYSTEM  DESCRIPTION  AND  EVENT  DESCRIPTION 


Figure  2.  Mission  Damage  Effectiveness  Model  Simulation  Summary. 
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Each  defensive  event  requires  the  input  of  the  effectiveness  of  the  enemy's  defensive 
system  against  each  aircraft  component  /subsystem  for  each  mode  of  operation  After  each 
defensive  event,  MISDHM  evaluates  the  various  mission-related  subsystems  to  detine  aircraft 
system  survivability  (for  that  event)  and  the  systems  capability  for  performing  at  the  next 
event.  The  survivability  and  capabilit’  methodology  also  includes  the  loss  of  functions  due 
to  reliability  factors. 

The  offensive  events  require  the  input  of  our  weapon  effectiveness  for  each  target  type 
of  each  possible  mode  of  deliver,  (a  degraded  mode  of  delivery  may  be  the  use  of  fixed 
sights  instead  of  the  tire  control  computer).  After  each  offensive  event,  target  kill 
probability  is  evaluated  to  define  aircraft  system  effectiveness  for  that  event;  this  number  is 
combined  with  those  ol  preceding  offensive  events  to  define  a cumulative  mission 
effectiveness  update. 

I Inis,  MISDHM  keeps  hook  on  the  health  and  effectiveness  of  an  aircraft  system 
throughout  the  course  of  its  mission.  Output  consists  of  system  and  subsystem  survival 
probabilities  and  the  probabilities  of  the  system  operating  in  each  mode.  The  important 
calculations  performed  within  the  model  include  the  treatment  of  multiple  warheads  (the 
ettectiveness  inputs  lor  the  defensive  events  are  one-on-one  simulation  results  from  a model 
other  than  MISDHM)  and  the  rigorous  accounting  of  effectiveness  contributions  from 
alternative  modes  of  operation. 

The  basic  structure  of  MISDHM  requires  the  aircraft  system  being  studied  to  be  divided 
into  two  complementary  elements:  the  electronic  and  the  vehicle  functions.  The  MISDHM 
program  is  run  separately  lor  the  electronic  and  the  vehicle,  utilizing  different  inputs  for 
each. 


MISDHM  computes  the  progressive  degradation  in  vehicle  flight  status  and  electronics 
mission  effectiveness  resulting  from  failures  due  to  reliability  and/or  enemy  air  defense 
actions,  flic  electronics  are  described  by  the  subsystem  networks  required  for  each  mode  of 
operation  for  each  offensive/defensive  subfunction.  The  status  of  the  electronics  is  measured 
in  terms  of  the  probability  of  each  systems  state,  where  each  state  represents  a specific  com- 
bination of  operational  and  non-operational  subsystems.  Target  kill  probabilities  and  mode 
operational  probabilities  tor  each  electronic  system  state  are  computed  and  output  as 
MISDHM  steps  from  one  mission  event  to  the  next  in  chronological  sequence.  This  consti- 
tutes the  electronics  mode  of  operation  of  MISDHM.  In  the  vehicle  mode,  there  is  but  one 
subsystem,  namely  the  vehicle.  The  program  output  consists  of  probabilities  of  three  vehicle 
modes:  normal  flight,  abort,  and  down. 

A typical  application  would  involve  both  the  electronics  and  the  vehicle;  the  entire 
program  would  be  run  separately  for  the  electronics  and  vehicle,  utilizing  the  different 
MISDHM  operating  modes  for  each.  Each  is  thus  flown  through  the  mission.  Results  of  the 
MISDHM  run  tor  the  two  cases  could  subsequently  be  combined  by  the  user  in  a higher 
order  analysis.  As  an  example,  the  expected  target  kill  probability  could  be  computed  as  the 
product  of  probability  of  arrival  (obtained  in  the  vehicle  mode)  and  expected  target  kill 
probability  (obtained  in  the  electronics  mode). 
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Execution  of  the  MISDEM  simulation,  which  consists  of  two  parts  (Program  I and 
Program  2),  required  different  interpretation  and  types  of  data  for  the  electronics  and 
vehicle  modes.  The  required  input  data  are  somewhat  different,  even  though  in  some  cases 
the  same  variable  names  are  used. 


MATHEMATICAL  MODEL 


THEORETICAL  OVERVIEW 

Basically,  the  function  of  the  model  is  to  propagate  the  probability  distribution  of  the 
states  of  a system.  This  is  done  in  a recursive  (step-wise)  fashion.  The  recursion  formula  for 
each  state  is  the  Law  of  Total  Probability: 


K 

Pj,n  “ *Pj/i,n)x<Pi.n^  * ' ' 

where 

Pj  n = probability  associated  with  the  state  of  interest  at  event  n 

Pj  n = probability  associated  with  one  of  K mutually  exclusive  states  at  event 
n-1  (preceding) 

Pj^j  n = conditional  probability  of  state  j at  event  n.  given  a prior  system  state  i 
In  vector  notation: 


>1" 

pl/l  pl/2  *•*  Pl/K 

V 

P: 

• 

p: 

Pj.n 

• 

• 

~ 

• 

• 

• 

• 

• 

• 

PK 

_PK/1  pK/2  •••  PK/K 

_pK_ 

(pj/i.nl  pi,n-l 


where 

Pj  n = vector  containing  probabilities  of  all  states  at  the  nth  event 

(Pj/i  n)  = transition  matrix  for  the  nth  event,  a matrix  of  transition  probabilities 
(conditional  probabilities) 

Pj  n-i  = probability  vector  at  event  n-1 
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Once  an  initial  probability  vector  is  defined,  the  probability  vector  at  any  later  event  is 
simply  determined  by  successive  multiplications  with  all  transition  matrices  in  the  time 
interval. 

The  major  portions  of  the  computer  simulation  are  for  the  calculation  of  the  transition 
probabilities  and  the  generation  of  the  states. 

There  are  two  different  sets  of  algorithms  in  MISDEM  for  transition  probabilities 
corresponding  to  two  definitions  of  the  system  state  vector  for  the  electronics  mode  and 
vehicle  mode,  respectively. 

The  main  characteristic  that  divides  system  components  into  each  mode  is  their 
longevity  after  being  damaged.  Electronics  mode  subsystems  include  those  components  that 
die  fast  (viz,  in  the  time  increment  before  the  next  event  can  happen),  while  the  vehicle 
mode  subsystems  are  those  that  die  slowly  after  being  damaged  (i.e.,  they  may  be  operable 
for  several  time  increments  before  failing).  It  is  possible  that  some  subsystems  can  be 
included  in  both  modes  depending  upon  the  level  of  kill  to  be  examined  by  the  user.  Thus, 
fuel  cells  can  be  an  electronics  mode  subsystem  when  catastrophic  KK  kills  are  considered 
and  could  be  included  as  a subsystem  in  the  vehicle  mode  when  enemy  damage  causes  a leak 
that  results  in  total  loss  of  fuel  several  time  increments  later.  The  user  must  define  his 
electronics  and  vehicle  subsystems  in  a manner  that  satisfies  the  situation  he  is  interested  in 
studying. 

The  common  motive  in  each  application  is  to  simulate  the  loss  of  significant  aircraft 
functions  as  a result  of  exposure  to  damage  possibilities,  described  in  probabilities  (discrete 
or  continuous)  associated  with  threat  encounters  (events)  during  a mission. 


ELECTRONICS  MODE 
General 

In  this  mode,  the  program  computes  an  average  value  of  the  system  payoff  variable 
(measure  of  system  effectiveness)  where  the  average  is  taken  over  all  system  states.  The 
measure  of  effectiveness  is  the  number  of  targets  killed.  The  average  target  kill  probability 
for  each  target  reflects  the  multiple  alternate  delivery  methods,  their  individual  effectiveness 
(Pk)  and  their  probabilities  of  use  (related  to  the  system  state  probabilities).  The 
probabilities  of  operation  with  various  combinations  of  subsystems  are  also  outputs  of  the 
model. 


The  expected  number  of  targets  killed,  given  N targets  assigned,  is: 
N 

eT(N)=  Pk  n 
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where  l\  „ is  the  expected  kill  probability  in  a single  event.  This  output  parameter  is  the 
primary  measure  of  mission  effectiveness.  The  expected  kill  probability  at  a given  event  is 

JCAP 

•V,,*  £ (■’> 

J = I 

where 

J = weapon  delivery  mode  sequence  number 


Pj  n = probability  of  the  Jth  mode  sequence  in  the  nth  event 


I’k  j = kill  probability  in  Jth  delivery  mode  sequence 
JCAP  - number  of  mode  sequences 

This  output  parameter  is  of  interest  as  an  intermediate  result.  The  probability  of  the  Jth 
mode  sequence  in  the  nth  event  is: 

Nnax 


P 


j.n  = 


E 

j=  i 


ri,J 


(4) 


where 


Pj  | = probability  of  state  j,  which  supports  mode  sequence  J 
Kmax  = number  of  states 

Modes  and  Equipment 

The  behavior  of  the  aircraft  system  at  each  event  is  dictated  by  the  function/ 
equipment  list  that  is  input  by  the  user.  The  function/equipment  list  defines  the  equipment 
(subsystems)  necessary  for  the  execution  of  each  function.  It  may  be  possible  to  carry  out  a 
single  function  with  several  different  combinations  of  subsystems.  Tor  example,  there  max 
be  three  ways  to  deliver  unguided  iron  bombs  onto  a target  by  using  different  combinations 
of  radar  information,  barometer  readings,  on  board  computer,  and  iron  sights  depending  on 
the  information  supplied  by  the  subsystems  that  are  working.  Each  combination  of 
subsystems  that  allows  a mission  function  to  be  performed  is  considered  a mode  of 
operation  (as  opposed  to  the  more  general  electronics  mode  and  vehicle  mode  applied  to 
ways  of  using  the  MISDEM  model).  Thus,  the  user  specifies  for  each  mission  function  the 
possible  modes  of  operation  and,  where  needed,  the  effectiveness  of  the  system  when 
operating  in  each  mode.  It  is  assumed  that  the  user  will  identify  the  best  mode  for  a 
function  as  the  one  with  the  highest  effectiveness  and  all  other  modes  available  for 
performing  that  function  will  be  of  lesser  effectiveness  (i.e.,  degraded  modes  of  operation). 
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Electronics  response  to  a single  event  generally  consists  of  a number  of  functions  (e.g.. 
detect  threat,  employ  ECM  (electronic  countermeasure),  etc),  and  for  every  function,  one  of 
several  modes  is  employed  (e.g,,  jamming,  decoy,  chaff,  etc).  The  selection  of  the  mode  is 
determined  by  the  system  state.  Each  response  to  an  event  consists  of  a mode  sequence  (one 
mode  from  each  function).  The  ensemble  of  all  possible  mode  sequences  is  generated  at  each 
event,  and  a probability  is  computed  for  each,  based  upon  system  state  probabilities  as 
shown  in  Equation  (4). 

System  States 


The  system  state  is  defined  as  a viability  vector  which  is  a list  of  the  subsystems  that 
are  viable  (i.e.,  capable  of  being  used,  if  turned  on)  and  those  that  are  not.  For  example,  a 
simple  system  might  have  just  two  subsystems  (pilot,  engine);  then  at  zero  (takeoff)  time 
both  subsystems  are  working  (i.e,.  viable)  and  the  initial  system  state  would  be  defined  b\ 
the  viability  vector  (1,1).  The  subsystems  are  identified  by  their  position  in  the  vector 
according  to  the  ordering  (pilot,  engine).  If  the  pilot  is  lost,  the  system  state  is  defined  by 
the  vector  (0,1).  In  all,  four  vector  configurations  (system  states)  are  possible  for  this 
system:  (1,1),  (1,0),  (0,1),  and  (0,0).  Inherent  in  the  input  function/equipment  list  is  a 
priority  which  ensures  that  each  viability  vector  gets  associated  with  the  highest  priority 
mode  among  those  whose  equipment  requirements  are  satisfied.  A probability  for  every 
system  state  is  computed  based  upon  subsystem  survival  probabilities  and  the  mathematical 
rules  tor  combing  them.  The  probability  assigned  to  a mode  sequence  is  simply  the  sum  of 
the  probabilities  of  system  states  associated  with  that  particular  mode  sequence.  Since 
exactly  one  mode  sequence  can  occur  at  each  event,  the  sum  of  the  mode  sequence 
probabilities  is  unity  (except  as  the  user  may  sacrifice  some  accuracy  for  computer  running 
time  reduction). 

The  computation  of  the  system  state  probabilities  constitutes  the  basic  function  of  the 
program.  They  are  computed  at  each  event  in  a recursive  fashion  (see  Figure  3)  utilizing  the 
state  probabilities  from  the  previous  event,  and  computing  the  new  ones  based  upon 
conditional  probabilities  of  transition  ("TRANS”)  related  to  reliability  failure  rates  and 
subsystem  Ps  (probabilities  of  survival)  associated  with  threat  warhead  encounters,  for  a 
given  value  of  CEP  (circular  error  probable)  of  guidance. 

In  the  program,  the  viability  vector  is  represented  as  a binary  number. 

[0  I 100  I l| 

where 

“I”  denotes  a viable  subsystem 

"0"  denotes  a subsystem  that  is  not  viable 

“M”  equals  number  of  subsystems 

Viability  is  defined  (implicitly  in  the  program  logic)  as  "capable  of  supporting  a system 
function,  if  turned  on”. 
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is  this 

SUBSYSTI m 
. OOWN  ' _ 


MUHIPLY  RESULT  BY  PREVIOUS  RESULT 


MORE  SUBSYSTI  MS 


NO  MORI  SUBSYSTEMS 


MORE  STATES  I 


NO  MORE  STATES  I 


MORE  STATES  J 


NO  MORE  STATES  J 


Figure  3.  Conventional  Damage  Transition  Algorithm 


1 1 I I I I I 1 1 (the  perfect  state  for  a 7 component  system) 


in  which  c, 
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where  the  first  state  (at  the  top)  is  the  perfect  state  and  states  of  2 through  k represent 
system  operation  with  various  components  down  or  off.  The  maximum  value  of  K for  this 
system  would  be  27. 

The  only  way  a subsystem  becomes  not  viable  is  that  it  is  damaged  to  failure  by  natural 
(reliability)  or  induced  (enemy  warhead)  damage. 

Transition  Algorithms 

In  a perfect  system  with  no  failures,  the  transition  matrix  is  the  identity  matrix 
(“ones’*  on  the  main  diagonal  and  “zeros”  elsewhere)  and  the  state  probability  vector  does 
not  change  throughout  the  mission.  (This  can  be  demonstrated  easily  bv  applying 
Equation  I .)  The  method  for  the  computation  of  the  transition  probability  is  an  application 
of  combinatorial  analysis  to  the  probabilities  of  failures  of  the  several  subsystems. 

UNCORRELATE1)  EAI  LURES.  When  the  subsystem  failures  are  uncorrelated,  the 
state  transition  probability  is  the  product  of  the  subsystem  failure  probabilities.  As  an 
example,  consider  the  following  transition  and  its  probability: 


Example  of  Uncorrelated  Transition  Probability  Generation 


Failure  probabilities 

M )M 

Subsystem  ordinal  number  k 

1 2 3 4 5 6 7 

System  state  i 

10  0 1 10  1 

System  state  j 

1 0 0 0 1 0 0 

Transition  probability  Pj/j 

Pj  •M*q4*p5*l-q7 

where 

Pk  = l-qk 

(Once  a subsystem  has  failed,  the  probability  of  recovery  is  zero  because  no  repair  is 
assumed  during  the  mission.) 

This  process  is  represented  by  an  algorithm  that  is  documented  by  Air  f orce  Systems 
Command  (footnote  1 ) (the  logic  flow  is  shown  in  Figure  3). 

The  transition  algorithm  for  this  case  has  subsystem  reliability  (in  the  transition)  as  a 
multiplying  factor.  Thus.  Pk  is  composed  of  two  factors:  dk  (reliability)  and  vk  (survival  of 
threat  damage). 
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RELIABILITY.  File  factor  dk,  representing  subsystem  reliability  in  the  transition,  is 
obtained  from 

dk  = exp  (-At  * MTBF)  = I - At/MTBF  + At2/2MTBI  2 (5) 

where 

At  is  the  elapsed  subsystem  time  between  the  last  event  and  the  current  event 

MTBF  is  the  mean  time  between  failures 

KILT.  PROBABILITY.  The  factor  vk  represents  subsystem  survivability  in  the 
transition,  which  is  the  same  as  the  subsystem  conditional  survival  of  an  encounter  occurring 
in  that  interval  starting  after  tn.|+  and  ending  prior  to  tn-.  Damage  that  develops  completely 
during  the  time  interval  At  is  referred  to  as  quick  damage.  This  is  characteristic  of  electronic 
damage.  [ Damage  that  does  not  develop  completely  until  after  subsequent  damage  exposure 
is  termed  slow  threat  damage  and  is  characteristic  of  vehicle  damage  (crack  propagation, 
fluid  leakage,  and  fire  propagation).  This  type  of  damage  mechanism  is  handled  by  the 
vehicle  mode  of  MISDFM.  | 

niRI  AT  WEAPON  GUIDANCE  FRRORS.  The  kill  computations  in  the  program 
provide  tor  variations  in  threat  weapon  lethality  caused  by  aircraft  countermeasures  by 
having  the  user  input  the  degraded  weapon  effectiveness. 

As  shown  in  Figure  4,  the  modes  in  a defensive  event  are  identified  with  input  values  of 
threat  CEP.  These  CLP  are  used  to  assign  probabilities  to  the  data  base  of  component 
survival  probabilities  generated  by  the  one-on-one  effectiveness  models. 

The  data  base  is  generated  externally  to  MISDFM  by  selecting  a simulation  that  can 
handle  a representative  array  of  fragment  or  projectile  sources.  For  the  proximity-ftized 
missile  warhead,  a typical  approach  for  modelling  missile-target  encounter  conditions  would 
be  to  select  an  array  of  offset  trajectories  (see  Figure  5)  centered  about  a representative 
mean  trajectory.  Each  of  these  trajectories  is  then  analyzed  by  an  end-game  fuzing  model  to 
determine  the  exact  burst  point  location  in  aircraft  coordinates. 

Figure  6 shows  the  offset  distribution  for  a CEP  of  100  feet,  and  the  probability 
associated  with  a 100-  to  200-foot  offset  that  is  applied  to  the  calculation  of  system  state 
transition  probability  in  the  “TRANS”  block  of  each  event  (Figure  3). 

The  probability  of  aircraft  damage  is  computed  by  integrating  (over  the  limits  of  each 
zone)  the  bivariate  normal  distribution  with  zero  mean.  The  standard  deviation  (sigma)  is 
derived  from  the  threat  weapon  CEP  as  follows. 

a * CEP/I.1 78  (0) 

LI 


1 1 
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The  result  of  the  integration  of  the  bivariate  normal  distribution  over  limits  Rkand  Rk+j 
<Rk+|  > Rk)  is: 

Pmiss(k)  = exp  <Rk-/2  o-)- exp  (Rk+|  -/2  a-)  (7) 


The  transition  probability  associated  with  such  an  event  is  then  computed  as: 


.I/' 


L 


Pj/i  (given  burst  point  k)  x P (burst  point  k) 


(8) 


over  all  L burst  points  described  in  Figure  4.  One  factor  in  P(burst  point  1 ) is  Pnljss  which  is 
applied  to  all  bursts  in  a given  zone.  The  other  factor  is  the  probability  of  the  burst 
occurring  at  one  of  the  L points  in  that  zone.  This  probability  is  just  1 + L,  reflecting  an 
assumed  circular  symmetric  trajectory  distribution. 

C'ORRFLATKD  FAILURES.  Another  major  type  of  damage  mechanism  that  is  treated 
by  MISDFM  is  associated  with  nuclear  weapons  (neutron  and  gamma  radiation,  blast  and 
thermal  effects).  In  this  case,  the  mutual  shielding  of  one  component  by  another  is  assumed 
negligible  and  the  predominant  correlation  of  subsystem  failure  is:  if  component  “A”  is 
killed,  all  subsystems  with  lesser  nuclear  hardness  are  also  killed.  The  transition  probability 
algorithm  makes  use  of  a lethal  radius  as  a measure  of  vulnerability  of  each  subsystem,  and 
uses  the  probability  distribution  of  weapon  miss  distance  normally  distributed  with  zero 
mean  to  deduce  the  probability  of  a given  set  of  subsystem  failures.  Figure  7 shows  the  miss 
distance  density  function  in  r (radius)  (the  Rayleigh  marginal  density).  The  subsystem  lethal 
radii  are  indicated  on  the  abscissa:  the  subsystem  ordinal  numbers  have  been  assigned  in 
decreasing  order  of  lethal  radius.  The  probability  of  kill  of  exactly  the  first  k subsystems  is 
the  probability  that  the  miss  distance  exceeds  rk+j  but  does  not  exceed  rk  (equal  to  the 
integral  of  the  miss  distance  density  function  with  lower  limit  rk+|  and  upper  limit  rk. 
designated  qk).  An  implicit  assumption  here  is  that  the  warhead  damage  effect  is  a 
monotonic  decreasing  function  of  distance  from  the  warhead.  The  following  is  an  example 
of  the  computation  of  Pj/j  in  a case  where  the  least  hardened  subsystem  fails. 

Failure  probabilities: 


1.  Natural  (reliability) 

2.  Unnatural  (threat) 
Subsystem  ordinal  number  k = 
System  state  i 

System  state  j 

Transition  probability  Pj/j  = 


q7q6(l5‘4tl3<l2tM 

qViociVi4‘iViViV'() 
7 6 5 4 3 2 1 
I 0 0 I 1 I I 
I 0 0 0 I 10 

P7I  • '•  <14P3P2<<1  i q'o 


: 

l 
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• XXNT  NORMAL  ERROR  DISTRIBUTION 

• ZERO  MEAN 

• CEP  - 100  FT 


The  rationale  for  the  final  step  is  as  follows: 

1.  Since  subsystem  2 survived  the  unnatural  (threat)  hazard,  the  probabilities  of  all 
subsystems  of  greater  hardness  (3  through  7)  surviving  are  determined  by  reliability  only. 
Therefore,  in  the  example  P7  = I-M7,  subsystems  6 and  5 are  not  viable  and  are  assigned  a 
transition  probability  of  1.  Subsystem  4 has  failed  naturally ; consequently,  the  transition 
probability  is  noted  as  one  of  failure.  Subsystems  3 and  2 are  treated  the  same  as 
subsystem  7. 

2.  Subsystem  I (the  least  hard)  has  a failure  probability  composed  of  the  sum  of  the 
probabilities  of  two  mutually  exclusive  events. 

a.  The  probability  of  threat  induced  kill,  expressed  as  q'|  (the  probability  of  the 
miss  distance  being  within  the  lethal  radii  for  subsystem  I but  more  than  the  radius  for 
subsystem  2)  times  the  probability  of  the  natural  failure  or  survival  of  the  subsystem, 
qj  + P|.  This  may  be  expressed  as  q’  • (qj  + P|)  but  since  Pj  = I - qj.  this  may  be 
shortened  to  q’j(  1 ) or  simply  q'j . 
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b.  The  probability  that  the  subsystem  failed  from  natural  causes,  qj;  and  the 
miss  distance  was  greater  than  the  greatest  lethal  radius,  qg.  Since  these  are  indepen- 
dent events  their  probability  must  be  expressed  as  (qiqg). 

As  a more  complicated  example,  consider  the  following: 

Subsystem  ordinal  number  k=  7654  3 21 

System  state  i 1 I I I I 0 I 

System  state  j I I 0 I 0 0 0 

Transition  probability  Pj/j  = 1 (1I|(i1())+mV 

-Hpl-Ki^)  +0)  +0)  +0)+0) 

In  this  form,  a computer  algorithm  for  Pj/j  starts  to  become  evident.  The  qg  damage 
term  is  always  present,  successive  reliability  factors  are  incorporated  until  all  subsystems  are 
accounted  for  and  successive  damage  terms  tq'j.)  are  added  until  the  leading  zeros  in  the  j 
state  are  exhausted. 

The  algorithm  is  presented  in  Figure  8.  A flag  is  set  to  I after  the  leading  zeros  are 
exhausted,  in  order  to  change  the  algorithm  to  eliminate  the  additive  q\  In  addition,  if  the 
last  event  was  not  nuclear,  the  algorithm  reverts  to  a reliability-only  transition. 

The  determination  of  the  subsystem  lethal  radii  is  based  on  two  factors,  illustrated  in 
Figure  ( I)  weapon  damage  effect,  a characteristic  of  the  warhead  and  geometry  alone, 
measured  in  pounds  per  square  inch  (for  blast),  calories  per  square  centimeter  (for  thermal), 
rads  per  second  (for  gamma),  and  neutrons  per  square  centimeter  (for  neutrons);  and 
(2)  subsystem  damage  threshold,  a characteristic  of  the  subsystem  alone,  measured  in  the 
same  units  as  weapon  damage  effect. 

The  definition  of  subsystem  lethal  radius  is  the  maximum  value  of  the  lethal  radii 
computed  for  all  damage  mechanisms,  and  differs  in  general  from  one  subsystem  to  the  next 
due  to  differences  in  subsystem  damage  thresholds.  The  algorithm  which  selects  this  unique 
lethal  radius  for  each  and  every  subsystem  is  diagrammed  in  Figure  10.  The  lethal  radius  is 
updated  on  every  Mth  pass  through  the  last  block  in  the  figure,  where  M is  the  number  of 
subsystems. 
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\DOwN,  OR  yt 

. UP?  / 


(Subsystem  killed; 

AFFIX  RU  I ABILITY  FACTOR) 


W«N 


SUBSYSTEM  RESURRECTED 
IMPOSSIBLE  event 
SET  P , - 0 


SET  FLAG  * I 
SUBSYSTEM  SURVIVES) 


jjr  i rriorX 
S.RSrSTfM  DOWN 


(Subsystem  still 

KILLED) 


IS  FLAG  - O' 


US]  ADD  NUCLEAR  SUM  i X 


NUCLEAR  KILLS  STILL  POSSIBLE 
ADO  NUCLEAR  SUFFIX 


(NUCLEAR  KILL 
NOT  POSSIBLE) 


COST INUE 


LAST  SUBSYSTEM 


REMAINING  SUBSYSTEMS  ARE  AFFECTED  ONLY  BY  RELIABILITY 


* P ) ALL  STATE  PAIRS 


Figure  8.  Nuclear  Transition  Algorithm 
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Figure  10.  Compute  Lethal  Radius. 


The  next  step  in  reordering  subsystems  is  illustrated  in  Figure  1 1.  The  output  of  the 
algorithm  is  a lethal  radius  RM(k),  and  an  order  number  IZ(k)  for  each  subsystem  number 
(k)  assigned  in  order  of  decreasing  lethal  radius.  The  program  selects  the  first  subsystem, 
then  tests  the  next  subsystem  to  determine  if  it  has  a larger  lethal  radius.  If  it  does,  the 
lethal  radius  and  the  order  number  of  the  first  subsystem  are  assigned  to  argument  2,  and 
the  lethal  radius  and  the  order  number  of  second  subsystem  are  assigned  to  argument  1. 
Therefore,  the  order  number  of  subsystem  2 is  now  1,  because  it  is  larger;  and  the  order 
number  of  subsystem  1 is  now  2.  Having  selected  the  largest  lethal  radius  for  argument  1, 
the  computer  then  tests  the  third,  fourth,  etc.,  (to  the  end)  against  the  first;  and  if  any  arc 
larger  than  the  first,  additional  switches  take  place.  This  process  of  pair-wise  comparisons  is 
repeated  (using  the  outer  do-loop)  to  load  RM  (2)  and  1Z  (2)  with  the  subsystem  having  the 
largest  lethal  radius  among  those  other  than  RM  (1)  (which  is  never  tested  again),  and  is 
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llwrc (ore  second  largest.  I he  process  is  repeated  until  each  position  in  RM  and  1/  arrays  arc 
filled-  H'1'  array  contains  the  vulnerability  order  number  for  each  original  subsystem 
number.  I bus.  the  subsystem  viabilities  are  reordered  by  means  of  the  I/,  array  before  the 
transition  probability  is  computed. 


1 he  weapon  miss  distance  density  function  is  a product  of  probabilities  of  two  random 
variables  (Pr/t|  x Pt|). 

1.  Pr/()Ut = conditional  miss  density  (conditioned  on  the  event  that  the  warhead  was 
delivered  to  the  target  neighborhood)  being  a function  of  the  weapon  guidance  accuracy 
(including  deceptive  countermeasure  effects)  in  terms  of  OIP 

*>r/d,r)  = ~^c'\p  (-r-/2tf -)  - ( l .4r/(TP-)exp(-.7r“/0  P-) 

The  probability  that  the  miss  exceeds  r^  is: 

CO  00 

r/d(r)  l*r  = j ~^exp  (-r~/2o‘)  dr 
rlc  rk° 

Rearranging  the  limits  to  introduce  a negative  sign  under  the  integral  to  put  the  integrand  in 
the  standard  form  of  eudu  produces: 

fk 

^r/d*1'  /-(r /a - ) exp (-r“/2o~)  dr  = exp  (-r~/2o') 

a? 

2.  IM  = probability  ol  delivery,  being  a product  of  probabilities  of  (assumed) 
independent  events. 

a.  Probability  of  weapon  launch  (a  function  of  detection  and  launch  delays,  as 
affected  by  denial  and  confusion  countermeasures) 

b.  Probability  of  in-flight  reliability  (guidance,  control  and  fining,  independent 
of  countermeasures) 

c.  Probability  ot  no  catastrophic  guidance  denial  due  to  destruction  of  a portion 
of  the  guidance  system. 

The  expected  number  of  rounds  delivered  is  SPd.  where  S is  the  maximum  number  of 
rounds  possible.  The  expected  number  is  called  N. 

To  reduce  the  number  of  program  iterations  required  for  simulating  a mission,  the 
effects  of  several  (N)  weapons  encountered  in  a single  threat  encounter  can  be  combined 
into  a single  effect,  assuming  the  single  weapon  miss  density  functions  are  identical. 
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Figure  1 1.  Assign  Vulnerability  Index  Numbers 


The  derivation  is  as  follows:  N shots  are  placed  in  such  a way  that  all  N lie  outside  i^*  | 
(called  region  AUB)  and  at  least  one  lies  inside  i|^  (called  tvgton  A)  Under  these  conditions, 
(only  I subsystems  numbered  k ot  less  are  killed  I'he  probability  of  this  occurring  is  I*  (all  in 
AUB  and  at  least  one  in  A) 

■ I*  (all  in  AUB)  \ I’ (at  least  one  in  A jtiven  all  in  AUB) 

■ I’ (all  in  AUB)  \ 1 1 • I*  (none  in  A given  all  in  AUB)| 

* P (all  in  AUB)  s.  1 1 - 1*  (all  in  B (liven  all  in  Al'BII 

- P (all  in  AUB)  \ | I • P(all  in  B and  all  in  AUB)  P (all  in  AUB)| 

- P (all  in  AUB)  - P (all  in  B and  all  in  AUB) 

- P (all  in  AUB)  - P (all  in  B) 
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But 

SPd  , 

P (all  in  AUB)  = IIP^  (ith  round  in  AUB)  = exp  (-SPdr£+|  /2  a-) 

and 

SPd 

P (all  in  B)  = IIP  (ith  round  in  B)  = exp  (-SPdrj:  /2 a~) 
i ”*  1 

Thus 

q'k.i  ® exP  (-SPdrfc+|/2  Oj~)  - exp  (SPdrg  / 2oj~) 

(9) 

' 

I 


The  approximation  lies  in  the  use  of  the  expectation  of  the  number  of  rounds. 

The  index  i denotes  the  dependence  of  these  parameters  on  the  state  of  the  system 
trom  which  a transition  is  to  be  made.  Of  course,  the  significant  states  of  the  system  in  this 
ease  are  those  with  different  defensive  avionics  subsystems  viable.  Figure  12  summarizes  the 
elements  of  the  subsystem  nuclear  kill  probability  computation. 

State  Generation 


The  total  probability  associated  with  the  system  state  is  1.0  (i.e.,  the  system  must  exist 
in  some  state)  consequently: 

K n 

£ P;  „ = 1 tor  all  n 

j=  1 

where 


K is  the  total  number  of  mutually  exclusive  and  exhaustive  system  states  of  interest 
(i.e.,  whose  probability  is  of  interest) 

n is  the  event  number 

Pj  is  the  probability  associated  with  the  state  j 

Therefore,  a requirement  for  defining  states  in  the  model  is  that  they  be  mutually 
exclusive  and  exhaustive. 

In  the  electronics  mode,  each  state  is  a binary  M-vector  representing  subsystem 
viabilities.  Exclusivity  is  apparent  because  each  state  is  a different  binary  number.  To  ensure 
an  exhaustive  set  be  identified,  it  is  only  necessary  to  generate  all  possible  permutations  of 
ones  and  zeros  in  the  (ordered)  array.  This  is  done  by  starting  with  all  ones,  then 
introducing  a zero,  moving  it  from  one  end  to  the  other  in  the  array,  then  introducing  a 
second  leading  zero,  letting  the  other  run,  then  a third  zero,  etc.,  until  M zeros  have  been 
used. 


For  each  state  generated,  two  data  items  are  developed  in  the  program:  (l)its 
probability  and  (2)  the  mode  sequence  that  the  system  would  exhibit,  given  that  the  system 
had  that  state. 
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M«k1c  Sequence  Identification 


A sample  system  configuration  is  presented  in  Table  I.  The  system,  as  defined  for  an 
electronics  mode  MISDEM  run,  is  composed  of  thirteen  subsystems  or  components.  Some 
of  the  components  listed  may  not  always  be  subject  to  quick  killing  (vi/.  engines,  pilot,  fuel 
systems  and  hydraulic  components).  However,  assume  for  this  illustration  that  we  are 
interested  only  in  assessing  the  systems  susceptibility  to  KK-kill  so  that  the  questionable 
components  do  become  considered  in  the  electronics  mode. 


Table  2 is  the  input  which  defines,  for  the  four  functions  to  be  performed  for  event  I. 
the  modes  and  equipments  used  by  the  system.  Four  functions  are  possible  in  the  example: 
flight,  detection,  evasion,  and  countermeasures.  For  the  flight  function,  one  of  three  modes 
must  be  performed,  in  order  of  priority: 

1.  The  normal  flight  control  system  must  operate  (the  preferred  mode). 

2.  The  backup  tlight  control  system  must  be  activated  (the  second-best  mode). 

3.  The  aircraft  is  disabled. 
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For  the  normal  flight  control  system  mode,  subsystems  I,  2.  and  4 through  12  must  be 
viable:  however,  if  this  condition  is  not  met,  the  next  mode  is  examined  by  the  computer, 
for  which  only  the  electromechanical  subsystems  are  required.  The  first  two  number 
columns  arc  ordinal  numbers  for  the  indicated  function  and  mode.  The  second  two  numbers 
are  employed  by  the  user  to  specify  which  of  the  other  functions  and  modes  are  to  be 
performed,  given  the  success  ( i.e. . viability)  of  the  first  function  and  mode.  The  program 
perforr  s a cataloging  operation  in  which  all  significant  candidate  systems  states  (one/zero 
configurations)  are  generated  one  at  a time  and  each  is  tested  against  the  mode  structure  to 
identify  the  best  mode  sequence  ascribable  to  each  and  every  state.  In  series  with  the 
viability  test  is  the  mission  descriptor  test.  Each  function/mode  may  require  certain 
environmental  or  arbitrary  program  controls  placed  upon  it  by  the  user.  This  requirement,  a 
I or  F (True  or  False),  shows  up  in  the  last  column.  If,  as  in  function  2,  mode  2.  passive 
(e  g.,  infrared)  surveillance,  a clear  air  mass  is  required  (indicated  by  a T in  the  last  column), 
this  function  cannot  be  performed  in  event  1,  since  the  corresponding  value  in  the  event 
description  (F  in  the  last  column)  indicates  that  the  weather  does  not  cooperate.  Conse- 
quently, the  computer  program  is  directed  to  another  function/mode  by  means  of  the  last 
two  number  columns  to  the  left  of  the  function/mode  descriptions. 

Table  3 presents  the  first  two  mode  sequences  deduced  for  event  I by  the  computer. 
The  first  (best)  mode  sequence  is  generated  by  the  perfect  system  state  (all  ones  in  the  state 
vector  (not  shown)).  All  distinct  flow  paths  are  normally  generated.  As  a check  on  the  user's 
logic,  the  corresponding  required  equipments  are  also  output.  The  second  mode  sequence  is 
degraded  and  results  from  a loss  of  subsystems  (the  jammer  being  out). 


Table  I.  Sample  System  Configuration. 


Equipment 

MTBF 

1 Left  engine 

1000.00 

2 Right  engine 

1000.00 

3 Sensors/  ECM 

1000.00 

4 Comm/nav/eomputer/instr/displays 

1000.00 

5 Pilot 

1000.00 

6 Fuselage  fuel  system 

1000.00 

7 Wing  fuel  system 

1000.00 

8 Left  master  cylinder 

1000.00 

9 Right  master  cylinder 

1000.00 

10  Rudder  actuator 

1000.00 

1 1 Horizontal  stabilizer  actuator 

1000.00 

1 2 Aft  hydraulic  power  pack 

1000.00 

13  Electro-mechanical  backup  night  control  system 

1000.00 

u 

I 
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Table  2.  Mode  Sequence  Logic  For  Kvent  1. 


Function/mode 

Equipment 

description 

Mission 

descriptor 

1 

1 

1 

Flight  functions 

0000000000000 

0 

0 

F 

1 

1 

Normal  (light  control  system 

1 101 1 1 1 1 1 1 1 10 

0 

0 

F 

1 

3 

1 

Back-up  (light  control  system 

1 101 1 1 1001 101 

0 

0 

F 

1 

4 

5 

5 

Aircraft  down 

0000000000000 

0 

0 

F 

0 

99 

0 

0 

0000000000000 

0 

0 

F 

1 

*> 

Threat  detection 

0000000000000 

0 

0 

F 

3 

1 

Passive  surveillance 

00 111  00000000 

2 

4 

T 

3 

3 

1 

Visual  mode 

0000100000000 

4 

T 

*> 

4 

4 

1 

No  detection  and  no  maneuver 

0000000000000 

0 

0 

F 

0 

99 

0 

0 

0000000000000 

0 

0 

F 

3 

1 

3 

2 

Evasive  maneuver 

0000000000000 

0 

0 

F 

3 

2 

4 

i 

Normal  mode 

1 101 1 1 1 1 1 1 1 10 

0 

0 

F 

3 

3 

4 

i 

Limited  g-mode 

110111 1001 101 

0 

0 

F 

0 

99 

0 

0 

0000000000000 

0 

0 

F 

4 

1 

4 

Countermeasures 

0000000000000 

0 

0 

F 

4 

*> 

5 

5 

Jam  mode 

001 1000000000 

0 

0 

F 

4 

3 

5 

5 

Jammer  out 

0000000000000 

0 

0 

F 

0 

99 

0 

0 

0000000000000 

0 

0 

F 

999 

0 

0 

0 

0000000000000 

0 

0 

F 

NOTE:  Evant  description — avant  1 is  defensive;  event  occurred  0.19  hours  after  takeoff;  event  description 
is  SAM  miesile  encounter  F. 


Each  mode  sequence  in  Table  3 is  tested  and,  if  unique,  is  assigned  a unique  number. 
This  is  done  in  order  that,  later  in  the  program,  capabilities  may  be  associated  with  these 
sequences  by  number.  The  process  of  identifying  and  numbering  the  mode  sequences  is 
repeated  for  each  state  supplied  by  the  state  generator.  The  process  is  shown  schematically 
in  Figures  13  and  14. 

Figure  13  shows  how  the  state  is  tested  against  mission  and  subsystem  requirements  of 
each  subfunction,  starting  with  the  first  mode  in  each  case.  When  satisfaction  is  found  in  a 
mode  (for  each  subsystem),  its  required  subsystems  are  loaded  into  a “subsystems  used” 
array,  and  the  current  subfunction  and  mode  indices  are  recorded  (for  the  current  state)  in  a 
mode  sequence  array. 

Figure  14  shows  how  each  new  mode  sequence  array  is  tested  against  its  predecessors 
and,  if  new,  a new  mode-sequence  number  is  assigned  to  the  current  state.  If  it  is  duplicated 
by  a predecessor  (i.e.,  is  not  unique)  the  old  mode  sequence  number  is  ascribed  to  the 
current  state.  In  either  case,  the  appropriate  mode  and  subsystem  lists  are  associated  with 
the  number  as  indicated  (Table  3). 
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Tabic  3.  Sample  Mode  Sequences  Identified. 


Mode  sequence 

Subfunctional  flow 

1 

Flight  functions 

Normal  flight  control  system 

Threat  detection 

L 

* 

No  detection  and  no  maneuver 

Countermeasures 

Jam  mode 

i 

Subsystems  used 

Left  engine 

Right  engine 

Sensors/ ECM 

Comm /nav /com  puter/instr/displays 

Pilot 

Fuselage  fuel  system 

i 1 

Wing  fuselage 

Left  master  cylinder 

Right  master  cylinder 

1 

Rudder  actuator 

Horizontal  stabilizer  actuator 

Aft  hydraulic  power  pack 

*> 

Flight  functions 

i 

Normal  flight  control  system 

Threat  detection 

No  detection  and  no  maneuver 

c 

Countermeasures 

Jammer  out 

1 | 
i 

J 

Etc. 

VEHICLE  MODE 

Outputs.  States,  and  Mode  Logic 

In  this  mode,  the  program  computes  the  probability  of  arrival  of  the  vehicle  at  each  of 
a succession  of  events  (threat  encounters,  targets,  home  base,  etc.).  The  probability  of 
arrival  is  defined  as  the  probability  that  the  aircraft  did  not  go  somewhere  else  (as  a result  of 
an  abort)  nor  go  down  prior  to  reaching  that  event. 

0 • " ' 

u 
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Figure  13.  Define  Mode  Sequence  and  Subsystems  Used. 
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DEFINE  (NEXT)  MODE  SEQUENCE,  TENTATIVE  PATH 
NO.,  AND  SUBSYSTEM  LIST  FOR  CURRENT  STATE 


DEFINE  (NEXT)  EARLIER  MODE  SEQUENCE 


IS  THIS  THE  > 
1ST  MODE  SEQUENCE? 


SELECT  TENTATIVE 
MODE  SEQUENCE  NO 


IS  CURRENT  FLOW 
PATH  IDENTICAL  TO 
EARLIER  MODE  SEQ? 


ASSIGN  CURRENT  STATE  NUMBER  TO 
THIS  MODE  SEQUENCE  NUMBER 


Figure  14.  Assign  Mode  Sequence  Numbers. 


JTCG/AS-76-S-004 


The  vehicle  mode  is  primarily  used  for  the  study  of  slow  damage  mechanisms  and  the 
determination  ol  the  aircraft  aborts  that  return  to  base  for  repairs.  The  mission  time/event 
series  needs  to  be  reorganized  at  this  point  into  relative  mission  time  units.  Events  which 
will  influence  the  probability  of  arrival  of  the  vehicle  are  the  defensive  events  (viz.  the 
exposure  ot  the  aircratt  to  AAA,  SAM,  or  air-to-air  threats).  The  number  of  mission  time 
units  in  a mission  corresponds  to  the  number  of  defensive  events  (e.g.,  a mission  that 
contains  four  defensive  events  has  four  mission  time  units).  Each  mission  time  unit  contains 
a constant  amount  of  clock-time.  As  an  example,  a mission  that  contains  four  defensive 
events,  at  mission  times  ot  0.1, 0.26,  0.43,  and  0.62  hours,  and  is  1.2  hours  in  length  would 
have  four  mission  time  units  each  0.3  hours  in  length  (1.2/4  = 0.3).  Therefore,  the  actual 
defensive  events  would  be  lumped  into  appropriate  mission  time  unit  events.  These  events 
would  be  separated  by  a constant  time  difference.  At.  Due  to  the  coding  of  the  program,  the 
analyst  is  restricted  to  using  either  four,  eight,  or  sixteen  of  these  reorganized  events.  If  the 
mission  has  a number  of  defensive  events  that  exceed  four  or  eight,  the  next  larger  number 
ot  mission  time  units  should  be  utilized.  For  example,  for  five  defensive  events,  eight 
mission  time  units  would  be  used  (three  of  the  “empty”  events  would  be  assigned  zero 
effectiveness). 

Figure  15  illustrates  a mission  with  eight  events,  with  two  of  them  possibly  leading  to 
mission  aborts.  Normally,  point  8 would  represent  recovery  at  an  airbase.  For  example, 
points  I through  3 and  5 through  7 might  be  threats,  and  point  4 might  be  the  target  and 
turnaround  point.  The  usual  calculation  of  arrival  probability  is  accomplished  for  the 
normal  path  by  a single  run  of  MISDFM.  Any  of  the  points  1 through  7 may  also  be  used  as 
starting  points  in  separate  MISDFM  calculations  of  arrival  probability  along  abort  routes.  By 
examining  the  output  for  the  initial  events  composing  the  mission,  the  user  may  want  to.  on 
the  basis  of  a pre-selected  threshold  value  of  the  probability  of  arrival,  declare  that  the 
mission  would  have  been  selected  to  result  in  the  normal  mission  being  aborted.  The 
time/event  series  designated  la.  lb,  etc.,  described  abort  routes  A and  B.  The  user  defines 
events  and  modes  of  equipment  functioning  for  these  time/event  abort  routes  in  a manner 
similar  to  that  utilized  for  the  initial  events  comprising  the  normal  mission  (i.e.,  event  2a  for 
abort  A might  be  attacking  an  alternate  target  because  the  aircraft  lost  its  radar  bombing 
capability,  needed  for  event  4,  at  event  3).  The  abort  route  simulation  differs  from  the 
normal  path  simulation  only  in  the  initial  values  of  the  aircraft  state  probabilities. 

The  logic  reflected  in  the  model  is  as  follows.  At  each  event  the  aircraft  may  be 
deemed  to  be  operating  in  one  of  three  modes  (normal,  abort,  or  down)  based  upon  threat 
damage  accumulated  (and  propagated)  prior  to  the  time  of  the  current  event.  Unless  it  is  in 
the  normal  mode,  it  is  not  normally  counted  as  arriving.  The  central  task  of  the  computer 
program  is  to  compute  the  arrival  probabilities  associated  with  the  normal  mode  at  each 
event.  The  probabilities  of  abort  and  down  modes  are  also  outputs  of  the  model.  The  model 
probabilities  are  dependent  on  the  probabilities  of  damage  subsequent  to  exposure.  Slow 
damage  mechanisms  are  assumed  to  operate,  so  that  the  effects  of  encounters  may  occur 
many  At  later.  The  time  convention  used  is  that  a defensive  event  is  the  time  at  which  the 
threat  exposure  is  initiated  (not  concluded),  and  is  the  time  from  which  damage  effects  are 
referenced. 
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Figure  15.  Schematic  Scenario  for  Vehicle  Mode,  Fight  Events. 
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For  each  event,  there  is  a cumulative  probability  distribution  of  possible  times-of-flight 
( t > remaining  relative  to  the  start  of  the  exposure.  These  times  are  expressed  in  mission  time 
units.  These  probabilities  are  sometimes  known  (empirically)  for  specific  times,  such  as 
5 seconds  (KK-kill).  30  seconds  (k-kill),  etc.,  and  may  be  used  to  construct  the  probability 
distribution  (a  required  input)  throughout  the  mission  duration.  Depending  upon  the  user 
(input)  mode  logic,  an  aircraft  may  be  classified  as  down  whenever  t — : 1 (t  is  measured  in 
mission  units  aiu!  is  always  an  integer). 

In  addition,  the  aircraft  may  be  classified  as  aborted  when  a potentially  mission- 
limiting  failure  is  detected.  The  time  that  this  detection  occurs  is  labeled  r and  is  called 
abort  detection  time.  It  is  distributed  analogously  to  t and  is  similarly  associated  with  a 
threat  encounter  and  measured  in  mission  units.  Figure  lo  illustrates  that  logic  for 
classifying  aircraft.  An  example  of  the  two  distributions  for  a single  encounter  is  presented 
in  Table  4.  The  probability  of  aircraft  failure  as  a function  of  mission  time  units  is  obtained 
by  cumulating  the  effects  of  the  defensive  events  as  the  mission  progresses.  The  abort 
probability  as  a function  of  mission  time  units  following  a defensive  event  is  supplied  by  the 
user.  Interpretation  of  the  abort  probabilities  is  as  follows:  considering  the  third  time  unit 
following  exposure,  the  probability  of  the  aircraft  failing  during  this  time  increment  due  to 
accumulated  damage  is  0.05,  and  the  probability  of  detecting  the  need  to  abort  the  mission 
in  the  same  time  unit  following  exposure,  is  0.4.  However,  the  damage  build  up  that  results 
in  a failure  time  t = 3 would  produce  a probability  of  0.3  that  the  crew  would  detect  the 
need  to  abort  in  the  first  mission  time  unit  following  exposure;  abort  detection  probability 
for  time  2 is  0.6,  and  by  time  4 the  crew  would  always  detect  the  need  to  abort  the  mission, 
since  the  aircraft  failed  in  three  time  units  following  exposure. 
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FOR  EVENT 


AIRCRAFT  DOWN 


AIRCRAFT  NORMAL 


AIRCRAFT  ABORT 


Figure  16.  Simple  F.xample  of  Mode  Logic  for  Vehicle  Case. 


As  shown  in  Table  4,  the  abort  detection  time  distribution  has  been  made  a function  of 
t,  to  enable  simulating  any  correlation  between  damage  and  its  detectability.  The  flight 
failure  time  and  abort  detection  time  distribution  function  are  the  only  means  used  in  the 
vehicle  mode  to  define  the  effects  of  damage.  Subsystem  failures  are  simulated  only  in  these 
terms. 
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Table  4 Might  Failure  Time  and  Abort  Detection 
Time  Probability  Functions. 


Failure 
time,  t 

Probability 
for  t 

Abort  dete 

ction  time,  r 

1 

3 

4 

Probability  for  r given  t 

1 

0.02 

1.0 

1.0 

1.0 

1.0 

0,04 

0.5 

1.0 

1.0 

1.0 

3 

0.05 

0.3 

O.b 

0.9 

1.0 

4 

1.00 

0.2 

0.4 

0.7 

1.0 

The  two  times  (t  and  r)  form  a doublet  that  is  used  as  the  state  variable  to  compare  to 
failure  and  abort  time  requirements.  I lie  MISDFM  simulation  allocates  a state  probability 
for  each  event  lor  all  combinations  of  t and  r.  The  state  probability  is  then  distributed  to 
modes  according  to  the  user's  mode  logic  which  defines  the  failure  and  abort  time 
requiivments  for  each  functional  mode.  A representative  example  of  mode  logic  is  given  in 
Table  5. 

Utilizing  a mission  with  four  mission  time  increments,  only  two  binary  bitts  are 
required  to  allow  description  of  the  four  time  requirements.  The  binary  numbers  and  their 
representative  constraints  are 

1 . 00  represents  no  constraint  on  the  system. 

2.  01  indicates  that  one  time  increment  must  remain  in  the  mission  before  system 
failure  or  the  need  to  abort  is  detected. 

3.  10  represents  two  time  increments  must  remain. 

4 I 1 is  three  remaining  time  increments. 

Therefore,  a 01  in  the  detection  time  column  (rl  of  Table  5 means  that  one  time  unit 
remains  before  the  need  to  abort  is  detected,  and  a 1 1 in  the  failure  time  column  (tl 
indicates  three  units  of  flight  time  remain  before  system  failure.  The  MISDFM  programs 
generate  all  combinations  of  the  state  vectors  and  compares  them  for  each  event  with  the 
mode  requirements  defined  by  the  user.  I ach  state  vector  that  meets  or  exceeds  the  mode 
requirement  becomes  associated  with  that  mode.  I ach  event  must  end  with  all  zeros  as  the 
requirement  li  e.,  the  most  degraded  mode)  so  the  probability  density  properties  for  each 
event  can  be  maintained.  The  probabilities  of  occurrence  for  each  state  are  accrued  to  each 
mode  and  event  where  they  satisfy  the  requirement.  Table  5 is  further  explained  as  follows 

Tvent  1.  There  are  no  failure  or  detection  time  requirements,  so  no  ones  appear  in 
the  time  requirements.  The  location  of  ones  in  this  array  is  analogous  to  the  location  of 
ones  m the  output  format  for  the  electronic  sample  case,  as  in  Table  I.  where 
Tquipment  Description”  is  replaced  by  “Time  Requirements”. 
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Event  2.  The  first  two  normal  modes  require  that  the  time  of  detection  exceeds 
the  previous  event  duration,  so  that  the  vehicle  will  have  arrived  at  the  current  event 
without  detecting  a need  to  abort.  All  time  requirements,  from  the  most  to  the  least 
stringent,  must  be  specified  for  each  mode  to  keep  out  probability  mass  belonging  to 
more-degraded  modes.  The  next  mode  is  also  a normal  mode,  since,  even  though  the 
damage  is  detected  within  the  previous  event  (detection  time  remaining  = zero),  the 
crew  determines  that  the  flight  time  remaining  equals  the  mission  time  remaining 
(three  units).  In  the  next  two  modes,  the  flight  time  is  less  than  the  mission  time,  so 
the  crew  elects  to  abort.  When  the  flight  time  remaining  after  the  event  is  reduced  to 
zero,  the  aircraft  has  landed  or  crashed. 

The  other  events  are  similar,  and  moving  toward  the  end  of  the  mission,  the  criterion 
for  continuing  to  fly  the  mission  when  damage  is  discovered,  is  relaxed.  When  the  aircraft,  in 
event  5,  has  no  better  abort  route  than  the  normal  mission  route,  the  abort  mode  is 
abolished. 


Table  5.  Representative  Example  of  Mode  Logic  for  Vehicle  Case. 


Event 

Event  description 

Mode 

Time  requirements 

Failure 

time, 

(binary) 

t 

Detection 

time, 

(binary) 

r 

1 

First  defensive  event 

Normal 

0 0 

0 0 

2 

Second  defensive  event 

Normal  A 

0 0 

1 0 

Normal  B 

0 0 

0 1 

Normal  C 

1 1 

0 0 

Abort  A 

1 0 

0 0 

Abort  B 

0 1 

0 0 

Down 

0 0 

0 0 

3 

Third  defensive  event 

Normal  A 

0 0 

1 0 

Normal  B 

0 0 

0 1 

Normal  C 

1 0 

0 0 

Abort  A 

0 1 

0 0 

Down 

0 0 

0 0 

5 

Fourth  defensive  event 

Normal  A 

0 0 

0 1 

Normal  B 

1 0 

0 0 

Normal  C 

0 1 

0 0 

Down 

0 0 

0 0 

6 

Landing  event 

Normal  A 

0 0 

0 1 

Normal  B 

1 0 

0 0 

Normal  C 

0 1 

0 0 

Down 

0 0 

0 0 

i 


j 
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Transition  Algorithm 

bach  element  of  the  transition  matrix  is,  as  in  electronic  mode,  a conditional 
probability  of  transition  from  a prior  state  (now  a doublet  tj,  rj)  to  a current  state  (tj.  r j). 
I he  cause  of  a transition  is  not  necessarily  subsystem  failure,  but  rather  progressive  loss  of 
flight  time  or  abort  detection  time  (relative  to  the  current  event  time). 

An  example  will  show  a transition  that  is  not  due  to  decrease  of  absolute  flight  time  or 
abort  time,  and  hence  is  not  due  to  any  threat  action.  Consider  the  following  transition: 


t r 

(binary)  (binary) 

System  state  i 0 1 0 1 

System  state  i 0 0 0 0 

Transition  probability  I’j/j  1 times  I 


The  transition  probability  is  unity  because  the  prior  remaining  flight  and  detection  time 
(due  to  all  prior  damage  exposures)  is  reflected  in  state  i,  and  must  be  decremented  by  one 
mission  time  unit  at  the  time  of  state  j.  There  are  no  other  possibilities. 

A more  complex  example  will  show  a transition  that  reflects  a decrease  in  absolute 
High!  and  detection  time: 


t r 


(binary) 

(binary) 

System  state  i 

1 0 

1 0 

System  state  j 

0 0 

0 0 

I ransition  probability  P,  j 

Pft  1 ) 

times 

Pdd.  I> 

where 

P,<l)  = P(t  = l)  = P(t-l  =0) 

Pd(l,  1 ) = P(  r = l/t=  1 ) = P(  r - 1 = 0/t  -1=0) 

t.  r are  the  times  of  (light  and  detection  associated  with  the  preceding  event  (and 
measured  from  that  time) 

Pj<  I ) is  the  probability  that  t = I 

P I ( 1 , 1)  is  the  conditional  probability  that  r = 1 given  that  t = 1 

The  transition  cannot  occur  as  a result  of  the  time  advance  alone,  since  it  would  result  in  a 
transition  to  the  state  (tj.  r j ) = ( 1 , I).  The  only  way  this  transition  can  take  place  is  that  the 
damage  exposure  just  concluded  resulted  in  a t = 1 and  a r = I.  so  that  t - l = r - l = 0 And 
the  probability  P(t  = I and  r = I ) = P(t  = I )l’(r  = 1/t  = I ) by  the  definition  of  conditional 
probability. 


I 
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The  whole  transition  algorithm  is  shown  in  Figure  17.  A state-pair  is  introduced  at  the 
top,  and  tested  to  reject  unwanted  states  and  save  computation  time.  Unwanted  states 
include:  ( 1 ) states  for  which  r > t (considered  invalid,  since  if  an  aircraft  is  forced  down,  the 
crew  is  sure  to  have  detected  the  need  by  that  time),  (2)  states  for  which  either  tj  or  Tj  <;0, 
in  which  case  the  aircraft  was  aborted  or  down  in  the  last  event  (there  is  no  need  to  carry 
them  forward),  and  (3)  states  for  which  the  current  (light  or  detection  time  has  increased 
(an  impossible  transition).  All  of  these  state  pairs  that  survive  are  subjected  to  one  of  four 
transition  probability  formulae  depending  on  the  magnitudes  of  the  decreases  in  t and  r in 
the  transition.  The  derivations  of  these  transition  probabilities  are  given  in  Table  6. 
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Table 
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List  of  Abbreviations  anil  Symbols 
(Mathematical  Model) 


Abbreviation 
or  symbol 

Equivalent  in 
simulation 
model 

— 

1 

Definition 

llnils 

CEP 

( EP(J) 

Circular  error  probable  (for  threat  weapon) 

dk 

PMMtM) 

Reliability  ofklh  subsystem  in  transition 

At 

DELTA! 

Elapsed  subss  stem  tune  since  last  event 

hours 

Ey  (N) 

ET 

Expected  number  ol  targets  killed  (N  assigned) 

l/.(k> 

I/ll  ) 

Vulnctability  index  associated  with  subsystem  k 

JCAP 

JCAP 

Number  of  mode  sequences 

K 

JCOUNT 

Number  of  states 

k 

M 

Subsystem  ordinal  number 

MLTM 

Mission  length  after  first  exposure 

Um'.r) 

CURVE(IJ) 

Designation  for  nuclear  damage  mechanism  type 

Depends 
on  m • 

MTBF 

MTBE(M) 

Mean  time  between  failure 

hours 

N 

NONE 

Number  of  events 

Pi 

PI(IJK) 

Stale  probability  for  prior  event 

pi.n 

PJ(I) 

Slate  probability  for  current  event 

Pj/i.n 

TRANS 

Transition  probability,  state  i to  state  j 

P)J 

PJ(JCOUNT) 

Probability  of  state  j (associated  with  mode 
sequence  J ) 

Pj.n 

PCAP(J) 

Probability  of  Jth  mode  sequence  in  nth  event 

Pk.J 

PK(J) 

Kill  probability  in  Jilt  delivery  mode  sequence 

PK.n 

SLIM 

Expected  kill  probability  (all  mode  sequences)  at 
event  n 

P|niss(k) 

PMISS(K.IJK) 

Probability  of  threat  miss  in  /.one  k (conven- 
tional) 

4k 

PCKILUM.L.K) 

Complement  of  (conventional  threats) 
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List  of  Abbreviations  and  Symbols  (contd) 


(Mathematical  Model) 


Abbreviation 
or  symbol 

Equivalent  in 
simulation 
model 

Definition 

Units 

Mk 

QPRM(M.IJK) 

Probability  of  threat  miss  in  kth  /.one  (nuclear) 

rk 

RM2 

Lethal  radius  for  kth  subsystem  (nuclear) 

feet 

rK 

R(K) 

Inner  radius  of  /.one  k (conventional) 

feet 

RM(k) 

RM(J) 

(Maximum)  lethal  radius  for  subsystem  k 
(nuclear  events) 

feet 

SPil 

FA(J) 

Expected  number  of  threat  weapons  arriving  in 
the  neighborhood  of  the  target  essentially  simul- 
taneously (nuclear  case) 

a 

SIGMA 

Standard  deviation  of  miss  distribution 

feet 

t 

NONE 

Tiroe-of-flight  due  to  a particular  threat 
exposure 

r 

NONE 

Abort-detection  time  due  to  a particular  threat 
exposure 

‘j 

LMAT(L) 

Current  time  of  flight  remaining 

rj 

LMAT(L) 

Current  abort-detection  time  remaining 

Vk 

PCSURV(M.L.K) 

kth  subsystem  survival  probability 

V*1:  famm* dot  IraWwc) 
n»  -J:  nmmnt  (mutremAq  ew>> 
m-ii  MM  <lW«q  in) 
m ”4:  MfiM  IcatorWtq  in) 
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SCHEMATIC  FLOW  DIAGRAMS 


This  section  describes  the  program  architecture,  logic,  and  arithmetic  sequence. 


OVERVIEW 

As  shown  in  Figure  18,  the  M1SDEM  program  lias  two  parts  called  Program  I and 
Program  2.  The  reason  tor  separating  the  programs  is  that  Program  1 is  needed  to  help  the 
user  set  up  Program  2.  Program  1 defines  all  the  possible  mode  sequences  that  can  be 
obtained  from  the  input  mode  sequence  logic.  (This  would  be  very  tedious  for  the  user  to 
do  manually.)  Program  l also  defines  a unique  mode  sequence  number  (tape  output). 
Program  1 processes  one  event  at  a time  until  event  data  are  exhausted.  The  user  reviews  the 
output  and  assigns  capabilities  (offensive  and  defensive)  to  each  mode  sequence  number  for 
every  event.  Program  2 first  computes  state  probabilities  reflecting  survivability  parameters 
from  the  last  event.  These  are  selectively  summed  (using  the  tape  input)  to  obtain  mode 
sequence  probabilities.  These  are  then  used  to  compute  expected  target  kill  probability  (all 
modes)  or  system  survivability  parameters  (depending  on  event  type)  based  on  the  input 
capabilities  for  the  current  event. 

Program  2 also  supplies  the  state  probability  distribution  at  each  event,  which  may  be 
used  to  initialize  the  state  probability  for  an  abort  flight  path  (in  a separate  application  of 
MISDEM). 

In  a typical  application  to  a mission.  Programs  I and  2 would  be  run  first  in  the 
electronics  mode,  and  second  in  the  vehicle  mode.  Probability  of  target  kill  from  the  first 
run  would  be  multiplied  by  arrival  probability  from  the  second  run  to  obtain  the  expected 
target  kill  probability. 


PROGRAM  ORGANIZATION 

Figures  19  and  20  show  the  principal  computation  and  logic  steps  in  the  program  in 
somewhat  greater  detail  than  Figure  IS,  and  list  all  input  and  output  parameters,  but  leave 
out  all  initialization  and  tape  control  steps  in  the  interest  of  simplicity  and  clarity. 


DETAILED  FLOW  DIAGRAMS 

Figures  21  through  49  am  detailed  schematic  diagrams  of  the  simulation  program. 
Figures  21  through  29  describe  Program  I,  and  Figures  30  through  49  describe  Program  2. 
Explanatory  text  is  provided  in  the  “Simulation  Model”  section  for  each  block  of  code 
shown  here. 
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[user] 


(ALL 

EVENTS) 


INPUTS:  EQUIPMENT  LIST  (EACH  MISSION) 

MODE  SEQUENCE  LOGIC  (EACH  EVENT) 

EVENT  DESCRIPTION  (EACH  EVENT) 

JJ 

^ PROGRAM  (select  NEXT  events ;l 

(determine  MODE  SEQUENCES  FOR  ALL  STATESH  ! 

I 1 ^ □ 

MODE  SEQUENCES  6 SUBSYSTEMS  USED  (ALL  EVENTS) 


PRINT 

OUTPUTS: 


) 


[user] 


1 


INPUTS: 

(EACH 

EVENT) 


CAPABILITIES  (WEAPON/TARGET  KILL  PROBABILITY 

(EACH  MODE  SEQUENCE ):]  THREAT  EFFECTIVENESS  PARAMETERS 


PROGRAM  I I 


SELECT  NEXT  EVEKtU- 


COMPUTE  STATE  PROBABILITIES  (EACH  EVENT) 


□ COMPUTE  MODE  SEQUENCE  PROBABILITIES  (EACH  EVENT) 


DETERMINE  SYSTEM  SURVIVABILITY  PARAMETERS 
(EACH  DEFENSIVE  EVENT) 


COMPUTE  OFFENSIVE  EFFECTIVENESS  PARAMETERS 
(EACH  OFFENSIVE  EVENT) 


pr1 


PRINT  OUTPUTS  (ALL  EVENTS): 

ELECTRONICS  CASE: 

MODE  SEQUENCE  PROBABILITIES 

TARGET  KILL  PROBABILITIES  (CURRENT  AND  CUMULATIVE) 
VEHICLE  CASE: 

AIRCRAFT  ARRIVAL,  ABORT,  AND  DOWN  PROBABILITIES 
CUMULATIVE  ARRIVAL  PROBABILITIES 


Figure  18.  MISDEM  Program  Simplified  Flow  Chart. 
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Figure  19.  Program  I Organization. 
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REFERENCE 

FIGURES 


DETERMINE  MODE  SEQUENCE 


ASSOCIATED  WITH  EVERY  SUBFUNCTION  AND  MODE  (MM)  IN  THE 
MOOE  SEQUENCE  LOGIC  IS  ANOTHER  SAM  WHICH  IS  TO  BE  TESTED 
NEXT  IF  THE  MISSION  DESCRIPTION  DOES  NOT  SATISFY  THE 
ENVIRONMENTAL  USE  RESTRICTIONS.  OTHERWISE.  THE  STATE 
IS  TESTED  AGAINST  THE  SUBSYSTEMS  REQUIREMENTS  OF  THE 
CURRENT  SAM.  IF  IT  PASSES  THIS  TEST.  THE  STATE  IS  NEXT  TESTED 
AGAINST  THE  HIGHEST  MODE  IN  THE  NEXT  SUBFUNCTION.  IF  IT 
FAILS  THIS  TEST.  THE  STATE  IS  NEXT  TESTED  AGAINST  THE  NEXT 
MODE  IN  THE  CURRENT  SUBFUNCTION.  THIS  PROCEDURE  IS 
CONTINUED  UNTIL  THE  MOOE  SEQUENCE  HAS  BEEN  COMPLETED. 


ASSIGN  A MODE  SEQUENCE  NUMBER 


AS  EACH  MODE  SEQUENCE  IS  GENERATED.  IT  IS  TESTED  FOR 
UNIQUENESS.  A NEW  (OR  OLD)  MODE  SEQUENCE  NUMBER  IS 
ASSIGNED  TO  EACH  STATE  NUMBER  AND  STORES  ON  TAPE  UNIT  4 
TO  BE  USED  IN  PROGRAM  TWO. 


IS  THIS  THE 
LAST  STATE? 


IS  THIS  THE 
LAST  EVENT? 


Figure  19.  Program  1 Organization  (contd.). 
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Figure  20.  Program  2 Organization. 
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Figure  20.  Program  2 Organization  (contd.). 
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Figure  20.  Program  2 Organization  (contd.). 
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COMPUTE  MODE  SEQUENCE  PROBASIUTICS 


ALL  THE  STATE  PROBABILITIES  ASSOCIATED  WITH  A MODE 
SEQUENCE  NUMBER  ARE  SUMMED.  THIS  IS  REPEATED  POR  ALL 
MODE  SEQUENCES.  EACH  OF  THESE  SUMS  REPRESENTS  THE 
PROBABILITY  OF  A MODE  SEQUENCE  IN  THE  CURRENT  EVENT. 


COMPUTE  MISSION  EFFECTIVENESS  (OFFENSIVE 
EVENTS  ONLY) 


COMPUTE  THE  EXPECTED  TARGET  KILL  PROBABILITY  BY  MULTI- 
PLYING THE  MODE  SEQUENCE  PROBABILITY  BY  ITS  ASSOCIATED 
(INPUT)  TARGET  KILL  PROBABILITY  (POR  THE  CURRENT  EVENT) 
AND  SUMMING  OVER  ALL  MODE  SEQUENCES.  THIS  STEP  IS 
EQUIVALENT  TO  EQUATION  2. 


DETERMINE  SYSTEM  SURVIVABILITY  PARAMETERS 
(DEFENSIVE  EVENTS  ONLY) 


SELECT  ONE  OF  THREE  KILL  ALGORITHMS  DEPENDING  ON  MISOEM 
MODE  AND  EVENT  TYPE 

ELECTRONICS  MODF:  QUICK  CONVENTIONAL  KILLS  OR  QUICK 
NUCLEAR  KILLS 

VEHICLE  MODE  SLOW  KILLS 


IS  THIS  THE 
LAST  EVENT? 


STOP 


Figure  20.  Program  2 Organization  (eonUl.). 
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I M 1 I *«i  I ! ! PROGRAM  vARIARll 


hi  .i(«  i Wu!  omoNs 


Rl  -l'  sr.ll"  !'!'•*  Hill  ION 


l 1 


WKI  H SVSUM  01  SCRIPT  ION 


MOO  MAN* 

inputs  OH 

It  HO 

siirs  i Kin 


ii  mis  suhsvstim  minim  numrtr  is 
CRI  All  R 1 MAN  IIIRKINI  COUNT,  RlPlACl 
W Mm  THAT  MVNBIR 


II  SURSYS T I M IS  I XCIUOI 0 I RON 
IONS  I 01  RAt  ION.  HAG  IT 


STORI  H>  SVSIIN 

OAMAuf  THRISHOLO  INFORMATION 
SUBSYSTlM  NAMl 


TOTAl  NOMRtR 


Of  SUBSYSTCMS < 22 


TOTAl  NKMRTR  Of  SOBSYSTINS  - 


Rl  AO  Nt  XI  CARO 


I NO  Of 
SYSTI M 
01  st  RIPT  ION 


WRITI  IRROR 
M|  ssAfi( 


FIQURE  22 

Figure  21.  MISDEM  Program  1 System  Description  Inputs. 
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Figure  22.  MISDEM  Program  1 - Event  Description  Inputs. 
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FIGURE  22 


FOR  A SU6FUNCTJ0W 


WAITE  EAAOA 
MESSAGE  


FIGURE  24 


Figure  23.  MISDEM  Program  1 - Read  and  Write 
Subfunctional  Flow. 
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StT  NUMBER  OF  ZEROS 


NUMBER  Of 

/EROS  < I T 


SPECIFY  IOC  AT  ION  OF 
LEFT  MOST  ZERO  AS  NO 


SET  ZERO  AT 
INDICATED  LOCATION 


NUMBER  OF 
ZEROS  <2? 


SPECIFY  LOCATION  Or 
ZERO  NO.  2 NEXT  TO 


SET  ZERO  AT 
INDICATED  LOCATION 


NUMBER  Or 
ZEROS  <8? 


SPECIFY  LOCATION  OF 
ZERO  NO.  8 NEXT  TO  , 


SET  ZERO  AT 
INOICATEO  LOCATION 


FI  CUBE  » 

Figure  24.  M1SDEM  Program  I - Initialize  for 
Mode  Sequence  Identification 
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are  mode  x 

MISSION 
REQUIREMENTS 
SAT  I SE I E D?  ^ 


^ ARE  \ 
r SUBSYSTEM 
REQUIREMENTS 
MET  BY  CURRENT 
\ STATE?  / 


IS  DESIGNATE  D 
SUBEUNCT I ON 
SAME  AS  CURRENT 
SUBFUNCTION? 


FIGURE  27 


Figure  26.  MISDEM  Program  1 
and  Subsystems  Used. 


- Define  Mode  Sequence 
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FIGURE  25 


FOR  A SUBFUNCTION 


FOR  A MODE 


STORE  DESIGNATED 
SUBFUNCTION  f.  MODE 


YES 


RECORD  COMFONENT  REQUIREMENTS 
IN  "SUBSYSTEMS  USED"  ARRAY 


RECORD  CURRENT  SUBFUNCTION  AND  MODE 
INDICES  IN  "MODE-SEQUENCE"  ARRAY 


STORE  DESIGNATED  SUBFUNCTION 
AND  MODE 


MORE 

SUBFUNCTIONS 


NO  MORE 
SUB FUNCTIONS 


MORE  MODES  FOR  THIS 
SUBFUNCTION 


NO  MORE  MODES 

FOR  THIS  SUBFUNCTION 


DEFINE  (NEXT)  MODE  SEQUENCE,  TENTATIVE  MODE  SEQ 
NO.,  AND  SUBSYSTEM  LIST  FOR  CURRENT  STATE 


IS  THIS  THE  X 
1ST  MODE  SEQUENCE? 


SELECT  TENTATIVE 
MODE  SEQUENCE  NO 


IS  CURRENT  FLOW 
PATH  IDENTICAL  TO 
EARLIER  MODE  SEQ? 


ASSIGN  CURRENT  STATE  NUMBER  TO 
THIS  MODE  SEQUENCE  NUMBER 


FIGURE  28 

Figure  27.  MISDFM  Program  I Assign  Mode  Sequence  Numbers. 
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CUAAINT  NUNBfA  01  l IKOS 


INC  Af  HINT 
LOCATION  Of 
ZlAO  NO  I 


/HA!>  zi  aoN 

NO.  I CONI 
AS  (AA  TO 
I NO  Of  UNi 
as  possitur 


/mas  a ao\ 

NO.  ; CONI 
AS  f AA  TO 
(NO  or  l INI 
As  AOSSlillT 


' ANT  NO A l > 
Zt  AOS  TO  A 007 


Figure  28.  M1SDEM  Program  1 - Reposition  Zeros  in 
State  Binary  Vector,  (KW). 
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FIGURE  28 


Figure  29.  MISDEM  Program  1 - Write  Mode  Sequence  and  Subsystems  Used. 
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FIGURE  31 


Figure  30.  MISDEM  Program  2 - System  Description  Inputs. 


1 


! I 


i 


! 


.1 


i 

Ji 

I 


. 


FIGURE  32 

Figure  31.  MISDEM  Program  2 - Event  Description  Inputs. 
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FIGURE  31 


FOR  EACH  DAMAGE  MECHANISM 


COMPARE  NEWLY  OBTAINED  LETHAL  RADIUS 
WITH  STORED  VALUE.  SELECT  MAXIMUM 
VALUE  AND  STORE 


MORE  SUBSYSTEMS  FOR 
THIS  DAMAGE  MECHANISM 


NO  MORE  SUBSYSTEMS  FOR 
THIS  DAMAGE  MECHANISM 


MORE  DAMAGE  MECHANISMS 


NO  MORE  DAMAGE 
MECHANISMS 


FIGURE  33 


Figure  32.  MISDEM  Program  2 - Compute  Lethal  Radius. 
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FIGURE  32 


FIGURE  34 


Figure  33.  M1SDEM  Program  2 - Assign  Vulnerability 
Index  Numbers. 
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FiaUMC  ft 


Figure  34.  MISDEM  Program  2 Compute  the  Reliability  of  Every 
Subsystem  and  Reshuffle  Reliability  Order. 
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INimuH  IN)  T I At  5t*ll  V|CI0«  . I 


Figure  35.  MISDFM  Program  2 Generate  a 
Prior  State  of  the  System. 


Figure  36.  MJSDEM  Program  2 - Test  for  Suppressed  Systems  in  Prior  State  and  Reorder  Subsystems. 
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FIGURE  36 


INITIALIZE  CURRENT  STATE  VECTOR 


NUMBER  OF 
ZEROS  <1? 


SET  ZERO  NO.  I AT 
INDICATED  LOCATION 


NUMBER  OF 
ZEROS  <2? 


SET  NUMBER  OF 
ZEROS  - 2 


SET  ZERO  NO.  2 AT 
INDICATED  LOCATION 


NUMBER  OF 
ZEROS  <8? 


SET  NUMBER  OF 
ZEROS  - 8 


SET  ZERO  NO.  8 AT 
INDICATED  LOCATION 


FIGURE  3« 


Figure  37.  MISDEM  Program  2 - Generate  Current 
State  of  the  System. 


IS  C*IT 


Figure  38.  MISDEM  Program  2 - Test  for  Suppressed  Subsystems  in 
Current  State  and  Reorder  Subsystems. 
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FIGURE  38 
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Figure  39.  MISDEM  Program  2 - Quick  Conventional  Damage. 
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. UP?  / 
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.SUBSYSTEM  DOWN? 


JTCG/AS-76-S-004 


REORDER  SUiSYSTEM  STATE 


INITIALIZE  Pj/i  • 1 
SET  FLAG  - I 


FOR  EACH  SUBSYSTEM 


•i/I  'j/1  * (,'V 

(SUBSYSTEM  KILLED; 

AFFIX  RELIABILITY  FACTOR) 


SUBSYSTEM  RESURRECTEO. 
IMPOSSIBLE  EVENT 
SET  Pj7|  - 0 


SET  FLAG  - I 
(SUBSYSTEM  SURVIVES) 


(SUBSYSTEM  STILL 
KILLED) 


YES  ADO  NUCLEAR  SUFFIX 


NUCLEAR  KILLS  STILL  POSSIBLE 
ADD  NUCLEAR  SUFFIX 
P - P 


(NUCLEAR  KILL 
NOT  POSSIBLE 


* REMAINING  SUBSYSTEMS  ARE 
AFFECTEO  ONLY  BY  RELIABILITY 


MORE  SUBSYSTEMS 


NO  MORE  SUBSYSTEMS 


FIGURE  41 


Figure  41.  MISDEM  Program  2 - Quick  Nuclear  Damage  and 
Reliability  and  Non-Defensive  Transitions. 
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FIGURE  41 


IS  NUMBER  OF 
ZEROS  . I OR  -8? 


CURRENT  NUMBER  Of  ZEROS 


INCREMENT 
LOCATION  OF 
ZERO  NO.  8 


INCREMENT 
LOCATION  OF 
ZERO  NO.  7 


/HAS  ZEROES 
NO.  7 CONE 
AS  FAR  TO 
END  OF  LINE 
AS  POSSIBLE ? 


ADO  ANOTHER 
ZERO 


' ANYMORE  ^ 
ZEROS  TO  ADO? 


FIGURE  43 


Figure  42.  MISDEM  Program  2 - Compute  State  Probabilities 
and  Reset  Current  State  Vector  Zeros. 
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FIGURE  42 


ABORl 

CAS£* 


CORAINT  NOMBtR  OF  ll  ROs 


INCRfMCNT 

location  OF 

JE  RO  NO  8 


/ *AS  ll  RO\ 
NO . / GONE 
AS  FAR  TO 
END  OF  UNt 
AS  POSSiPlf? 


A DP  ANOTMfR 


' AN>  MORI 
ZEROS  TO  App 


Figure  43.  MISDEM  Program  2 - Reset  Prior  State  Vector  Zeros, 
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FIGURE  43 


FIGURE  46 


Figure  44.  M1SDEM  Program  2 - Compute  and  Write 
Mode  Sequence  Probabilities. 
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FIGURE  44 


Figure  45.  M1SDEM  Program  2 - Compute  Offensive 
Events  Mission  Effectiveness. 


73 


READ  AND  WRITE  SUBSYSTEM  KILL  PROBABILITY 
FOR  ALL  SUBSYSTEMS  AND  All  OFFSET  TRAJECTORIES 


LAST  STATE 
NUMBER? 


SELECT  A STATE  NUNBl  R 


INITIALIZE  MISS  PROBABILITY  FOR  LARGEST  RAOIUS  TO  1.0 


CONVERT  CEP  TO  STANOARO  DEVIATION 


COMPUTE  TEST  VALUE  FOR  RADIUS  OF  ZONE  I DEN  I ?$' 


SELECT  A ZONE  RADIUS 


Figure  46.  MISDEM  Program  2 - Quick  Conventional  Kills. 
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FIGURE  45 


FIGURE  40 


FIGURE  47 


COMPUTE  ZONE  MISS  PROBABILITY  - FIRST  TERM 
MINUS  SECOND  TERM 


DECREMENT  MISS  PROBABILITY  FOR  LARGEST  RADIUS  BY 
CURRENT  MISS  PROBABILITY 


SET  SECOND  TERM  IN  MISS  PROBAR ' L I T Y EQUAL  TO  ZERO 
TO  AVOID  UNDERFLOW 


COMPUTE  SECOND  TERM  IN 
MISS  PRORABIl I T Y 


REINITIALIZE  FIRST  TERM  . SECOND  TERN 
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FIGURE  46 


AA 

Figure  47.  MISDEM  Program  2 - Slow  Kills. 
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SUBROUTINE  TWODNZ 


CALCULATE  SUBSCRIPT  OF  FIRST 

VALUE  IN  INDEPENDENT  VARIABLE  COLUMN 


/ S ^ 

* GIVEN  VALUE^ 

LESS  THAN 
LOWER  LIMITS 

VF  table? 


DEFINE  SUBSCRIPTS  FOR 
UPPER  AND  LOWER  VALUES 
IN  RESULT  INTERVAL 


SET  RESULTS 
EQUAL  TO 
LOWER  LIMIT 
OF  TABLE 


GIVEN  VALUE 
IS  LESS 
THAN  THIS 
NUMBER  IN 
TABLE 


GIVEN  X 
UPPER  AND  LOWER 
LIMITS  SEARCH 
FDR  INTERVAL 
IN  WHICH 
GIVEN  VALUE 
X LIES 


GIVEN  VALUE  IS 
EQUAL  TO  THIS 
NUMBER  IN  TABLE 


GIVEN  VALUE  LIES 
WITHIN  TABLE, 

CALCULATE  INTERPOLATION 
FACTOR 


GIVEN  VALUE 
IS  GREATER 
THAN  THIS 
NUMBER  IN 
TABLE 


GIVEN  VALUE  IS 

GREATER  THAN 

UPPER  LIMITS  OF  TABLE 


SET  RESULTS  EQUAL  TO 
UPPER  LIMIT 


INTERPOLATION  FACTOR  * 0 


DO  LINEAR  INTERPOLATION 
FOR  ALL  PARAMETERS 
REQUESTED 


Figure  49.  MISDEM  Program  2-Linear  Interpolation  Routine  (Subroutine  TWODNZ). 
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SIMULATION  SOURCE  DECK 


MISSION  DAMAGE  EFFECTIVENESS  MODEL  - PROGRAM  ONE 

MISDEM  - PGM  1 


DIMENSION  CNI3) tCNAMfl  8,23),DEF(3> , DNAME 18) , FNAME I 10,27,  10  It 

1 H3MTIA)  ,JES(  10,2  7,25  1 , JUI 23 . 25) ,KS 1 10 , 27 ) , KU < 23 ) ,KW( 23 ) , LEE!  10)  , 

2 LF ( 1 0) , LLF (27) , L L0( l 3 ,27 ) , LMAI 10 , 27 ,23 ) , LMAT ( 10, 27, 23 ) ,1  MATT  ( 23  I 

3 ,LQ< 10,2  7) ,MO( 10,271  , MI | 23 ) , MMQI 10,27 ) , MQ< 10, 27 ) ,OFF ( 3 1 , wEAPNt 7 ) 
INTEGER*2  IJNI  256) 

LOGICAL  MD, MOT 

EQUIVALENCE  (LMAI 1 ,1,1  I ,LMAT(1 ,1,1 ) ) 

DATA  BLK/AH  / , OFF  /A  HOFF  E , AHNS  I V , AHE  / .DEF/AHDEFE, AHNSIV, 

•AHE  / , OUT / 2HNO/ 

1 INITIALIZE 

LI  =27 
LCN8  *8 
LAO* 10 
LCEI *23 
L 72=  7 
I NCE 1*0 

00  20  1*1, LCEI 
MUI  1*1 

00  20  J*  1 ,LCN8 
CNAMEI J, I I =BLK 
20  CONTINUE 

RE  AD  (5  ,950)  NZT  ,MCR,MPR  , MAV 


2 READ  AND  WRITE  SYSTEM  DESCRIPTION 

WRITE! 6,9601 
WR ITEI 6* 970) 

DO  AO  1*1, LCEI 

RE AO  1 5, 980)  M,  (ONAME(J) , J* 1 , LCN8 ) , T 8FM, ( HDMT I M PRM > ,MPRM* 1 , A) ,ONOF F 
I F ( M .E  Q. 999  I GO  TO  60 

WR ITE ( 6, 990 ) M, I DNAME ( J I , J*1 , LC N8 ) , T 0FM, ( HDMT (MPRM  ) ,MPRM«  l, A) .ONOFF 
IFIM.GT.LCEI.OR.M.LT.l  > GO  TO  50 
IF(M.GT.INCEI) INCEI-M 
IF(ONOFF.EQ.OUT)MI ( I »*0 
DO  AO  J*  1 , LC N8 
CNAMEI J, I )*DNAMEI Jl 
AO  CONTINUE 

RE  AD ( 5 ,980)  M 
IFIM.EQ.999)  30  TO  60 
50  WR I TE ( 6, 1000) 

STOP 

60  CONTINUE 
REWIND  A 
70  CONTINUE 


1 
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READ  AND  WRITE  EVENT  DESCRIPTION 

RE AO (5,1010, END-8401  I E VENT , T2 , NC, (WEAPN(I I, 1-1,71, MOT 
WRITE  I 6, 1270) 

DO  71  N-1,3 
CN(Nt -OFF (Nl 

IF(NC.NE.O)  CNINI-DEF(N) 

71  CONTINUE 

WRITE(6, 10201  IE  VENT « ( CN( Nl ,N*l ,3  I , T2 , ( WEAPN( I ), 1*1,71,  MOT 
IHNC.NE.7)  WR  I TE  ( 6,1360) 

IF(NC.EQ.T)  WRI TE (6, 1 065  I 

READ  SUBFUNCTIONAL  FLOW  FOR  EVENT  AND  GENERATE  SUBSYSTEM 
REQUIREMENTS  VECTOR 

DO  250  1*1, L40 
DO  240  K -l » LI 
MD( I ,K)>. FALSE. 

RE AD (5*1070) LF( II ,LLF( K I ,LQ (I , K) ,L LQ( I , K ) , ( FNAME ( I , K, J I , J* 1 , 10  I , 

* MQ(I , K) ,MMQ( I , K ) ,MD( I,K), (LMAT ( I,K,JI ,J-l,  INCE I ) 

L-0 

DO  90  J* 1 , I NCE I 

IF(LMAT{ I ,K, J) .EQ.OIGO  TO  100 
0 L-L*  1 

IF(LMAT( I ,K, J) .EQ.L1GD  TO  90 
LHATTIU  -0 
GO  TO  80 
0 LMATTCLl-l 
GO  TO  140 
00  J-L*l 

IF ( J.GT.INCEI1G0  TO  140 
DO  110  L*J,INCEI 
10  LMATTILI-0 
40  00  150  L-l.INCEI 

50  LMA( I , K, L)  -LMATTILI 

WRITE  SUBFUNCTIONAL  FLOW 

IF(NC.NE.T) 

1WRITE(6,1080)LF( 1 1 , LLF I K) , L0( I ,K),LLQ( I,K), ( FNAME ( I ,K , J ) , J- 1 , 1 0)  , 

* MQ ( I ,K) ,MMQ( I , K ) ,M0( I * K ) , (LMAT ( I , K , J ) , J* 1 , INCE I) 

IF (NC.EQ.7)  WR I TE ( 6, 1085)  LF ( I ) , LLF (Kl , LQ( I , K ) , U 0( I , K ) , 

* I FNAME (I, K, J I, J- 1,10) , (LMATI I,K, J),J-1,INCEI ) 

IFILLF(K)  .EQ.99)  GO  TO  230 

IF (LF ( I ) .EQ.999)  GO  TO  260 

GO  TO  240 
30  LEEdl-K-1 
GO  TO  250 
40  CONTINUE 
50  CONTINUE 
I-L40M 
REA0(5,1070)M 
IF(M. EQ.999)  GO  TO  260 
WRITE! 6, 10901 
STOP 
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C 

C6  INITIALIZE  FOR  MODE  SEQUENCE  IDENTIFICATION 

C 

260  LE-I-l 
LQSVL-0 
LLQS-0 
NZ*0 

DO  2 90  J*1 * LAO 
DO  290  K*l,Ll 
KSI J.KI-0 
DO  290  I *1 *L72 
JESt  J.K, I 1*0 
290  CONTINUE 

DO  300  1 *1* LCE I 
KWII  1*0 
300  CONTINUE 

DO  310  JCOUNT *1*256 
IJNI  JCOUNT )*0 
310  CONTINUE 
JCOUNT  *0 
JC  AP  =0 
C 

C 7 GENERATE  FINAL  STATE  OF  THE  SYSTEM  - BINARY  VECTOR  KW 

C 

330  CONTINUE 

DO  3*0  1*1. I NC El 
KWI I 1-1 
3*0  CONTINUE 

IFINZ.LT. 1IG0  TO  *00 
LLl-l 

350  KW(LL1»*0 

IFINZ.LT. 2IG0  TO  *00 
LL2-LL 1+ 1 
360  KMILL2 1*0 

IFINZ.LT. 3IG0  TO  *00 
LL3-LL2*! 

370  KMILL3I-0 

IFINZ.LT. *)G0  TO  *00 
LL**LL3*1 
380  KWILLA  )*0 

IFINZ.LT. 5IG0  TO  *00 
LL5-LLAM 
390  KMILL5I-0 

IFINZ.LT. 61  GO  TO  *00 
LL6*LL5*1 

391  KMILL6I-0 
IFINZ.LT.?)  GO  TO  *00 

LL  7*  LL6* 1 

392  KWILL7I-0 
IFINZ.LT. 8)  GO  TO  *00 
LL  8*LL  ?♦ 1 

393  KWILL8I-0 
*00  CONTINUE 
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8 TEST  TOR  SUPPRESSED  SUBSYSTEMS 

DO  410  MM-l « I NCE ( 

IFCMIIMMI.NE.OI GO  TO  410 
IFIKWIMM1.EQ.11G0  TO  720 
410  CONTINUE 

IF INC.EQ.7JG0  TO  405 
IFIMCR.EQ.0IG0  TO  402 
IFIKWI MCR I.EQ.OIGO  TO  720 
402  IFIMPR.6Q.OIGO  TO  404 

IFIKWIMPRI.EQ.OIGO  TO  720 

404  IFIMAV.EQ.OIGO  TO  405 
IF|KW(MAV).EQ.3)G0  TO  720 

405  CONTINUE 
JCOUNT  * JCOUNT  »1 

9 OEFINE  MODE  SEQUENCE  AND  SUBSYSTEMS  USED 

DO  530  l-l.LCEI 
53C  KUIII-0 
L TEST*  1 

DO  640  L * l,LE 
LEQ-LEEILI 
DO  540  K-  l.LEQ 
540  KSI L # K ) *0 

I F ( L .EQ.LQSVL.OR.LTEST.EQ. II  GO  TO  550 
GO  TO  640 

550  00  630  K = 1 . LEQ 

IFIK.GE.LLQS. OR. LTEST. EQ.l ) GO  TO  560 
GO  TO  630 

560  IFI.NOT.MOT. AND.MOIL.Kl I GO  TO  580 
GO  TO  590 
500  LQSVL*MQ< L,KI 

llqs-mmqu.ki 

LTEST-0 

IFILQSVL.EQ.LI  GO  TO  630 
GO  TO  640 
59D  LTEST*1 

00  600  M*1»INCEI 

IF(KW< M| .LT.LMA<L,K,MI  I GO  TO  620 
630  CONTINUE 

00  610  M*1»INCEI 

, IF ( KW( Ml .GT.LNAILtKfM)  . OR. KM (Ml . EQ.O I GO  TO  610 
IFIKMIMI.EQ.LMAIL.K.MI  .AND.KMIMI.EQ.il  KU(M|*1 
610  CONTINUE 
KSIL.KI-1 
LQSVL»LQ(L,KI 
LLQS-LLQIUKI 
LTEST-0 

620  IFILQSVL.GT.LIGO  TO  64D 
630  CONTINUE 
640  CONTINUE 
J-JCAP*! 
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C 

CIO  ASSIGN  MODE  SEQUENCE  NUMBERS 

C 

650  J=J-l 

IF(J.NE.O)  GO  TO  680 

JCAP*JCAP*1 

J^JCAP 

DO  660  L * I > IE 
LEQ=LEE( L ) 

DO  660  K* 1 * LEQ 
JES(  L,K,  J)  *KS(  L , K ) 

660  CONTINUE 

DO  670  I U = l • I NCE  I 
JU( IU,J)=KU( IU) 

670  CONTINUE 
GO  TO  700 

680  DO  690  L * I,LE 
LEQ=LEE(L) 

DO  690  K*l,  LEQ 

IE( JESIL.K, J).NE.KS(L,K)  I GO  TO  650 

690  CONTINUE 
700  CONTINUE 

I JN( JCOUNT ) = J 

C 

Cil  REPOSITION  Z EROS  IN  STATE  BINARY  VECTOR  (KWI 

C 

720  Ic ( NZ . LT. I. OR. NZ.  GT .8) GO  TO  790 

GO  TO  (770,760,750,743,730,728,726,724),  NZ 
724  KW(LL8>=1 
LL8=LL8+L 

IF((LL8-B).LE.  (INCEl-NZH  GO  TO  393 
726  KW(LL7)*l 
LL  7=LL  7*  l 

IF(( LL7-7I.LE. ( INCEI-NZ) ) GO  TO  392 
728  KW( L L6 ) * 1 
LL6=LL6*l 

IF<<  U6-6I.LE.  (INCEI-NZ))  GO  TO  391 
730  KW( LL  5 ) * 1 
LL  5*LL5* 1 

IF((LL5-5).LE.(INCEI-NZ))GO  TO  390 
740  KW(LL4)*1 

LL4*LL4» 1 

IF((  LL4-4KLE.  (INCEI-NZ)  )G0  TO  380 
750  KW  ( L L3 ) * I 
LL  3*LL3* 1 

IF(( LL 3-3). LE. (INCEI-NZ)) GO  TO  370 
760  KW(LL2 1*1 
LL2*LL2*-1 

IF(( LL2-2).LE. (INCEI-NZ) )G0  TO  360 
77C  KM( LL I ) = l 
LL1*LL1«-1 

lF(CLLl-I).LE.(INCEI-NZ))GO  TO  350 
790  CONTINUE 
NZ*NZ*l 

IF(NZ.LE.NZT)30  TO  330 
NZ  *0 
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C 

C 12  WRITE  MOOE  SEQUENCE  AND  SUBSYSTEMS  USEO 

C 

DO  820  J=1,JCAP 

WRITE16, 1170) 

WRITEI6,  11*0)  J 
DO  800  L * l *LE 
LEQ=LEEIL) 

00  800  K*1,LEQ 

IF! JESIL.K, J).EQ.O>  S3  TO  800 
WRITEI6, 1150)  IFNAMEIL*K,I  » • I *1 , 10> 

800  CONTINUE 

IFINC.NE.7)  WRITE(6,1180) 

IF(NC.EQ.7)  WRITEI6,!! 85) 

DO  810  1=1 » I NCEI 

IF( JU(M, J).NE.  1)  GO  T]  810 

WRITE  (6.  1190)  ICNAMEU.M),  K=lf  ICN8) 

810  CONTINUE 
820  CONTINUE 

13  WRITE  MODE  SEQUENCE  ON  TAPE  FOR  USE  IN  PROGRAM  2 

WRITE  I 4# 1280)  II JN( IJK) , l JK= 1 , JCOUNT ) 

GO  TO  70 
8*0  END  FILE  * 

REWIND  * 

STOP 

I*  DEFINE  FORMATS 

950  FORMAT (1013) 

960  FORMAT  ( 26H1  S YSTE  M CONFIGURATION) 

9 70  FORMAT  1 1H0, 8X, 9HEQUIPME NT, 25X,  *HMTBF,7X,7HTH0MI 1 ),6X, 7HTHDMI2), 

♦ 7X»7HTHDM(3) »8X, 7HTHQW I * I/56X • 3HG  0,1 IX, IHN, 12X, 1HB, 15X, 1HT//IX) 

9 80  F0RMATU3,IX,8A*,F10.2/*E12.5,T65,  A2) 

990  FORMAT ( 1H  , IX , 1 3 , 2 X , 8A * , FI 0. 2 . 2X ,* I E12. 5 ,2 X ) , 3X, A2 ) 

1000  FORMAT  ( 50H0  TOO  MANY  INPUTS  OR  M IS  OUTSIDE  ALLOWABLE  RANGE) 

1010  FORMAT (I  3, IX, Ft. 2, IX, I 3, IX, 7 A*, 2X, Li) 

1020  F0RMATI1HU,5X,9HEVENT  NO., 12, 1X,3HIS  , 3AA/ 5X, 15HE VENT  OCCURRED  , 
*F6.2, 

*20H  HOURS  AFTER  TAKEOF F/5X , 21HEVENT  DESCRIPTION  IS  ,7AA/5X,L1) 

1060  FORMAT  ( IHO , 7X , 1 6HSUBF  UNCT I QN7  MODE , 32X, 2 1HEQUI PMENT  DESCR I PT ION, 1 2 

* X, 1 8HM I SSI  ON  OESCRIPTOR/1X) 

1065  F0RMATI1H0,  7X , 16HSUBFUNCT I ON/ MODE, 32X , 35HFL IGHT  AND  DETECTION  TIM 

*E  REMAINING) 

1070  FORMAT I*I3,1X,10A*,1X,2I3,1X«L1/23I3) 

1080  FORMATIlX,2I3,2X,2I3,2X,l0A*,T86,2I3,12X,Ll,T59,23U) 

1085  FORM ATIIX,2I3,2X, 213, 2X,10A*,  T59.23U) 

l C90  FORMAT  I29H0  TOO  MANY  INPUTS  FOR  F ARRAY) 

11*0  F0RMATI1H0, 11X , I 3) 

1150  FORMATIIH  , *5X , 10A* I 

1170  FORMATIIH  ,5X,16HMOOE  SEQUENCE  NO* 16X, 18HSUBFUNCT IONAL  FLOW) 

1180  FORMAT  1 1H0.58X ,15H SUBS YSTE MS  USED/1X) 

1185  FORMAT I1H0,58X ,31H ACTUAL  AND  APPARENT  FLIGHT  TIME) 

1190  FORMATIIH  ,60X,8A*) 

1270  FORM  AT (2  3H1  EVENT  DESCRIPTION) 

1280  FORMATI25I3) 

END 
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C «M*M**M»Mt*M*MtM**M»M*n*»*M**M*MMM***M**MM**MMMM 

c 

C MISSION  DAMAGE  EFFECTIVENESS  MODEL  - PROGRAM  TWO 

C MISDEM  - PGM  2 

C 

C ************************** *****************«************************** 


DIMENSION  CEP  I 50) ,CN(3» .CURVE! A, 30  I , DEF ( 3) , ONAMEI 8 ) , FAI 50 ) , 

1 FTDNN(16fI6) *FTFNNI  16 ) , GRAPH C 30 ) ,H0MT (A ) , ISUBI A ) , IT  I 231 , ITT  I 23  I 

2 ,IV<23) ,12123) ,JZI23) ,KW(23) ,KWW(23) ,MI( 23)  , MT BFI 23) .OFF  I 3) , 

3 PC API  50)  .PC KILL! 23,  10 ,27 > .PCSUR VI 23, 10. 27 ) , P 1 1 256 ) ,P J< 2 56 ) , 

A PKI50) ,PM(23) ,PMM(2  3),PMISS(27,256),PTDN(16, 16  I , PTFN 1 16  ) , 

5 00 ( 256), QPR MI  23, 256) ,RI 10) , R A 123 ) , RESUL 1 2 ) , RMI 2A ) , THDMI 23, A), 

6 TIME  FI  23) ,TI MENI 23) » WEAPN ( 7 ) 

INTEGER*2  IJNI256) 

LOGICAL  MOT 

OATA  OFF /AHOFF E , AHNSI V , AHE  / , OEF /AHOEFE , AHNS I V, AHE  /.0UT/2MN0/ 

C 

Cl  INITIALIZE 

C 

LCN8-8 
LCEI*23 
L 72* AO 
INCE  1*0 

00  30  I*1,LCEI 
MTBFII 1*0.0 
PM( I »*0.0 
Mill  ) * 1 

ivm»i 

IZI I )*I 

00  20  MPRM* 1 ,A 
THOMI I ,MPRM) *0 
20  CONTINUE 
30  CONTINUE 

RE  AO  (5,1110)  NZT , MCR, MPR , MAV , MLTH.N ABORT 

C 

C2  READ  ANO  WRITE  SYSTEM  DESCRIPTION 

C 

WRITE  (6,1120) 

WRITE  (6,1130) 

DO  50  I-l.LCEI 

RE  AO  I 5, 11A0)  M, ( ONAMEI J), J-l , LCN8 ) ,TBFM, (HDMT (MPRM), MPRM* 1,A),TMN 
1 , TMF ,ONOFF 

IF  I M. EQ. 999)  GO  TO  70 

WRITE  (6,1150)  M, I ONANE I J ) » J*1 » LCN8 I ,T  BFM, (HOMT (MPRM ) »MPRM* 1 , A ) 

IF  (M.GT.LCEI.OR.M.LT. 1)  GO  TO  60 
IF  (M.GT.INCEI)  INCE I* I 
IF  (ONOFF.EQ.OUT)  MI  1 1 ) -0 
00  AO  MPRM* 1 ,A 
THOMI I ,MPRM) *HDMT( MPRM I 
AO  CONTINUE 

MTBFI I )*TBFM 
TIMENl  D-TMN 
TIMEFI I ) *TMF 
50  CONTINUE 

READ  (5,1 1A0)  M 
IF  (M.EQ.999)  GO  TO  70 
60  WRITE  (6,1160) 

STOP 

70  CONTINUE 
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3 INITIALIZE 

ET-0.0 
' 1 1 -0 
T1*0.0 
REMIND  A 
REWIND  3 

4 INITIALIZE 

80  CONTINUE 

00  85  HM'INCEI 
RN(MI*0.0 
85  CONTINUE 

5 READ  ANO  WRITE  EVENT  DESCRIPTION 

RE AO  (5, l 1 70 » END* 1 100)  I EVENT , T2 ,NC , I WEAPNI I > . 1*1 , 7» , MOT 
WRITE  (6,12301 
00  90  N*1 ,3 
CN( N )*0FF IN) 

IF  (NC.NE.O)  CN(N)*OEF(N) 

90  CONTINUE 

WRITE  (6,11801  IE VENT, (CN(N),N*1,3),T2« (WEAPNI I),  1*1, 7),  MOT 
00  100  L*l, I NCEI 
IZ(L)*L 
100  CONTINUE 

IF  (NC.EQ.O.OR.NC.EQ.S  .OR.  NC.EQ.7)  GO  TO  200 

6 COMPUTE  LETHAL  RA01 l ANO  ASSIGN  INDEX 

00  120  I * 1 , NC 

REAO  (5,1260)  NPOINT 

N2-2*NPOINT 

READ  15,1220)  I SUB  1 1 ) , I CURVE! I , J ) , J-l ,N2 ) 

IS-ISUB(I) 

DO  120  N*i, I NCEI 
00  110  IJ-1.N2 
GRAPH! IJ)*CURVE(  I, I J) 

110  CONTINUE 

CALL  TWODNZ  (THOM! M, IS ) , 1, GRAPH, NPOINT, 2.RESUL ) 
RA(MI»RESUL(2I 
RM(M)*ANAXl(RA(NI ,RM(M) I 
120  CONTINUE 

INE-INCEU1 

RMIINEI-0.0 

00  160  LOO*1 , I NCE I 

tNIT«LOO*l 

00  130  J*INIT»INE 

IF  I RM(J).LE.RNILOO))  GO  TO  130 

STORE-RMI J) 

RM( JI*RN( LOOI 
RM (LOO) « STORE 
l TORE* l Z ( LOO) 

IZ!LOO)-IZ( J) 

IZ ( J ) * I TORE 
13C  CONTINUE 
160  CONTINUE 

WRITE  (6,1250)  (RM( M) , M*1 , I NE ) 
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C 

C?  COMPUTE  SUBSYSTEM  RELIABILITY  IN  TRANSITION 

C 

200  CONTINUE 

DO  2 30  M=l, INC  El 
PM|MI=1.0 

IE  (MMMI.EQ.OI  GO  TO  230 
IF  ( T2.LE.TIMENIMI I GO  TO  230 
IE  I T2.GE.TIMEFIMI » GO  TO  210 
DELTAT=T2-T1 

IE  I Tl.LE.TIMENIMI  1 0ELTAT*T2-T IMENIMI 
GO  TO  220 
210  DELTAT=0.0 

IF  I TIMEFIMI.GT.Tlt  DE LTAT=T I MEE IM » -T l 
220  PM(MI=1.0-0ELTAT/MTBF( M I *0E LTAT**2/ (2.0*MTBF(MI**2 I 
230  CONTINUE 
T1  = T2 

DO  240  L = 1 » INCEI 
JV  = IVILI 
PMMI  L ) =PMI  JV) 

240  CONTINUE 
C 

C8  INITIALIZE  STATE  AND  STATE  PROBABILITIES 

C 

NZ  *0 
NZ  1 = 0 

DO  250  I *1 f LCE I 
KM( I 1=0 
250  CONTINUE 

DO  260  J = 1 « L72 
PC API J I =0  >0 
260  CONTINUE 
JCAP=0 

IF  I IEVENT.NE.il  GO  TO  290 
ICOUNT  = 1 

DO  210  M=l, INCEI 

IT(MI=1 

ITT! Ml =1 

IE  IMIIMI.EO.il  GO  TO  270 
IT(M)=0 
ITT! Ml =0 
2 70  CONTINUE 
C 

C9  ABORT  FLOW  CONTROL 

C 

IE ( N ABORT .NE .0  I GO  TO  275 
GO  TO  278 
C 

CIO  INITIALIZE  ABORT  STATES 

C 

275  lF  (I  EVENT. NE.l)  GO  TO  230 

REA0I3IJEVENT, JCOUNT.I PJILI ,L=l, JCOUNTI 
I F I JE VENT .NE .N ABORT  I GO  TO  275 
DO  2 77  I * 1 * JCOUNT 
Pill »=PJII» 

277  PJII 1=0.0 
GO  TO  305 
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11  INITIALIZE  NORMAL  STATES 

278  DO  280  L *1 1 2 56 
PI  I L ) - 1*0 
PJIL >*0.0 
280  CONTINUE 
GO  TO  310 

290  00  300  l =1 « JCOUNT 
pin  ) =p j 1 1 ) 
pj< i ) *0.0 
300  CONTINUE 
305  CONTINUE 
ICOUNT  *0 

Cl2  GENERATE  A PRIOR  STATE  OF  THE  SYSTEM  (ITI 

C 

310  I FI NABORT . EQ.O )G0  TO  315 
GO  TO  320 
315  CCNTINUE 

IFIIEVENT.EQ.IIGO  TO  660 
320  CONTINUE 

00  330  I*1,INCEI 
IT<I)*1 
330  CONTINUE 

IF  (NZ.LT.1I  GO  TO  613 
1 1 1- 1 

360  IT(tIl>=0 

IF  (NZ.LT.2I  GO  TO  613 
112=111*1 
370  ITI! 12  ) = 0 

IF  (NZ.LT.3)  GO  TO  610 
113=112*1 
380  I Tl 1 1 3 1=0 

IF  (NZ.LT.6I  GO  TO  613 
1 1 6=  1 1 3*  l 
390  l T| I 16  » =0 

IF  INZ.LT.5)  GO  TO  610 
115=116*1 
600  I TC I 15  1*0 

IFINZ.LT. 6IG0  TO  610 
116=115*1 

605  I T( 1 16  1*0 
IFINZ.LT. 7)G0  TO  610 
1 1 7*  1 1 6*  1 

606  ITIII71*0 
IFINZ.LT. 8IG0  TO  610 
118*117*1 

607  1 T C 1 1 S I *0 
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C 

Cl3 

C 

410 


420 


422 

424 

425 

C 

C 14 

c 


430 

440 

C 

C 15 
C 

450 


460 


• 490 


500 


510 


520 


530 


535 


536 


537 

540 


TEST  FOR  DEIETEO  OR  CRITICAL  SUBSYSTEMS 
CONTINUE 

00  420  MM*1 • I NCE I 
IF  (MI (MM1.EQ. II  GO  TO  420 
IF  ( ITIMMI.EO. 11  GO  TO  840 
CONTINUE 

IF(NC.E0.7I  GO  TO  425 
IF (NCR .EO.OI GO  TO  422 
IF( I T( MCR1.EQ.0)  GO  TO  840 
IF ( M PR .E  Q.O I GO  TO  424 
IF(ITIMPR).EQ.O)  GO  TO  840 
IF ( MAV  .EQ.O IGO  TO  425 
IFK  TIMAVl.EQ.OI  GO  TO  840 
1C  OUNT  x I C OUNT  «•  1 

RESHUFFLE  SUBSYSTEM  ORDER 

00  430  L 3 1 » I NC  E I 
JV*I VILI 
ITT(L) *IT(JV> 

CONTINUE 
JC  OUNT  =0 

GENERATE  A CURRENT  STATE  OF  THE  SYSTEM  (KW) 


CONTINUE 
DO  460  1*1, INCEI 
KW( 11=1 
CONTINUE 

IF  ( N2I.LT. 11  GO  TO  540 
LL  1 * 1 
KWILL1 1*0 

IF  (NZI.LT. 21  GO  TO  540 
LL2  = LLU1 
Kw(LL2 1*0 

IF  ( NZI.LT.3)  GO  TO  540 

LL3=LL2*l 

KW(LL3 1*0 

IF  (NZI.LT.4I  GO  TO  540 

LL4*LL3H 

KW (L  L4 1 = 0 

IF  (NZt.LT. 5)  GO  TO  540 

LL5*LL4*1 

KW(LL51*0 

IFINZI.LT.6IG0  TO  540 

LL6=LL5*1 

KW ( L L6 1*0 

IF(NZI .LT.71G0  TO  540 

LL7*LL6«-1 

KWILL7 1*0 

IFINZI  .LT.81G0  TO  540 

LL8*LL7M 

KW ( L L 8 1 *0 

CONTINUE 


555 

C 

C17 

c 


560 

C 

C18 

C 
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TEST  FOR  SUPPRESSED  SUBSYSTEMS 

00  550  MM*1 , INCE I 
IF  (MI(MM).EQ.ll  GO  TO  550 
IF  ( KWIMMt.EQ. 1)  GO  TO  T10 
CONTINUE 

IFINC.EQ.7IG0  TO  555 
IFIMCR.EQ.OIGO  TO  552 
IFIKW(HCR>.EQ.O)GO  TO  710 
IFIMPP.EQ.OIGO  TO  556 
IFIKWIMPRI.EQ. 0) GO  TO  710 
IFIMAV.EQ.OIGO  TO  555 
IF(KU(MAVI.EQ.O)GO  TO  710 
JCOUNT  aJCOUNT* 1 

RESHUFFLE  SUBSYSTEM  ORDER 

DO  560  L*1,INCEI 
JV-IVIU 
KMW(U*KM(  JVI 
CONTINUE 

COMPUTE  THE  SYSTEM  STATE  TRANSITION  PROBABILITY 
IF  I II.NE.6)  GO  TO  1500 

QUICK  CONVENTIONAL  THREAT  DAMAGE  AND  RELIABILITY 

SUMTRK-0.0 

DO  620  K*1  * KMA X 

SUMTRA'0.0 

00  610  L*1,LMAX 

TRANS-1.0 

DO  600  M*l, INCEI 

PC  SUR VC  M, L , K)  * 1.0  - PCKI LL(M,L,K) 

IF  ( ITT(M)-KMWIM) I 570,580,590 

TRANS*0.0 

GO  TO  620 

IF  I ITT! MI.EQ.O.OI  GO  TO  600 
TRANS*PCSURVIM,L,M*PMMIM)*TRANS 
GO  TO '600 

TRANS=( l.O-PCSURVI M,L,K)*PMM(MI l*TRANS 
CONTINUE 

SUMTRA*  TR ANS+SUMTR  A 
CONTINUE 

SUMTRKM  SUMTRA/LMAX)*PMISSIK,I COUNT  I ♦SUMTRK 
CONTINUE 
TRANS*  SUMTRK 
GO  TO  700 

IFtII.NE.7l  GO  TO  630 
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C 

C20  SLOW  THREAT  DAMAGE 

C 

TRANS  - 1.0 
NTFN«0 
I TFN*0 

1 1 NC  F l *1  NCE  1 12 
00  1510  L-1.1INCEI 

IFIKMMlLI.EQ.il  NT  FN*  2*  Mil  NCEI/21 -LI *NTFN 
1510  IF  I 1TT1L l.EQ.l I I TFN-2** 111 NCE I /2»-L IMTFN 
NTDN-0 
I TDN*0 

JINCEi  -1  I NC  E I ♦ 1 

00  1520  L-J1NCEI • 1 NCE I 

IF (KMM1L l.EQ.l INTON-2* M INCEI-L I *NTON 
l ‘20  IF (I TT(L l.EQ.l I ITDN»2*M INCE l-L I ♦! TON 
NT  TF  -NTFNM 

1 T TF  ■ I TFNA 1 
N T TO  *NTON*  l 
I TTO»I TDNM 

IF( ITTD.GT. 1TTF.0R.NTTD.GT.NTTFIG0  TO  1605 
IF| I TTF.LE. l.OR. I TTO.LE.ll  GO  TO  1605 
IFINTTF.GT.1ITTF-1I.OR.NTTO.GT.IITTO-1HGO  TO  1605 
GO  TO  1610 
16C5  TRAN  SI -0 

GO  TO  1720 
1610  TRANS1 *0 

IF (NTTF-I I TTF- II  1 16  70.  1620 • 1620 
162  0 IF  I NTTD-I I TTO- 11  11630,1625,1625 

1625  DO  1627  l»NTTF,MLTH 
SUMPTD-0 

DO  1626  J-NTTD.I 

1626  SUMP  TO *P TON  I I , J I ♦SUMPT  0 

162  7 TRANSl-PTFNIl) *SUMPTO* TRANS  1 
GO  TO  1720 

163C  00  1631  !*NTTF,MLTH 
1631  TRAN  SI 3PTFN ( I I APTDNI I , NTTOI +TR ANSI 
GO  TO  1720 

16  7C  IF1NTTD-1 ITTD-ll  1 1680,  1675, 1675 

1675  DO  1676  J»NTTD ,NTTF 

1676  TRAN  SI *PTDN(NTTF ,J|*TRANSl 
TRANSl-PTFNINT TF I* TRAN  SI 
GO  TO  1720 

1 68C  TRAN  SI »PTFNI NT  TF I *PTDN I NTTF , NTTOI 
1720  TR ANSaTRANSl *TRANS 
GO  TO  700 
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21  QUICK  NUCLEAR  DAMAGE  AND  RELIABILITY 

I ALSO  NON-DEFENSIVE  TRANSITIONS! 

630  TRANS-1.0 
NF*1 

IFII I.EQ.OI  GO  TO  6*0 
TRAN S»QO 1 1 COUNT  I 
NF-0 

6*0  CONTINUE 

DO  690  M-l.INCEl 
IF  I ITTIMI-KMWtM)  I 6S3.660.6T0 
6 SO  TR AN S*0. 0 
GO  TO  TOO 

660  IF  I ITTIMI.EQ.O)  GO  TO  680 
NF«1 

TRANS*TRANS*PMM|M| 

GO  TO  690 

6 TC  TRANS»TRANS*<  1 .O-PMMIM)) 

68C  IF  INF.EQ.Ol  TRANS  »TRA  N$*QPRM I M, I COUNT ) 

69C  CONTINUE 

22  COMPUTE  THE  STATE  PROBABILITIES 

TOC  PJt  JCOUNTI- TRANS  *P1 1 ICOUNTMPJt  JCOUNTI 

23  RESET  CURRENT  STATF  I KM  I ZERO  LOCATIONS 

T10  IF  INZI.LT. l.OR.NZl.GT.8)  GO  TO  820 

GO  TO  (800*780.760<7*0.720.T15«71*«713I.NZI 
113  KMILLBI-l 
LL8-LL8»1 

IFII LIB-81. LE. I INCEI-NZIIIGO  TO  S3  T 
11*  KWILL7I-1 
LL 1*LL  7* 1 

IFIILL  7-7I.LE.IINCEI-NZI  MGO  TO  536 
115  KWILL6IM 
LL6-LL6*  l 

IFIILL6-6I.LE. (INCEI-NZII) GO  TO  535 
12 C KMILL5I-1 
LL  5*LL5* 1 

IF  II LL5-5I.LE. I INCE I-NZ l I IGO  TO  530 
T*C  KWILL*I-1 
LL*«LL*»1 

IFII  LL*-*I.LE. I INCEI-NZI II  GO  TO  520 
160  KMILL3I-1 
LL3*LL3*1 

IF  ||  LL3-3I.LE.  I INCEI-NZ I IIGO  TO  510 
180  KUILL2I-1 
LL  2»LL  2+ 1 

IF  (I  LL2-2I .LE.  I INCEl-VZ 1 • I GO  TO  500 
800  KMILL1 1*1 
LL1-LL1M 

IFII  LLl-ll.LE. lINCEt-NZIllGO  TO  *90 
820  CONTINUE 
NZ l-NZ 1*1 
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C 

C24  ABORT  FLOW  CONTROL 

C 

IFINABORT.EQ.OIGO  TO  835 
GO  TO  840 
835  CONTINUE 

IF  I IEVENT.EQ.il  GO  TO  970 
C 

C25  RESET  PRIOR  STATE  IIT»  ZERO  LOCATIONS 

C 

840  IF  (NZ.LT.l.OR.NZ.GT.BI  GO  TO  950 

GO  TO  <930. 910. 890, 870. 850.845,844, 8431 . NZ 

843  ITIII8I-1 
1 1 8-  1 1 8*  I 

IF|| 1 18-81. LE. < INCEI-NZ »»G0  TO  407 

844  IT<II7I>1 
1 1 7»  1 1 74 1 

IFK  I I 7-71. LE. I INCEI-NZIIGO  TO  406 

845  I T < I 1 6 I > 1 
I I 6-  1 1 64 1 

IFK  I 16-61. LE.  UNCEl-NZ  > IGO  TO  405 
850  ITIII5I-1 
115-11541 

IFK  II5-5I.LE.  ( INCEI-NZ ) IGO  TO  400 
870  ITII 141-1 

114-11441 

IFK  I 14-4). LE.  IINCEI-NZI)  GO  TO  390 
890  I Til I3I-1 
in-  1 1 34  1 

IFK  113-31.  LE.  I INCEI-NZII  GO  TO  380 
9 LO  ITIII2I-1 

1 1 2- 1 1 24 1 

IFK  I 12-21. LE. I INCEI-NZII  GO  TO  370 
S30  ITIim-l 
II1-II  141 

IFK  II1-II.  LE.IINCEI-NZ1IGO  TO  360 
S50  CONTINUE 

NZ-N241 

IF  < NZ.LE.N2Tl  GO  TO  320 
NZ-0 

S 70  CONTINUE 
C 

C 26  WRITE  OUTPUT  STATE  PROBABILITY  TAPE 

C 

IFINABORT.NE.OIGO  TO  979 
JE  VENT-IEVENT 

WRITE! 3) JEVENT . JCOUNT, I PJILI ,L-1 .JCOUNT ) 

C 

C 27  COMPUTE  THE  MOOE  SEQUENCE  PROBABILITIES 

C 

S 79  READ  (4,12801  ( l JN( I JK » , I JK-l , JCOUNT I 
JCAP-0 

00  980  I JK- 1 .JCOUNT 
J-l JN( l JKI 

PC  API J»-PCAP(J»4PJ< IJKI 
JCAP-NAXOI JCAP.Jl 
S8C  CONTINUE 

PARI VE »PCAPIl»4PCAP(2»  4PCAPI3I 
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C 

C2S  R E AO  ANO  WRITE  OUTPUT  MODE  PROBABILITIES  ANO  C ARABIC  I TE  S 

C 

00  1000  J-l.JCAP 

RE  AO  IS, 12601  CEPI  JI,FA(  J),PK(  J) 

IF  INC.EQ.OI  GO  TO  990 

IFINC.NE.7I  WRITE  (6.12T0I  J,  PC  API  J I , CEP  I J I , FA  ( J I 
IFINC.EQ.7I  WRI TEI6.12 751  J.PCAPIJI 
GO  TO  1000 

990  WRITE  (6,1190)  J.PCAPI  J I ,PM  J) 

1C00  CONTINUE 

IFINC.NE.7I  GO  TO  9030 
WRITE  I 6, 9020 1 PAR  I VE 
9C20  FORMAT!  IX, 'PARI  VE  • NE12.9I 
C 

C29  IF  EVENT  IS  OFFENSIVE .COMPUTE  ET  - MISSION  EFFECTIVENESS 

C 

903C  11-0 

IF  (NC.NE.O)  GO  TO  10)0 
00  1010  J-l.INCEI 
IV! Jl-J 
1310  CONTINUE 
SUM-0.0 

00  1020  J-l.JCAR 
SUM-SUM* PCAPI J)*PKI J» 

1020  CONTINUE 

WRITE  16,12001  (EVENT, SUM 
ET-E  T*  SUM 
WRITE  (6,1210)  ET 
GO  TO  80 


IF  EVENT  IS  OEFENSIVE  DETERMINE  SYSTEM  SURVIVABILITY  PAR  AMETFRS 


10)0  IFINC.EQ.7)  GO  TO  2900 

I 1-1 

00  1040  IQ-1, INCEI 
IV(IQ)-miQ) 

1040  CONTINUE 

IF  (NC.NE.6)  GO  TO  1070 
11-6 
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C 

C 3 1 QUICK  CONVENTIONAL  KILLS 

c 

RE  AO  ($.12901  LMAX.KMAX.IRIKI.K-l.KHAX) 

WRITE  (6.13001  (RIKI.K-l.KNAX) 

00  1050  M-l.INCEl 
00  1050  L-l.LNAX 

RE  AO  (5.1310)  (PCKILLI M.L.KI.K-l.KMAX) 

1050  WRITE  16.1320)  M.L . (PI K ILL  I M, L ,K ) , K- 1 ,KN AX ) 
XMAX1*KNAX-1 
00  1060  I JKM  . JCOUNT 
PMISSI KHAX. I JKI-1.0 
J-I JNI IJK) 

signa-cepi  Ji/i.ua 

OE  NUN»2.0*SIGMA**2 
DENI )5-l T5.0*0ENUM 
EX1-1.0 

00  1060  K-1.KNAX1 

EX2-0.0 

RM2>R(K)**2 

IF  (RN2.LT. OFN175)  EX2*EXP( -RN2/0ENUN) 

PMISSI  K.I JKI -EX1-EX2 

PNISS(KNAX,I JKI-PMI SSI KMAX. I JK)-PMISS(K, IJK) 
EX1-EX2 
1 C 60  CONTINUE 
GO  TO  80 

32  SLOW  KILLS 

2500  1 1 - T 

00  2122  l-l.NLTH 

RE  AD (5. 2 120)  F TFNNI l ) , ( FTONNI l , J ) , J- l.NLTH ) 
2120  FORMAT  I 17F3.2) 

IFF- I- l 

IFil.EQ.il  GO  TO  2125 
PTFN(I)  ■ ( F TFNNI 1 I -F TFNNI IFF) ) 

GO  TO  2129 

2125  PTFNI I )-F TFNNI I ) 

2129  CONTINUE 

00  2130  J-l.NLTH 
JFF-J-1 

IFU.EQ.IIGO  TO  2160 

PTONM  » J )*F  TONNI  I , J) -F TONNI I , JFF I 

GO  TO  2130 

2160  P TON  I I * J ) - FTONNI I , J) 

213C  CONTINUE 
2122  CONTINUE 
GO  TO  80 
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QUICK  NUCLEAR  KILLS 


IC7C  CONTINUE 

OQ  L 090  IJK-UJCOUNT 
QO(t JKl-0.0 
j« I JNI  I JK  ) 

SIGMA*CEPt Jl/l.178 

OENUM*2.0*SIGNA**2 

0ENIT5»175.0*DENUM 

QOI I jk»-exi**ea«J> 

00  1080  M-l.INCEl 
QPRMIM.t JKI -0.0 

RNM2»RM|M)**2 

RNM12«RMIM»ll**2  rF  DFML79I  GO  TO  10S0 

IF  IRMM2.LE.RMN12.0R.RMM12.GE.DEN1I9I  v.u 

E XMl*E  XP I -R MMl  2 /OE NUMI 

a.  * Nii5 . * 

1080  CONTINUE 
IC90  CONTINUE 
GO  TO  80 
1100  REMIND  * 

REW1N0  3 
STOP 


3<« 


DEFINE  FORMATS 


( 5U“U  iuu  III 

format  1 13.  IX, F*. 2 . IX.  I , 3M/5X.  15HEVENT  OCCURRED 

?i«S;5I;!5;h£ve»t  o*,c.i»t«»  .5  .«»  /«•“ 

H90^0RMAT  ;;»?‘\;j,^E;;JJ;?5E^Js*FoJ’oFFENsi5E2ECENT  NUMBER,  13.  IX, 
120C  FORMAT  ( IHO,  10X*W>HEFFfeCHXcnca.> 

1„J*?Si2i6T,M«..««.»HCU;uUTm  .'NCTI'MII  ,S  .IU.T. 

U!0  fO«K»T 

I’lHO  //»*.»««» 


1230 

1260 

1230 

1260 

1270 


FORMAT 

format 

FORMAT 

FORMAT 

FORMAT 


I IHO , 10X  »3H J 


12X,7HFAIJ»  »,E12.5I 


% 
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127?  FORMAT ( 1H0.10X.3HJ  », I 3,2X .9MPCAP I J)  -.E12.5I 

1280  FORMAT  125131 

1290  FORMAT  42t3.10F7.0l 

1300  FORMAT  4 1H0.2X .31HCOMPONFNT  PROBABILITIES  OF  Kill/ 
1 17H  COMP#  ELEV#  R-.10IF8.0, 1X1) 

1310  FORMAT  (8E10.4) 

1320  FORMAT  I2I6.2I4X.E10.AI I 
END 


SUBROUTINE  TWOONZ  1 XG,  N 1 2, 2 1 .NX , NZ  , ANSI  1 

1 TWO  OIMENSIONAl  LINEAR  INTERPOLATION  ROUTINF 

DIMENSION  21(11.  ANSltll 
IXl-INU-ll *NX*1 

C?  IF  GIVEN  VALUE  IS  LESS  THAN  LONER  LIMITS  OF  TABLE,  SET  RESULTS 

C EQUAL  TO  LONER  LIMIT 

IF  (XG-Zlllxm  70.70,10 
10  LL-tXL*l 
LU-I  XL  ♦•NX-1 

C3  SEARCH  FOR  INTERVAL  IN  WHICH  GIVEN  VALUE  LIES 

00  20  J-ll.LU 
IF  (XG-2KJ))  30,80,20 
20  CONTINUE 

C'*  IF  GIVEN  VALUE  IS  GREATER  THAN  UPPER  LIMIT  OF  TABLE,  SET 

C RESULTS  EQUAL  TO  UPPER  LIMIT 

J* I XL*NX -1 
GO  TO  80 

C&  GIVEN  VALUE  WITHIN  TABLE,  CALCULATE  INTERPOLATION  FACTOR 

30  RAT-I2U  JI-XGI /(  21 1 Jl  - 2 1 4 J- 1 1 1 
40  JP-J-INI2-1 »*NX 

C^  00  LINEAR  INTERPOLATION  FOR  RESULTS 

00  60  K«l,NZ 

C7  CHECK  FOR  ZERO  SUBSCRIPT 

IF  I JP.EQ.l 1 GO  TO  50 

AN  SI  I K 1 *2 1 1 JP) -RAT* 1 21  I JPI-21 I JP-1 II 
GO  TO  60 
50  ANSI  1 1 1 -2 1 1 1 1 
60  JP»JP*NX 
RETURN 
70  J-IXL 


SC  RAT-0.0 
GO  TO  40 
ENO 
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ELECTRONICS  MODE  TEST  DECK 


2 3 0 0 

1 ELECTRONICS  A 13. 

23.  10. 

2 ELcCTRONICS  R 20. 

10.  20. 

<599 

1 .10  6 23MM  QUAD,  POSITION  A T 

1112  URINARY  FUNCTION  a 

1221  NORMAL  mode 
l 

1 3 5 5 COMPLETE  FAILURE 

0 99  0 0 

2122  MILITARY  FUNCTION  B 

22  55  NORMAL  MODE 

2 

2355  OFGRAOFO  MODE 

0 59  0 0 

999  300 

2 .20  6 23MM  QUAD,  POSITION  5 T 

1112  MILITARY  FUNCTION  A 

1221  NORMAL  MODE 

1 

1355  COMPLETE  FAILURE 
0 99  0 0 

2122  MILITARY  FUNCTION  B 
2255  NORMAL  MODE 


OOF 
OOF 
C OF 

o 0 F 
OOF 
OOF 


OOF 
OOF 
0 OF 

OOF 

OOF 


2 

2 3 

0 99 

999  0 

3 

1 l 

1 2 
1 

1 3 
r 99 

2 l 

2 2 
2 

2 3 

0 99 

999  0 

1 l 

1 2 
1 

1 3 
P 99 

2 l 
2 2 
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•5  5 DFGRAOEO  MODF 

3 0 

0 0 

.30  6 23MM  QUAD,  POST  T I ON  4 T 

1 2 MILITARY  pUNr  T I ON  A 

2 1 NORMAL  MODE 

S 5 COMPLFTF  FAILU*F 

0 0 

2 2 MILITARY  FUNCTION  B 

5 5 NORMAL  MODF 

5 * DEGR APC  P MfjOF 

3 C 

0 0 

.30  n MARK  82  SNAKFYt  F 

1 2 MILITARY  FUNCTION  A 

2 1 NORMAL  MODF 

5 5 CDMPLFTE  F A I L 'J*  F 

0 0 

2 2 MILITARY  FUNCTION  B 

5 5 NORMAL  MODE 


0 0.  F 

OOF 
OOF 

OOF 
0 OF 
OOF 


OOF 

OOF 

OOF 

OOF 

OOF 


99 


— ■ — 
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2 

2 3 5 5 OEGRADEO  MODE 

n 99  ^ f» 

999  300 

5 .43  A 23MM  QUAD*  POSITION  4 T 

1112  MILITARY  FUNCTION  A 

1221  NORMAL  MOOE 
1 

1 3 5 5 COMPLETE  FAILURE 

0 99  O 0 

2122  MILITARY  FUNCTION  B 

2255  NORMAL  MODE 

2 

2355  DEGRADED  MOOE 
0 99  3 0 

999  000 

6 .50  0 LAND  AT  BASE  T 

1112  MILITARY  function  A 

1255  NORMAL  MODE 

0 99  0 0 

999  000 


OOF 

OOF 

OOF 

OOF 
OOF 
5 5 T 


0 OP 
OOF 


100 


.1 

J 
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2 

0 0 0 

0 0 

1 

ELEC  TRONIC  S A 

10. 

? 

ELEC  TRONIC  S ft 

20. 

999 

1 

.10  6 

23MM 

QUAD  » °n$I T l ON  4 

T 

10. 

1. 

10. 

1. 

10. 

1 . 

l 

2 100. 

0.  0. 

.1  0. 

l O00. 

2 

. 20  f> 

23MM 

QUAD, POSITION  4 

T 

10. 

1. 

10. 

1. 

10. 

1. 

1 

2 100. 

0.  0. 

.1  0. 

1000. 

1 

. 30  6 

23«M 

QUAD , POSI T I ON  4 

T 

10. 

1. 

1 0. 

1. 

10. 

1. 

1 

2 100. 

0.  0. 

.1  0. 

1000. 

4 

.30  0 

MARK 

82  SNAKEYE 

F 

o.o 

0.0 

0.6 

5 

.40  6 

23MM 

QUAD, POSIT  ION  4 

T 

10. 

1. 

10. 

1. 

10. 

1. 

1 

2 l OC. 

C.  0. 

. 1 0. 

100". 

6 

.40  " 

LA  sr 

AT  ft A$c 

T 

0.3 
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VEHICLE  MODE  TEST  DECK  - PROGRAM  I 


I 1 

ACTUALLY  GOOD  POP  2 OP  3 IT 


ACTUALLY  GOOD  FOR  DELTA  T 


APPARFNTLY  GOOD  FOR  2DR3  DT 


APPARFi^TLY  GOOD  FOR  DELTA  T 


7 2 3MM  OUAD,  POSITION  A 
FLIGHT  FUNCTION 


10000. 


1 3000 . 


10000. 


10000. 


5 5 

0 0 

0 0 

,20 

l 2 


5 5 


5 5 


5 5 


5 5 


NORMAL  MODE 


7 2 3MM  QUAD,  POSITION  A 
FLIGHT  FUNCTION 

NORMAL  MODE  A 

NORMAL  MODE  B 

NORMAL  MODE  C 


ABORT  MODE  A 


ABORT  MODE  B 


DOWN 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OF 


OOF 


5 * 


•>  •> 


r 99  o o 


999  3 o 0 
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7 23MM  QUAD v POSITION  A 
PLIGHT  FUNCTION 

NORMAL  MOOF  A 

NORMAL  MODE  R 

NORMAL  MOOE  C 

ABORT  MODE  A 

OnwN 


1 ’ 
* 

l 3 
1 

l * 


S 5 


5 5 


*»  S 


O 0 


999  000 


7 2 3MM  QUAD,  POSI  T l ON  A 
FLIGHT  FUNCTION 

NORMAL  MODE  A 

NORMAL  MODE  B 

NORMAL  MODE  C 


DOWN 


6 

l l 


, 5 D 
l 2 


S S 


5 5 


S •» 


7 LAND  AT  BASF 
FLIGHT  FUNCTION 

NORMAL  MODc  A 


normal  MOOF  B 


NORMAL  MODE  C 


DOWN 


0 99  0 0 


999  300 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 


OOF 
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VEHICLE  MODE  TESr  DECK  - PROGRAM  2 


4 1114  0 

l VFH  ACTUALLY  SHOO  FOR  2 0*  3 OT 


10000. 


2 VFH  ACTUALLY  0003  FOR  OFL  T A T 


iOOOO, 


3 VFH  APPARENTLY  GOOD  FOR  20R3  OT 


10000. 


4 VFH  APPARENTLY  5000  FOR  3ELTA  T 


10000, 


1 . 10  7 2 3 MM  QUAD,  POSIMON  4 

10.  1. 

02100 100 l 00 l 00 

04  50100100100 

05  30  60  90100 
100  20  40  70100 

2 .20  7 2 3MM  QUAD.  POSITION  4 

10.  1. 

10.  1. 

10.  1. 

10.  1. 

10.  1. 

10.  1. 

02100100100100 

04  50100100100 

05  33  60  90100 
100  20  40  70100 

3 .30  7 2 3 MM  QUAD,  POSITION  4 

10.  1. 

10.  1. 

10.  1. 

10.  1. 

10.  1. 

02100100100100 

04  50100100100 

05  33  60  90100 
100  20  40  70100 
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*5  .40 

10. 

1 D. 

10. 

10. 

02lCDlD01001or* 

04  5010)100100 

05  00  60  90100 
100  20  60  70100 

6 .50  7 LAND  AT  BASF 


7 23Mm  QUAD,  POSITION  4 
1. 

1. 

1. 

1. 


ID. 

10. 

ID. 

ID. 


1. 

1. 

1. 

1. 


021 DD 1001 0O ion 

04  50100100100 

05  3D  60  90100 
100  ?D  4 0 70 ICO 
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SIMULATION  MODEL 


This  section  discusses  the  manner  in  which  calculations  are  performed  within  the 
computer  routines.  The  blocks  of  code  are  headed  by  comment  cards.  These  comment  cards 
are  also  the  titles  of  the  schematic  diagrams  (Figures  21  through  49),  which  facilitate 
cross-referencing.  Following  the  discussion  of  the  code  is  a subsection  entitled  User 
Information,  which  describes  restrictions,  simulation  errors,  and  limitations  affecting  input. 
This  is  followed  by  a description  of  abbreviations  and  symbols  for  the  simulation  model. 


PROGRAM  I 

Program  I consists  of  only  one  routine.  The  purpose  of  this  routine  is  to  define  and 
number  all  possible  mode  sequences  for  all  events  resulting  from  all  possible  states  of  the 
system.  The  mode  sequence  numbers  are  an  array  that  is  stored  on  an  intermediate  device 
for  use  in  Program  2. 

statements 


MISSION  DAMAGE  EFFECTIVENESS  MODEL  - PROGRAM  ONE 

MISDEM  - PGM  l 


DIMENSION  CN(3I  .CNAMEt  8,23>,DEF(3) ,DNAMEI8),FNAME( 10,27, 10), 

1 H)MT(A) ,JES(  10,2  7,25  ) , JUt  23 , 25 ) , KS 1 10 . 27 » , KU < 23  ) ,KM( 23) . LEEt 10) . 

2 IF  1 101 ,LLF(2  7I .LL0I13 ,27) ,LMA( 10,27,23 »,LM AT t 10,  27, 23) .LKATTI  23 ) 

3 , L0( 10,27) ,MD( 10,27)  , MI  I 23 ) , MMQt 1 0, 27 ) ,MQ t 10 , 27 ) , OFF ( 3 ) , MEAPNI 7 ) 
INTEGERS  I JNt  256 ) 

LOGICAL  M0, MOT 

EQUI VALENCE  tLMAt 1 ,1,1 ) , LMAT (1,1,1) ) 

DATA  BLK/AH  / , OFF  / A HOFF  E , AHNS  I V , AHE  / .DEF/AHDEFE , AMNS IV, 

•AHE  /, OUT /2HN0/ 

are  used  to  allocate  storage  by  use  of  the  DIMENSION  statement,  to  reserve  two  integer 
storage  locations  for  each  of  the  256  INTE(1ER*2  words,  to  declare  the  variables  Ml)  and 
MDT  as  logical  with  four  storage  locations  allocated  for  each,  to  define  storage  that  is  to  be 
shared  by  two  or  more  entities  by  use  of  the  EQUIVALENCE  statement,  and  to  define  the 
initial  values  for  various  variables  by  use  of  the  DATA  statement. 

The  statements 
C 

Cl  INITIALIZE 

C 

LI  *2 7 
LCN8“8 
L A0* 10 
LCEI  -23 
L 72“  7 


C 
C 
C 
C 

c 
c 

I 
■ 
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are  used  to  define  variable  names  for  the  dimensional  constant.  LI  is  the  maximum  number 
of  modes  (in  all  subfunctions)  of  the  mode  logic  (see  Figure  21).  L40  is  the  maximum 
number  of  subfunctions.  LCN8  is  the  maximum  number  of  equipment  name  segments.  L72 
is  the  maximum  number  of  mode  sequences  (or  flow  paths).  LCEI  is  the  maximum  number 
of  subsystems  in  the  equipment  configuration,  used  as  a check  (later  in  the  program)  on 
program  input. 

The  statements 

1NCE  1*0 

00  20  l-ULCEt 
MIC!  )»l 

00  20  J-1.LCN8 
CNAMEt  J,l  >«RIK 
20  CONTINUE 

are  used  to  initialize  subsystems  counter  (INCEI)  to  zero,  to  set  the  on-off  (lag  to  "on”  and 
to  insert  blanks  in  all  the  subsystem  name  locations. 

The  statement 

RE  AD  I 5 ,950) N2T .MCR.MPP , MAY 


is  used  to  read  the  first  input  data  card.  These  variables  arc  used  for  control  of  state 
generation  (through  NZT)  and  state  probability  generation  (through  MCR.  MPR  and  MAV). 
These  tour  inputs  are  discussed  at  length  in  the  User  Information  subsection. 

The  statements 


C 

C?  READ  AND  WRITE  SYSTEM  DESCRIPTION 

C 

WRITE! 6,9601 
WR I TE ( 6, 970 ) 

00  A0  !■  1, LCEI 

READ (5,980)  M, (ONAMEtJ ) , J« 1 , 1C N8 ) • T BFM, ( HOMT ( M PRM ) , MPRM» 1 , A ) ,0N0F F 
IF(M .E  Q.999IG0  TO  60 

WR  |TE(  6,990)  M,  (DNAMEt  J ) , J*1  , LC N8  ) , T BFM,  ( l-DMT  < NPRM  ) , MPRM.  1 , A » , ONOF  F 

IFtM.GT.LCEI.OR.M.LT.U  GO  TO  50 

IFIM.GT. INCEI) INCEI-M 

IF (0N0FF .EQ.0UT )Ml 1 1)»0 

DO  AO  J*1,LCN8 

CNAMEt  J,  l )« DNAMEt  J> 

AO  CONTINUE 

RE  AO  t5,980>  M 
IFIM.E0.999)  30  TO  60 
50  MRI TE  t 6, 1000) 

STOP 

60  CONTINUE 
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are  used  for  input/output  of  the  system  description.  The  first  two  lines  write  the  headings  in 
preparation  for  the  output.  The  next  statement  starts  a loop  which  allows  for  the  maximum 
input  of  23  subsystem  definitions.  Each  subsystem  definition  consists  of  the  equipment 
name,  mean  time  between  failures,  nuclear  damage  thresholds'  and  a logical  flag  (ONOFF) 
which  can  be  used  to  suppress  a subsystem  (prevent  its  being  considered  further).  If  the 
system  number  that  is  input  is  equal  to  999,  the  program  switches  control  out  of  the  loop 
(i.e„  “60  CONTINUE")  signifying  that  there  are  no  more  inputs  of  this  type.  However,  if  a 
system  configuration  is  input,  the  program  proceeds  to  write  out  information.  The  next 
statement  is  used  to  check  for  blanks  (subsystems  left  out)  or  too  many  subsystems  (relative 
to  LCEI)  in  the  input  data.  If  either  one  of  these  conditions  exist,  program  control  is 
transferred  to  an  error  message  which  is  written  out  and  the  program  stops.  The  next 
statement  updates  the  counter  which  keeps  track  of  the  total  number  of  subsystems  input. 
The  next  statement  is  used  to  set  a flag  if  a subsystem  is  being  suppressed  The  next  three 
statements  store  the  subsystem  name  by  a reference  index  number  corresponding  to  LCEI. 
Program  control  is  then  transferred  to  the  beginning  of  the  system  configuration  loop  to 
read  the  next  subsystem  description.  Should  the  program  go  through  the  loop  23  times  and 
never  transfer  control  to  the  statement  which  says  "60  CONTINUE”,  the  program  will  read 
another  card.  If  this  card  contains  another  system  description,  the  program  will  write  an 
error  message  and  stop;  however,  if  the  input  contains  a “999"  for  the  subsystem  number, 
control  will  be  transferred  to  the  “60  CONTINUE"  statement  and  program  execution  will 
flow  as  described  in  the  next  set  of  statements. 

The  statement 

REWIND  * 

is  used  to  rewind  the  tape  to  the  beginning  on  which  the  mode  sequence  number  array  is  to 
be  written. 

The  statement 
70  CONTINUE 

begins  the  event  description  loop. 

The  statements 


C 

C3  READ  AND  WRITE  EVENT  DESCRIPTION 

C 

READ(5,1010,END«8*0I  I E VENT , T2 ,NC, (WEAPNI I » , !« I , 7 I ,MDT 
WRITE! 6» 1270) 

00  71  N-1,3 
CNIN»«0FE(N) 

IEINC.NE.0)  CNINI *0EF( NI 
71  CONTINUE 

WR  ITEI6, 1020)  IE  VENT  * I CNt  NI »N*1 .31.  T 2,  I WEAPNI I) , I « 1 , 7 ) , MDT 


J 
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arc  used  to  read  and  write  the  event  description  (one  of  many);  event  number,  time,  type 
(defensive  or  offensive),  weapons  description,  and  mission  descriptor  defining  the 
environmental  conditions  at  the  time  of  the  event. 

I F « NC . NE . 7 1 MRITEt6, 1)601 
tFtNC.E0.7l  WRITEt6, 1065) 


are  used  to  select  between  mode  logic  output  table  headings  corresponding  to  electronics 
modes  (NC  = 0 to  6)  or  vehicle  mode  (NC  = 7). 

The  statements 


C 

C4  REAO  SU8FUNCT 10NAL  FLOW  FOR  EVENT  ANO  GENERATE  SUBSYSTEM 

C REQUIREMENTS  VECTOR 

C 

00  250  1*1.140 
00  240  K*1,L1 
M0( I ,KI*. FALSE. 

REAP  15 ,10701 LF ( I I , LLFt K » ,1 0(1 , K I , L LQt 1 , K I , ( FNAME  t I , K , J I , J*  1 , 10  I , 
* MQtl ,K) ,MMQ< I , K ) ,M0t I.K), (LMAT ( I.K.JI »J*1»  INCE I » 

L * 0 

00  90  J* l, INCE I 
IFtLMATt I ,K, J» .EQ.01G0  TO  100 
80  L-L*l 

IFtLMATt 1 ,K,J» .EQ.LIG3  TO  90 
LMAT  T t LI *0 
GO  TO  80 
90  LMATTt  LI  * 1 
GO  TO  140 
100  J*L*l 

IF  ( J .G T.  I NCE  1 1 GO  TO  1<>0 
00  110  L-J.INCEI 
110  LMATTt  L I *0 
140  00  150  L*l,INCEl 

150  LMAJI.K.L)  -LMATTtLI 


are  used  to  read  the  subsystem  requirements,  specified  in  terms  of  subsystem  ordinal 
numbers,  and  then  to  derive  the  subsystem  requirements  binary  vector  for  each  subfunction 
and  mode,  and  store  it  in  the  LMA  array. 


The  statements 


C5  WRITE  SUBFUNCTIONAL  FLOW 

C 

IFtNC.NE.7l 

1 WRITE  (6,  1080  ILF  I 1 1 , LLC  f K) ,LQI I .KI.LLQt  I.K),  t FNAME  t I.K.JI. <1*1.  101. 

* MQtl ,K| ,MMQt I.Kt ,M0t I.KI.ILMATI I.K, JI.J-l, INCE 1 1 
IFtNC.EQ.7l  MR  I TE  1 6, 10851  LF til , LL F IKI , LQt I , K) , LLQt I , K I , 

• IFNAMEtl.K.JI.J- 1,131 .tLMATt I.K, J I, J* l, INCE I I 
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IFILLFUI  .E  Q.  99  ) GO  TO  230 
IFILFtll  .EQ.999I  GO  TO  260 
GO  TO  2*0 
230  LEE  ( I ) *K-1 
GO  TO  250 
2*0  CONTINUE 
250  CONTINUE 
I«L*0*l 
RE ADt 5 . 1 070 ) H 
IFIN.EQ.999)  GO  TO  260 
WRITE! 6« 1090) 

STOP 


are  used  to  write  the  input  subsystem  requirements  in  binary  form,  along  with  subfunction, 
mode  and  mode  logic,  as  shown  in  Figure  23. 


The  process  of  reading,  generating  the  LMA  vector,  and  writing  as  discussed  in  the 
previous  two  paragraphs  continues  until  one  of  two  things  happen.  If  the  event  input  list  has 
been  exhausted,  control  will  be  transferred  to  the  statement  ‘‘260  CONTINUE".  However, 
if  the  number  of  inputs  exceed  the  dimensional  values,  the  program  will  write  an  error 
message  and  stop. 

The  statements 


C 

CS  INITIALIZE  FOR  NODE  SEQUENCE  IDENTIFICATION 

C 

260  LE-I-i 
LQSVL-0 
LLQS-0 
NZ*0 

DO  290  J * I » L *0 
00  290  K*|,Ll 
KSt J.KI-0 
DO  290  I * L t L 72 
JES<  J,K, I )«0 
290  CONTINUE 

DO  300  l -l.LCEI 
KWt I 1*0 
300  CONTINUE 

00  310  JCOUNT *1«256 
IJNI JCOUNT 1*0 
313  CONTINUE 
JC OUNT  »0 
JC  AP  *0 


are  used  to  initialize  all  variables  used  in  the  mode  sequence  identification  process  which 
starts  in  the  next  block  of  code. 
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The  statements 

C 

C 7 GENERATE  FINAL  STATE  OF  THE  SYSTEM  - BINARY  VECTOR  KW 

C 

330  CONTINUE 

00  340  1*1*1 NC El 
KW( I 1*1 
340  CONTINUE 

IFINZ.LT.1IG0  TO  400 
LL1-1 

350  KWILL1 1*0 

IFtNZ.LT.2IG0  TO  400 
LL2*LL1*1 
360  KMILL2 1*0 

(F (NZ.LT .31  GO  TO  400 
LL3-LL2M 
370  KM(LL3I*0 

IF (NZ.LT *4)G0  TO  400 
LL4«LL3*1 
380  KW(LL4I*0 

tF(NZ.LT.5IG0  TO  400 
LL5-LL4M 
340  KWILL5I-0 

IF (NZ.LT .61  GO  TO  400 
LL6-LL5+1 

341  KM(LL6 1*0 

IFtNZ.LT. 7)  GO  TO  400 
LL7*LL6«-1 

392  KWtL  L 7 1*0 
IFtNZ.LT. SI  GO  TO  400 
LL8=LL7»1 

393  KWt  L L8 1*0 
400  CONTINUE 

work  in  conjunction  with  a later  block  of  code  to  generate  a set  of  current  state  (binary) 
vectors  KW,  limited  by  the  current  number  of  zeros  (NZ)  in  the  array. 


The  statements 

C 

C 8 TEST  FOR  SUPPRESSED  SUBSYSTEMS 

C 

00  410  MM*l * I NCE I 
I F ( M I ( MM  I . NE . 0 I GO  TO  VIO 
IFtKMtMMI.EQ.llGO  TO  720 
410  CONTINUE 

tF(NC.EQ.7IG0  TO  405 
IFtMCR.EO.OIGO  TO  402 
IF(KWtMCR) . EQ.OIGO  TO  720 
402  IFtMPR .EQ.OIGO  TO  404 

IF(KMtMPR).  EQ.OIGO  TO  720 

404  IFtMAV. EQ.OIGO  TO  405 
IFtKWtMAVI. EQ.OIGO  TO  720 

405  CONTINUE 
JCOUNT  * JCOUNT 
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are  used  to  reject  (i.e.,  to  withhold  from  further  use)  a state  if  it  requires  either: 

1.  A subsystem  be  viable  that  has  been  deleted  (by  the  Ml  flag) 

2.  A subsystem  (one  of  three)  be  failed  that  is  critical  to  all  significant  modes  (the 
subsystems  are  labeled  MCR,  MPR,  and  MAV). 

The  second  test  applies  only  to  the  electronic  mode  (NC  = 0 to  6).  Note  that  no  number  is 
assigned  to  a rejected  state. 

The  statements 


C 

C 9 DEFINE  MODE  SEQUENCE  AND  SUBSYSTEMS  USED 

C 

DO  530  I -1. ICE  I 
53C  KUIll-0 
LTEST-l 

DO  640  l * l.LE 
LEQ-LEEIL) 

DO  540  K«  1 ,LEQ 
540  KSfL.KI-0 

IF(L. EQ.LQSVL.OR.LTEST.EQ.il  CO  TO  550 
CO  TO  640 

550  00  630  K-l.LEQ 

IFtK.GE.LLQS. OR. LTEST. EQ.l » GO  TO  560 
GO  TO  630 

560  (FI. NOT. MOT. AND. MO (LtKI ) GO  TO  580 
CD  TO  590 
58C  LQSVL-MQtL.KI 
LLQS-MMQU.KI 
LTEST-0 

IF(LQSVL.EQ.L)  GO  TO  630 
GO  TO  640 
590  LTEST-l 

DO  600  N-l.INCEl 

IF(KM(M) .LT.LMAtL.K.MI | GO  TO  620 
600  CONTINUE 

DO  610  N-l.INCEI 

IFtKM(M) .GT.LNAf L.K.MI  .OR. KM  INI . EQ.O ) GO  TO  610 
IFlKWtMI.EQ.LNAIL.K.N)  .AND.KMIMI.EO.il  KUIMI-l 
610  CONTINUE 
KSIL.KI-l 
LQSVL-LQt L.KI 
LLQS-LLQ(L.K) 

LTEST-0 

620  IFILOSVL.GT.LIGO  TO  649 
630  CONTINUE 
640  CONTINUE 
J-JCAP+1 
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are  used  to  generate  the  KS  and  KU  binary  arrays  which  define  subfunction/mode  numbers 
used  and  the  subsystem  numbers  used,  respectively,  by  the  current  state.  LE  and  LEE  are 
the  number  of  subfunctions  and  modes  (per  subfunction)  defined,  which  limit  the  outer 
(subfunction)  and  inner  (mode)  do-loops  ending  at  640  and  630,  respectively.  LQSVL  and 
LLQS  are  the  subfunction  and  mode  numbers  which  are  assigned  for  subsequent  testing  of 
the  KW.  LQSVL  and  LLQS  are  initially  1 and  1.  They  are  later  assigned  the  values  MQ, 
MMQ  if  the  mission  descriptor  (MD)  test  fails,  and  are  assigned  the  values  LQ,  LLQ  if  the 
LMA  test  (subsystem  requirements)  is  successful.  The  LMA  test  takes  place  in  the  inner 
do-loop  on  M,  ending  at  600.  The  KU  array  and  KS  array  are  loaded  with  a I after  each 
successful  mode.  LTEST=0  is  a flag  denoting  failure  of  the  MD  test,  which  results  in 
rejection  of  further  tests  on  the  current  subfunction  and  mode.  When  LMA  test  fails,  the 
very  next  mode  is  required,  so  LTEST  is  set  to  1,  preventing  the  selection  of  the  next 
subfunction.  As  a result  of  these  processes,  the  KU  (M)  array  is  the  union  of  all  subsystem 
requirements  of  the  mode  sequence. 

The  statements 


C 

CIO  AS 

c 

*50  J*J-1 


ASSIGN  NODE  SEQUENCE  NUMBERS 


IE  < J.NE.0I  GO  TO  680 
JCAP-JCAP*! 

J*  JC  AP 

DO  660  L - 1 «LE 
LEQ-LEECL) 

DO  660  K*1,LEQ 
JESt L.K.JI-KSI  L.K) 

660  CONTINUE 

DO  670  I U*1 • I NCE I 
JU( IU. JI*KUI IU) 

670  CONTINUE 
GO  TO  700 

680  DO  690  L ■ l.LE 
LEQ-LEEtl) 

00  690  K«I,LEQ 

IF< JEStl.Kt JI.NE.KSI L.K) ) GO  TO  650 
690  CONTINUE 
7C0  CONTINUE 

UNI  JCOUNT  | >J 

are  used  to  assign  a mode  sequence  number  (J)  to  the  current  state  number  (JCOUNT).  This 
is  done  by  comparing  the  new  JS  array  against  its  immediate  predecessor,  identified  by  the 
variable  JES.  If  identical,  the  decremented  J number  is  assigned  to  the  current  JCOUNT.  If 
not,  the  next  earlier  KS  is  tested  (i.e.,  J is  decremented  again).  This  is  continued,  if  no 
match  is  found,  to  the  first  KS.  If  it  still  is  not  matched,  the  tentative  JCAP  is  adopted  as 
the  J value.  Thus,  in  the  testing  against  all  previous  KS,  any  match  stops  the  process  at  some 
J,  which  then  gets  ascribed  to  the  current  JCOUNT.  along  with  the  associated  JES  and  JU 
(the  final  array  representing  subsystems  used).  But  if  it  fails,  the  JCAP  value,  which  has  been 


113 


JTCG/AS-76-S-004 


saved  as  a potential  value  of  J from  the  beginning  of  the  test,  is  finally  authorized. 
Incidentally,  JCAP  tracks  the  largest  value  of  J,  and  is  used  as  a limit  on  do-loops  in 
subsequent  portions  of  the  program. 


The  statements 


CU  REPOSITION  ZEROS  IN  STATE  BINARV  VECTOR  (KWI 
C 

720  Ic(NZ.LT.i.0R.NZ.GT.8)G0  TO  790 

GO  TO  (770, 760, 750, TAJ, 730, 728, 726, 726),  NZ 
726  KW(LL8I>1 
LLS-LL8*! 

IF((LL8-8).LE. (INCEI-NZ))  GO  TO  393 
726  KW(LL7)*l 
LL  7*LL 7*1 

IF(( LL7-7I.LE. ( INCEI-NZI)  GO  TO  392 
728  KW(LL6) *1 
LL6*LL6* 1 

IF  ( ( LL6-6I.LE. ( INCEI-NZ 1 1 GO  TO  391 
730  KW(LL5)«1 

LL5-LL5*1 

IF(( LL5-5I.LE. (INCEI-NZI >GO  TO  390 
760  KW(LL4 I *1 
LL A*LLA* 1 

IF(( LLA-AI.LE.  (INCEl-NZ)IGO  TO  380 
750  KW(LL3I*1 
LL3-LL3*1 

I F ( ( LL3-3I.LE. ( INCEI-NZI »G0  TO  370 
760  KW(LL2)~l 
LL2*LL2*1 

IFKLL2-2I.LE.  (INCEI-NZMGO  TO  360 
77C  KW( LL 1 i = l 
LLI»LL1*1 

IF((LL1-1I.LE. (INCEI-NZMGO  TO  350 
790  CONTINUE 
NZ*NZ*1 

IFtNZ.LE.NZTISO  TO  330 
NZ=  0 


are  used  to  increment  the  locations  (LL1,  LL2,  etc.)  of  the  (up  to)  eight  possible  zeros 
within  the  KW  array  after  each  state  has  been  “pushed”  through  the  preceding  three  blocks 
of  code. 


The  statements 


C 

Cl?  WRITE  MODE  SEQUENCE  AND  SUBSYSTEMS  USED 

C 

DO  820  J*l, JCAP 
WRITE! 6, 11701 
WRITE(6, 1160)  J 
00  800  L ■ 1 , LE 
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LEQ*LEE<L ) 

00  800  K*1,LEQ 

I F I JESIL.K, JI.EQ.O)  33  TO  800 
WRITEI6, 1150)  IFNAMEIL.K.I  1,1*1,10) 
800  CONTINUE 

IFINC.NE.7)  MR  ITEI6.il 80) 

IFINC.EQ.7)  WRITEI6.ll 85) 

00  810  W-l.INCEI 

IF!  JU(M,  J),NE.  1)  GO  TO  810 

WRITE  (6.  1190)  (CNAMEU,M),K*1,LCN8) 

810  CONTINUE 
820  CONTINUE 


aiv  used  to  print  the  modes  and  subsystems  associated  with  each  mode  sequence  number. 
This  is  accomplished  by  means  of  the  JES  and  JU  arrays,  associated  with  J- 


The  statements 


C 

Cl  3 WRITE  MODE  SEQUENCE  ON  TAPE  FOR  USE  IN  PROGRAM  2 

C 

WRITE  I A,  1280)  IIJNtlJK) , I JK= l , JCOUNT ) 

GO  TO  70 
8 AO  END  FILE  A 
REWIND  A 
STOP 


write  the  mode  sequence  number  array  on  unit  4.  Control  is  transferred  back  to  statement 
"70  CONTINUE”  where  another  event  is  processed.  When  all  events  have  been  processed, 
the  tape  containing  the  mode  sequence  number  is  rewound  and  Program  1 stops. 


The  statements 


C 

Clt  OEFINE  FORMATS 

C 

950  FORMAT (1013) 

960  FORMAT  I26HI  STSTEM  CONFIGURATION) 

9 70  FORMAT  I IHO, 8X , 9HEQU l PME NT, 25X,  AHMTBF,7X,7HTHOM( 1 >,6X, 7HTHOMI 2), 

• 7X  » 7HTH0M( 3 ) , 8 X , 7HTH0W (A)/56X,3HG  0, 1 IX , 1HN, 12X, 1HB, 15X, IHT//IXI 
980  F0RMAT(I3,1X,8AA,F10.2/AE12.5,T65,A2> 

990  F0RMATI1H  ,1X,13,2X,8AA,F10.2,2X,A(E12.5,2X),3X,A2> 

1000  FORMAT  ( 50H0  TOO  MANT  INPUTS  OR  M IS  OUTSIDE  ALLOWABLE  RANGE) 

1010  FORMAT (I3,1X,F6.2,IX,I3,1X,7AA,2X«L1) 

1020  F0RMATI1HU,5X,9HEVENT  NO. , 1 2 , IX , 3H IS  , 3AA/ 5X, 15HE VENT  OCCURRED  . 

• F 6 • 2 , 

*20H  HOURS  AFTER  TAKEOF F /5X . 2 IHEVENT  DESCRIPTION  IS  ,7AA/5X,L1> 

1060  FORMAT  I 1H0 ,7X .16HSUBFUNCT I ON/ MODE ,32 X, 21HEQUIPMENT  DESCR I PT ION, 1 2 

• X, 18HMI SSI  ON  OESCRIPTOR/1X) 

1065  FORMAT ( IHO,  7X , 1 6HSUBF UNCT I ON/ MODE . 32X , 35HFL 1GHT  AND  DETECTION  TIM 
*E  REMAINING) 

1070  F0RMAT(A13,IX,10AA,1X,2I3,IX,L1/23I3) 

1080  F0RMAT(lx,2I3,2X,2I3,2X,l0AA,T86,2I3,l2X,Ll,T59,231l) 
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1085  FORMAT(IX,2I3,2X,2I3,2X,IOA9,  T59.2311) 

1C90  FORMAT  ( 29H0  TOO  MANY  INPUTS  FOR  F ARRAY) 

1190  FORM AT ( IHO, 1 IX , I 3) 

115C  FORMATUH  t 95X  1 1 0A9 1 

1170  FORM  A T( 1 H ,5X,16HM00E  SEQUENCE  NO* 16X» 18HSUBFUNCT 10NAL  FLOW) 
1180  FORMAT ( 1H0*  58X  * 1 5H SUBS Y STE MS  USED/1X) 

1185  FORMAT (IHO* 5 8X  *31H ACTUAL  AND  APPARENT  FLIGHT  TIME) 

1190  FORMAT I 1H  ,60X,8A9) 

1270  FORM  AT (2  3H 1 EVENT  DESCRIPTION) 

1280  F0RMATI25I3) 

END 


define  the  formats  used  by  Program  1 to  read  and  write  the  inputs  and  outputs. 


PROGRAM  2 

Program  2 consists  of  a main  routine  and  one  subroutine.  The  purpose  of  Program  2 is 
to  use  the  mode  sequences  defined  in  Program  1 in  conjunction  with  input  capabilities  (own 
system  and  threat)  to  compute  the  mode  probabilities  and  system  effectiveness  and 
survivability  parameters. 

Main  Routine 

The  purpose  of  the  main  routine  is  to  confiol  the  inputs  and  perform  all  logic  and 
computations  except  for  an  interpolation  that  is  performed  by  the  subroutine. 

The  statements 


**  ************************  * *******  ****************************  ******** 

MISSION  0AMAGE  EFFECTIVENESS  MOOEL  - PROGRAM  TWO 

MISDEM  - PGM  2 

** **** ********************  ************************************ ******** 
DIMENSION  CEP( 50) .CNI3I , CURVE! 9,30 ) .OEF ( 3> , ONAMEI 8 ) ,FA(50) , 

1 FTONN (16*16) ,FTFNNI 16 ) .GRAPH ( 30 ) .HOMT (9 ) , ISUB( 9 ) , IT< 23) , ITT ( 23 ) 

2 *IV(23),IZ(23)*JZ(23),KW(23), KWW( 23  I » M I( 23)  ,MT6F( 23) ,OFF< 3) , 

3 PCAP<  50) ,PCKI LL( 23,  10 ,27 ) .PCSURV ( 23, 10, 27 ) . P I( 256 ) ,P J( 2 56) , 

9 PKC  50) *PM(23) ,PMM( 2 3 ) .PMISS ( 27, 256) , PTONI 16, 16 ),PTFN( 16), 

5 Q0I256) ,QPRM( 23, 256) ,R( 10) , R A(23 ) ,R€SUL (2 ) ,RM( 29 ) , THDM< 23,9), 

6 TIMEFI23I »TIMEN(23I*  WEAPN( 7 ) 

INTEGER«2  IJNI256) 

LOGICAL  MOT 

OATA  0FF/9H0FFE , 9HNSI V , 9HE  /,0EF/9H0EFE,9HNS IV.9HE  /,0UT/2HN0/ 
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arc  used  to  allocate  storage'  by  use  of  the  dimension  statement  to  reserve  two  integer  stoiage 
locations  for  each  of  the  256  1 NTI-C • I R*2  words,  declare  the  variable  MDT  as  logical  with 
four  storage  locations  allocated  for  each,  and  define  the  initial  values  for  various  variables  by 
use  of  the  data  statements. 

The  statements 
C 

Cl  INITIAlUE 

C 

ICN8-8 
ICE! -23 
t 72-40 


are  used  to  define  variable  names  for  the  dimensional  constants.  L('N8  is  the  maximum 
number  of  equipment  name  segments.  LC'i-l  is  the  maximum  number  of  subsystems  in  the 
equipment  configuration,  used  as  a check  (later  in  the  program)  on  program  input.  1 72  is 
the  maximum  number  of  mode  sequences  generated  by  Program  l for  any  event. 

The  statements 


INCE 1-0 

00  30  l • V *t.C€  l 
MTBE ( I 1-0.0 
P»t I ) -0. 0 
*1 (Il-l 

tvm-i 
mi  >-i 

00  20  MPRM-1,4 
THOMt l .MPRNI -0 
20  CONTINUE 
30  CONTINUE 

12 

are  used  to  initialize  the  subsystems  counter  (INOEI)  to  zero,  to  set  all  MTBF  values  and 
probability  of  reliable  operation  values  equal  to  zero,  to  set  the  on-off  flag  for  all 
subsystems  to  “on”,  to  initialize  the  event  and  subsystem  counter  arrays  and  to  zero-out  all 
damage  threshold  storage  locations  for  all  subsystems. 

The  statement 

RE  AO  (5.11101  N2  T , NCR, NPR . MAV .Ml TH.N ABORT 

is  used  to  read  the  first  input  data  card.  These  variables  are  used  for  control  of  state 
generation  (through  NZT).  state  probability  generation  (through  MCR,  MPR.  and  MAV). 
mission  length  (MLTH),  and  event  abort  control  (N ABORT). 

y 


* 
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The  statements 


C 

Ci  READ  ANO  WRITE  SYSTEM  DESCRIPTION 

C 

WRITE  I b • 1120) 

WRITE  16,11301 
00  SO  1-l.lCEI 

READ  I S, 11601  M, I DNAME I J » , J- l , ICNR I , T8FM, ( HDMT (MPRM|,MPRM« 1,61 ,TMN 
l.TMF  ,ONOFF 

IF  (N.EQ.999)  CO  TO  70 

WRITE  I6.US0)  M,|DNARE(J),J-1,ICN8I,TBFM,  IMOMT  (MPRMl  ,MPRM- 1, 61 
IF  (M.GT.LCEI.OR.M.LT.  1)  GO  TO  60 
IF  (M.GT.INCEt!  INCEI'I 
IF  (ONOFF.EO.OUTI  Mill  1-0 
DO  60  MPRM-1,6 
THOM ( | ,MPRM| ■HDMTt  MPRN I 
60  CONTINUE 

MTBFI11-TBFM 
TINENI 1 1 «TMN 
TIMEFI 11-TMF 
SO  CONTINUE 

READ  IS, 11601  N 
IF  IM.EQ.999)  GO  TO  70 
60  WRITE  16,11601 
STOP 

TO  CONTINUE 


are  useil  lor  input/output  of  the  system  description.  The  first  two  lines  write  the  headings  in 
preparation  for  the  output. 


I he  next  statement  starts  a loop  which  allows  for  the  maximum  input  of  2.)  subsystem 
definitions,  lach  subsystem  definition  consists  of  the  equipment  names.  MTBF,  nuclear 
damage  thresholds  anil  a logical  flag  (ON  OFF)  which  can  housed  to  suppress  a subsystem 
(prevent  its  being  considered  farther).  If  the  system  number  that  is  input  is  equal  to  uoq, 
the  program  switches  control  out  of  the  loop  (i.e.,  "70  CONTINUF")  signifying  there  are  no 
more  inputs  of  this  type.  However,  if  a system  configuration  is  input,  the  program  proceeds 
to  write  out  information.  The  next  statement  is  used  to  check  for  blanks  (subsystems  left 
out)  or  too  many  subsystems  (relative  to  LCEI)  in  the  input  data.  If  either  of  these 
conditions  exist,  program  control  is  transferred  to  an  error  message  which  is  written  out  and 
the  program  stops.  The  next  statement  updates  the  counter  which  keeps  track  of  the  total 
number  of  subsystems  input. 


The  next  statement  is  used  to  set  a tlag  if  a subsystem  is  being  suppressed.  The  internal 
do-loop  which  follows  stores  the  damage  threshold  information  with  the  subsystem  index. 
The  next  three  statements  take  the  MTBF,  the  time  on  and  time  off  information  and  stores 
them  by  the  subsystem  index. 


c 


JTCG/AS-76-S-004 


l 

< 


Program  control  is  then  transferred  to  the  beginning  of  the  system  configuration  loop 
to  read  the  next  subsystem  description.  Should  the  program  go  through  the  loop  23  times 
and  never  transfer  control  to  the  statement  which  says  “70  CONTINUE”,  the  program  will 
read  another  card.  If  this  card  contains  another  system  description,  the  program  will  write 
an  error  message  and  stop.  However,  if  the  input  contains  a “999”  for  the  subsystem 
number,  control  will  be  transferred  to  the  “70  CONTINUE”  statement  and  program 
execution  will  flow  as  described  in  the  next  set  of  statements. 

The  statements 
C 

C3  INITIALIZE 

C 

ET-0.0 
11*0 
T 1 *0.0 


initialize  the  values  of  expected  number  of  targets  killed,  the  previous  event  type,  and  event 
time. 

The  statements 


REMIND  6 
REWIND  3 


rewind  the  tapes  on  units  three  and  four  so  that  they  will  start  at  the  beginning. 

The  statements 
C 

C4  INITIALIZE 

C 

80  CONTINUE 

00  85  M*I,INCEI 
RN(MI*0.0 
85  CONTINUE 

begin  the  event  description  loop  and  zero  out  the  lethal  radius  storage  locations  for  all 
subsystems. 

The  statements 


C 

C5  read  and  write  event  description 

c 

READ  (5*  1170*ENO*UOO>  IEVENT  t T2  .NC . I WEAPN(I)  , 1-1  > 71 , MOT 
WRITE  (6*12301 
00  90  N-1,3 
CN(N)-OFFtN) 
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IF  INC.NE.OI  CN(NI«DEF(NI 
90  CONTINUE 

WRITE  (6,11801  IE VENT,  (CN(N),N*1«3),T2« (WEAPNI  !)»I  = l,7),MDT 


are  used  to  read  and  write  the  event  description  (one  of  many):  event  number,  time,  type 
(defensive  or  offensive),  weapon  description,  and  a mission  descriptor  defining  the  environ- 
mental conditions  at  the  time  of  the  event. 


The  statements 


00  100  L*l,INCEt 
l Z ( L I *L 
100  CONTINUE 


are  used  to  copy  the  subsystems  ordinal  number  into  the  vulnerability  index  array  as  initial 
values  which  are  changed  later,  only  if  there  is  a nuclear  event. 

The  statement 

IF  (NC.EQ.O.OR.NC.EQ.S  .OR.  NC.EQ.7I  GO  TO  200 

transfers  control  to  a later  section  of  the  program  if  the  event  is  not  an  electronics  mode 
nuclear-definsive  event. 

The  next  set  of  statements 


C 

C 6 COMPUTE  LETHAL  RAOI  I AND  ASSIGN  INOEX 

C 

DO  120  I * 1 , NC 

READ  (5,12601  NPOINT 

N2-2*NPOINT 

READ  (5,1220)  I SUB (I) , ( CURVE( I , J I , J-I.N2 ) 

ts-isuem 

DO  120  M-l.INCEI 
00  110  IJ«1,N2 
GRAPH( IJI«CURVE( 1,1 J) 

110  CONTINUE 

CALL  TWODNZ  (THOM! M, IS  I, 1, GRAPH, NPOINT, 2, RESUL » 
RA(M)«RESUL(2I 
RN(MI«AMAXl (RAIMI ,RM(MI I 
120  CONTINUE 

INE-INCEIM 

RM(INE»«0.0 

DO  160  L00*1 , I NCE I 

INIT-LOOM 

00  130  J«IN!T,!NE 

IF  (RM(J).LE.RM( LOO) I GO  TO  130 

STORE-RMI J) 
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RM( JI*RM( LOOI 
RN( LOO) “STORE 
I TORE ■ l Z t LOO) 

IZ(LOO)-IZ(  J) 

IZI Jl-ITORE 
13C  CONTINUE 
140  CONTINUE 

WRITE  (6,1250)  (RM(M) ,M“1, INE) 


will  bo  discussed  in  two  parts  in  order  to  clarify  their  use.  These  statements  are  used  for 
nuclear  defensive  events  only. 

The  first  set  encompasses  statements  from  "DO  1 20”  through  “1  20  CONTINUE”.  The 
statements  through  “1 10  CONTINUE”  are  used  to  read  nuclear  weapon  effect  data  and 
store  the  information  in  GRAPH  for  four  damage  mechanisms.  The  CALL  TWODNZ  applies 
a linear  interpolation  routine  to  derive  the  lethal  radius  associated  with  each  subsystem  and 
each  of  the  four  damage  mechanisms. 

The  two  statements  following  the  CALL  statement  detect  the  maximum  value  of  the 
lethal  radius  of  the  subsystem,  over  all  four  damage  mechanisms. 

Pie  second  set  of  statements  beginning  with  INE  = INCF1  + 1 are  used  to  reassign  the 
lethal  radii  in  descending  order  and  to  assign  a vulnerability  index  number  ( 1 2.)  to  each 
subsystem,  where  I corresponds  to  the  maximum  lethal  radius.  The  lethal  radii  values  are 
subsequently  printed  out. 

The  next  statements 


C 

C7  COMPUTE  SUBSYSTEM  RELIABILITY  IN  TRANSITION 

C 

200  CONTINUE 

00  230  M “ l » INCEI 
PHIM)«1.0 

IF  (NKMI.EQ.OI  GO  TO  230 
IF  I T2.LE.TINENIM) I GO  TO  230 
IF  t T2.GE.TIMEFIM) ) GO  TO  210 
DELTAT-T2-T1 

IF  I Tl.LE.TIMENIM))  OE L TAT-T2-T I MENIM ) 

GO  TO  220 
210  0ELTAT»0.0 

IF  ( TIM£F(M).GT.T1I  OE LTAT-T IMEF (M) -T1 
220  PMt  M ) “ 1, O-OELT AT/MTBFI M ) ♦OE LTAT**2/ I 2»0*MTBF ( M I **2 I 
230  CONTINUE 
T1«T2 

00  240  L«l, INCEI 
JVM  VI  L) 

PMMI LI-PMt JV) 

240  CONTINUE 
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arc  used  to  compute  the  reliability  of  every  subsystem  which  is  time  dependent.  The 
reliability  equation  is: 


dk  = exp  (-  At  + MTBF)  = I - At/MTBF  + At^MTBF^ 


where 


dk  represents  subsystem  reliability  in  the  transition 
At  is  the  elapsed  subsystem  time  in  the  transition 
MTBF  is  the  mean  time  between  failures 


Tire  statements 


C8  initialize  state  ano  state  probabilities 

c 

NZ-0 
NZ  1*0 

00  250  I-I.LCEI 
KHI I 1-0 
250  CONTINUE 

00  260  J-I.L72 
PCAPI JI-0.0 
260  CONTINUE 
JCAP-0 

I F I IEVENT.NE.il  GO  T3  290 
I COUNT -1 

00  27 0 H-l.INCEI 

1 T I H I ■ l 

I TTt  HI *1 

IF  IHIIHI.EO.il  GO  TO  270 
t T(H 1*0 
ITT!  HI-0 
270  CONTINUE 


initialize  the  current  state  binary  vector  and  probability  of  all  mode  sequences  at  zero.  If  the 
event  is  not  the  first,  program  control  is  transferred  to  a later  part  of  the  program.  However, 
if  it  is  the  first  event,  the  prior  state  vectors  IT  and  ITT  are  initialized  to  all  “ones"  (except 
deleted  subsystems)  representing  a perfect  state  at  takeoff  and  no  other  prior  states  are 
considered  for  this  event. 


The  statements 


C 

C9  ABORT  FLOW  CONTROL 

C 

IFINAB0RT.NE.0IG0  TO  275 
GO  TO  278 
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differentiate  between  abort  events  (used  for  simulating  the  first  event  of  an  abort  mission) 
and  non-abort  events  and  directs  the  program  control  accordingly. 


For  those  first  abort  mission  events,  the  statements 


10  INITIALIZE  ABORT  STATES 

275  IF  (I  EVENT. NE.l)  GO  TO  290 

REAOmJEVENT,  JCOUNT, t PJ(U  ,L*1, JCOUNT) 
IFIJEVENT.NE .N ABORT  I GO  TO  275 
DO  2 77  I * 1 • JCOUNT 

phi  >*pjii> 

277  PJII l-P.O 
GO  TO  305 


read  the  state  probability  on  tape  unit  3 into  prior  state  probability  PI  for  the  event  number 
indicated  by  (equal  to)  NABORT.  In  addition,  the  current  state  probability  PJ  is  initialized 
to  zero.  Program  control  is  transferred  to  a later  section  of  the  program  which  generates  the 
initial  state  of  the  system. 


The  statements 


C 

C!1  INITIALIZE  NORMAL  STATES 

C 

278  00  280  L *1 , 2 56 
PIILI-1.0 
PJID-0.0 
280  CONTINUE 
GO  TO  310 

290  00  300  1*1, JCOUNT 
PI(I)<PJ(II 
PJt 11*0.0 
3C3  CONTINUE 
305  CONTINUE 
ICOUNT  *0 


initialize  the  non-abort  state  probabilities  for  the  prior  event  (PI)  equal  to  one  and  for  the 
current  event  (PJ)  equal  to  zero.  Since  the  event  is  not  an  abort  event,  the  prior  event 
probabilities  are  set  equal  to  the  current  event  probabilities  from  the  last  event,  and  the 
current  event  probabilities  are  reinitialized  to  zero.  The  prior  state  counter  is  also  set  equal 
to  zero. 
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The  statements 


C 

CI2 

C 

310 

315 

320 

330 

360 

37C 

380 

390 

400 

405 

406 

407 


GENERATE  A PRIOR  STATE  OF  THE  SYSTEM  UT» 

IFINA8ORT.EQ.0IG0  TO  315 

GO  TO  320 

CONTINUE 

IF< IEVENT.EQ.IIGO  TO  440 

CONTINUE 

00  330  I-WINCEI 

ITU  »-l 

CONTINUE 

IF  (NZ.LT.l)  GO  TO  413 
1 1 1-1 
IT( 1 111*0 

IF  INZ.LT.2)  GO  TO  413 
112*111*1 
ITU  121*0 

IF  ( NZ .IT. 3)  GO  TO  410 
1 1 3*  1 1 2*  1 
ITU  131*0 

IF  (NZ.LT.4)  GO  TO  413 
II 4* 1 1 3*  1 
ITU  141*0 

IF  (NZ.IT.5I  GO  TO  410 

115*114*1 

ITUI5»*0 

IFINZ.LT. 6IG0  TO  410 
1 16*115*1 
ITU  161*0 

IFINZ.LT. 7IG0  TO  410 
1 1 7*  II  6*  1 
IT(!I7I*0 

IF (NZ.LT .81  GO  TO  410 
1 1 8* 1 1 7* 1 
ITU  181*0 


are  used  to  place  exactly  NZ  zeros  into  the  prior  state  vector.  This  block  of  code  works  with 
the  block  of  code  just  prior  to  Compute  State  Probabilities,  to  generate  all  prior  states. 

The  statements 


C 

03  TEST  FOR  OELETEO  OR  CRITICAL  SUBSYSTEMS 

C 

410  CONTINUE 

DO  420  MM*l tlNCEI 
IF  (MUMMJ.EQ.il  GO  TO  420 
IF  ( ITIMMI.EQ. II  GO  TD  840 
420  CONTINUE 

IF(NC.EQ.T)  GO  TO  425 
IF (MCR .E 0*0 1 GO  TO  422 
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IF(| T(MCRI.EQ.O)  GO  TO  840 
422  IF (MPR .EQ.O ICO  TO  424 

IFIIT(MPRI.EQ.O)  GO  TO  840 

424  IFtMAV.EO.OIGO  TO  424 
|F(l T(HAV).EQ.O)  GO  TO  840 

425  I COUNT > {COUNT* I 


ore  used  to  test  for  deleted  and  critical  subsystems.  This  block  of  code  is  identical  to  the  one 
in  Program  I. 


The  statements 


4 RESHUFFLE  SUBSYSTEM  ORDER 

00  4 JO  L-l.INCEl 
JVlVtLI 
ITTm-IT'JVl 
430  CONTINUE 
440  JC OUNT *0 


reorder  the  subsystems  in  preparation  for  use  in  the  transition  algorithm. 
The  statements 


GENERATF  A CURRENT  STATE  OF  THE  SYSTEM  (KM I 


C 

CIS  GENERATF  A CURRENT  S 

C 

450  CONTINUE 

00  460  I«1,INCEI 
KMIII-1 
460  CONTINUE 

IF  (N2I.LT.I)  GO  TO  540 
LLl-l 

490  KMlLLI 1*0 

IF  IN2I.LT.2I  GO  TO  540 
LL2-LLU1 
500  KMILL2I-0 

IF  (N2I.LT. 3)  GO  TO  540 
LL3-LL2M 
510  KMILL3 1*0 

IF  (NZI.LT. 4)  GO  TO  540 
LL4-LL3M 
520  KM(LL4I*0 

IF  I N2I.LT. 51  GO  TO  540 
LL5-LL441 
530  KWILL5I-0 

I F(N2 I .LT.6IG0  TO  540 
LL6-LL5*! 
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5 35 

KNILL6 1*0 

IFtN7I.LT. 

TIGO 

TO 

550 

LL  7«LL6»1 

536 

KWILL7 1*0 

IFtN2I.LT. 

81  GO 

TO 

550 

LL8-LL7*1 

537 

KMILL8I-0 

550 

CONTINUE 

are  completely  analogous  to  the  initial  state  generation  and  are  identically  coded. 


The  statements 


CIS  TEST  FOR  SUPPRESSED  SUBSYSTEMS 

C 

DO  550  MM-l.INCEI 
IF  I MIIMMI.EQ.il  GO  TO  550 
IF  I KMtMMI.EQ. 1)  GO  TO  710 
550  CONTINUE 

IFINC.EQ.7IGO  TO  555 
IFtMCR.EQ.OIGO  TO  552 
IFtKMtMCRI.EQ.OIGO  TO  710 
552  IFtMPR.EQ.0IG0  TO  555 

IF |K  Ml MPR I .EQ. 0) GO  TO  710 
555  lFtM5V.EQ.OIGO  TO  555 

IFIKM(MAVI.EQ.O) GO  TO  710 
555  JCOUNT»JCOUNTM 


are  completely  analogous  to  the  prior  state  test  and  are  identically  coded. 
The  statements 


C 

Cl?  RESHUFFLE  SUBSYSTEM  ORDER 

C 

00  560  L-l.INCEI 
JV-IVtU 
KMMt  L I *KMI JVI 
560  CONTINUE 


are  completely  analogous  to  the  prior  state  subsystem  reordering  and  are  identically  coded. 

The  statement 
C 

C1^  COMPUTE  THE  SYSTEM  STATE  TRANSITION  PROBABILITY 

C 

IF  III .NE .61  GO  TO  1500 
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allows  the  first  transition  algorithm  to  be  used,  which  applies  to  electronics  mode  quick 
conventional  threat  type  events.  The  test  parameter  is  II.  the  prior  event  type. 

The  statements 


Cl°  QUICK  CONVENTIONAL  THREAT  DAMAGE  AND  RELIABIIITV 
C 

SUNTRK-0.0 

DO  620  K*l,KMAX 

SUMTRA-0.0 

00  610  L-l.LMAX 

TRANS*  1*0 

00  600  N-l, INCEI 

PCSURVIM.L.K)  - 1.0  - PCKILLtM.L.Kl 
IP  ( ITTtMI-KWMtM) I 570*580.590 
510  TRANS-0.0 
GO  TO  620 

580  IP  I I T T( M | . EQ. 0. 01  GO  TO  600 

TRANS*PCSURVtM,L,KI*PMMtM|*TRANS 
GO  TO  600 

59C  TRANS-t 1.0-PCSURVt M,L. KI*PMM|MI l*TRANS 
tCC  CONTINUE 

SUNTRA*TRANS»SUMTRA 
6 l C CONTINUE 

SUMTRK  *t SUMTRA/LMAX)*PMtSSIK, I COUNT  I ♦SUMTRK 
620  CONTINUE 

TRANS*SUMTRK 
GO  TO  700 

15CC  IF( I 1.NE.7)  GO  TO  630 


compute  the  transition  probability  TRANS  for  the  state-pair  under  consideration.  The  last 
statement  transfers  control  to  the  nuclear  transition  algorithm  if  appropriate.  Otherwise,  the 
prior  event  is  deemed  to  be  a vehicle  mode  application,  which  requires  the  slow  threat 
damage  transition  algorithm  which  follows. 

The  statements 


C 

C?3  SLOW  THREAT  DAMAGE 

C 

TRANS  « 1.0 
NTFN*0 
I TFN-0 

IINCEt  MNCEI/2 
DO  1510  l*l»IINCEt 

IFtKMMtLI.EQ.il  NT  FN*  2*  MIINCEI/21 -Li *NTFN 
1510  IF  t ITTtLI.EQ.il  ITFN*2**lt  tNCEl /21-LIMTFN 
NTDN*0 
I T0N*0 

J INCEI  *1  INCEt*I 

00  1520  L-JINCEI .INCEI 

IF (KMMtL I.EQ.l INT0N*2**t INCEI-Ll »NT0N 


127 





4 


JTCG/AS-76-S-004 


152  0 IF  ( I TTILI.EQ.lt  ITON*2**(  INCE  l-L  t MTON 
NT  TF  *N  TF  N*1 
I T TF  * I TFN*l 
NTTD=NTON«-l 
I T TO* I TDN*1 

IF(  ITTD.GT.  ITTF.OR.NTTD.GT.NTTF)GO  TO  1605 
IF(  ITTF.LE.l.OR.ITTO.LE.l)  GO  TO  1605 
!F(NTTF.GT.IITTF-l).0*.NTTD.GT.UTT0-inGO  TO  1605 
GO  TO  1610 
16C5  TRAN  SI  *0 

GO  TO  1720 
1610  TRANS l *0 

IF(NTTF-(ITTF-III16 70, 1620 ,1620 
1620  IFINTTD-I I TTD-ll 11630, 1625, 1625 

1625  DO  1627  I*NTTF,MLTH 
SUMPTD*0 

DO  1626  J*NTTD » I 

1626  SUMPTD*PTON( I « Jl ♦SUMPTD 

1627  TR ANS1*PTFN( II *SUMPTD*  TRANS1 
GO  TO  1720 

163C  00  1631  I*NTTF,MLTH 
1631  TRANSl=PTFN(  1 1 *PTDN( I , NTTOI *TR ANSI 
GO  TO  1720 

167C  I F INTTO-I I TTD- lit  1680*  1675,1675 

1675  DO  1676  J*NTTD,NTTF 

1676  TR  ANSI  =P TON ( NTTF  « J ) *TR  ANSI 
TRANS1=PTFN(NTTFI*TRANSI 
GO  TO  1720 

168C  TR  AN  SI =PTFN( NTTF I *PTDN ( NTTF , NT TO I 
1 720  TRANS=TRANS1*TRANS 
GO  TO  700 


compute  the  transition  probability  for  the  vehicle  mode.  The  statements  through  1520 
convert  the  binary  t and  r values  to  decimal  values.  The  next  four  statements  change  the 
time  reference  from  the  prior  event  to  the  current  event.  The  next  six  statements  eliminate 
unwanted  states  by  setting  the  transition  probability  to  zero.  The  remaining  statements  are 
as  described  in  the  flow  chart  (Figure  8).  The  four  transition  probability  equations  are  given 
in  Table  6. 

The  statements 


C 

C’l  QUICK  NUCLEAR  DAMAGE  AND  RELIABILITY 
C (ALSO  NON-DEFENSIVE  TRANSITIONS) 

C 
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630  TRANS^l.O 
NF»1 

I Fill. Eg. 0)  GO  TO  640 
TR AN  S*QG( ICOUNTI 
NF»0 

640  CONTINUE 

00  690  M*  1 , INC  El 
IF  ( ITTiMI-KMMiMtl  653.660.670 
650  TRANS-0.0 
GO  TO  700 

660  IF  ( ITTiMI.EQ.OI  GO  T1  680 
NF*1 

TRANS»TRANS*PMM<M> 

GO  TO  690 

67C  TR  ANS«TRANS*(  l .O-PNHiN)  | 

68C  IF  (NF.EQ.01  TRANS«TRANS*QPRM(M, ICOUNTI 
69C  CONTINUE 


compute  the  transition  probability  for  the  electronics  nuclear  events  and  non-defensive 
electronics  events.  The  flag  NF  = I in  the  third  statement  indicates  that  nuclear  damage  is 
not  considered.  This  is  revised  immediately  if  the  prior  event  was  electronics-nuclear,  but  is 
reset  to  1 when  the  first  surviving  subsystem  is  encountered.  For  non-nuclear  events,  only 
reliability  effects  occur. 

The  statements 
C 

C 22  COMPUTE  THE  STATg  PROBABILITIES 

C 

70C  PJ(  JCOUNT)-TRANS*PI(  ICOUNTI«-PJt  JCOUNT) 


accumulate  (add)  the  probabilities  for  all  transitions  into  the  same  current  state.  This  sum  is 
the  current  state  probability.  This  is  repeated  for  all  current  states  as  the  state  do-loop  is 
executed.  The  equation  is: 

K 

Pj.  n = s (Pj/j>n)x(Pi.n) 

The  statements 


C 

C23  RESET  CURRENT  STATF  IKN)  ZERO  LOCATIONS 

C 

710  IF  (N2I.LT. 1.0R.N2I.GT.B)  GO  TO  820 

GO  TO  (800.780.760*  740 .720. 715*7 14, 7131* N2 I 
713  KW(LL8)*l 
LL8«LL8*1 

IFtt  LLB-8).LE.( INCE 1-N2 1 1 1GO  TO  537 
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714  KW(LL7)=l 
LL  7*LL7*1 
IFCILL7-7I. 

715  KW( L L6 1*1 
LL6*LL6*l 

I F < ( LL6-6  > . 

72C  KW(LL5>»1 
LL5*LL5*l 
I F ( ( LL5-5 ) < 

74C  KW(LL4»=>1 
LL4=LL4*1 
£FU  LL4-4I, 

760  KW<LL3)=1 
LL  3*  LL3* 1 
I F ( ( LL  3-3  ) . 

760  K W ( L L2 ) = 1 
LL  2-LL  2*- 1 
IF(( LL2-2I. 

8C0  KW( L L 1 1 3 1 
LL  l=LL  1*- 1 
I F C < LLl-U. 

820  CONTINUE 
NZ  l*NZIH 
IF  (NZl.LE. 
NZ  1*0 


LE.IINCEl-NZI) ) GO  TO  536 


LE.  (INCEI-NZmCO  TO  535 


LE. (INCEI-NZIi)GO  TO  530 


LE.  (INCEI-NZniGO  TO  520 


LE. (INCEI-NZIi)GO  TO  510 


LE. ( INCEI-NZIIIGO  TO  500 


LE.  (INCEI-NZNIGO  TO  490 


NZT ) GO  TO  450 


work  in  conjunction  with  an  earlier  block  of  code  called  “Generate  Prior  State”  to  do  just 
that.  The  role  of  these  statements  is  to  increment  the  location  of  the  right-most  zero  in  the 
state  vector  until  it  has  gone  as  far  as  it  can  to  the  right  and  then  increment  the  next  zero 
location,  etc.  If  all  zeros  have  been  moved  as  far  as  possible  to  the  right,  another  zero  is 
brought  into  play.  However,  if  the  added  zero  exceeds  the  number  provided  by  the  user 
(NZT  test),  the  process  of  generating  states  is  arrested  and  the  program  continues. 

The  statements 


C 24  ABORT  FLOW  CONTROL 

C 

IF (NABORT.EQ.O I GO  TO  835 
GO  TO  840 
835  CONTINUE 

IF  I IEVENT.EQ. II  GO  TO  970 

control  state  generation  for  first  events.  If  the  event  is  not  an  abort  and  is  the  first  event, 
there  is  only  one  initial  state,  and  the  prior  state  generation  is  bypassed.  If  the  event  is  an 
abort,  the  normal  state  generation  process  is  continued. 
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The  statements 


c 

C?5  RESET  PRIOR  STATE  f I T • ZERO  LOCATIONS 

C 

840  IF  I NZ.LT.l .OR.NZ.GT.8 » GO  TO  930 

GO  TO  ( 9 30 1 9 10  «890 * 870  » 850 • 845 , 844 , 843) • NZ 
843  ITII181*1 
1 1 8*  1 1 8*  l 

IF((  1 18-8I.LE.  UNCEI-NZ  I IGO  TO  407 
344  IT( I 171*1 
II  7-  II  7*  1 

IFIII I7-7I.LE.I INCEI-NZIIGO  TO  406 
845  IT( 1 16  1*1 
II 6* II 6*1 

IF((  I 16-6). LE.  (INCEI-NZ)  IGO  TO  405 
850  I T < I 151*1 
115*115*1 

IFII 115-51. LE. (INCEI-NZ1IG0  TO  400 
870  ITII 141*1 
114*114*1 

IFII I I4-4I.LE. ( INCEI-NZ11  GO  TO  390 
890  I TU  131*1 
113*113*1 

IFII II3-31.LE.I INCEI-NZ II  GO  TO  380 
910  IT! 1 121*1 
112*112*1 

IFII 1 12-2 1 . LE. I INCEI-NZI)  GO  TO  370 
530  I T! I 1 1 1* l 
II  1*11 1*1 

IFIII  U-ll.LE.  UNCEI-NZMGO  TO  360 
550  CONTINUE 
NZ*NZ* l 

IF  I NZ .LE.NZT1  GO  TO  320 
NZ  *0 

570  CONTINUE 


arc  analogous  to  the  “Reset  Current  State  Zero  Locations”  block  of  code  and  is  identically 
coded. 


The  statements 


C 

C 26  WRITE  OUTPUT  STATE  PROBABILITY  TAPE 

C 

IFINABORT.NE.OIGO  TO  979 
JE  VENT *1  E VENT 

WRITE! 3) JEVENT , JCOUNT, < PJ(L I ,L-l , JCOUNT I 

prepare  a record  of  state  probabilities  on  a normal  mission,  which  may  be  employed  to 
initialize  an  abort  mission.  The  variable  JEVENT  is  used  to  denote  a normal  mission  event 
number  so  that  the  IEVENT  variable  may  be  used  on  an  abort  mission  without  ambiguity. 
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The  statements 


C 

C?7  COMPUTE  THE  MOOE  SEQUENCE  PROBABILITIES 

C 

<.79  READ  t A t 1280 ) ( l JN(  t JK I , I JK»l , JCOUNT  I 
JCAP*0 

DO  980  I JK* l » JCOUNT 
J- I JNt I JKt 

PC  API J) -PC API J)*PJI UK) 

JCAP-MAXOI JCAP.J) 

S8C  CONTINUE 

PARI  VE  »PCAP(1)»PCAP12I  *PCAP (3) 


are  used  to  read  in  the  mode  sequence  number  array  from  program  one  and  then  distribute 
all  state  probability  into  variables  PCAP(J)  (representing  probability  of  the  Jth  mode 
sequence).  The  distribution  is  made  in  accordance  with  the  mode  sequence  number  (J) 
array,  whose  subscript  is  the  state  number.  PAR1VE  is  the  sum  of  probabilities  of  the  first 
three  mode  sequences.  In  a vehicle  mission,  this  sum  is  printed  out  (later).  The  mode 
sequence  probability  equation  is: 


jmax 


The  statements 


C 

C28  READ  AND  WRITE  OUTPUT  MODE  PROBABILITIES  ANO  CAPABILI TES 

C 

DO  1000  J-l.JCAP 

READ  15,12601  CE  PI  J ) ,F  A ( J)  , PK(  J ) 

IF  INC.EQ.O)  GO  TO  990 

IF(NC.NE.T)  WRITE  16,12701  J, PC  API J I , CEP  I J I , FA  I J I 
IFINC.EQ.7t  WRITE (6, 1275)  J.PCAP(J) 

GO  TO  1000 

990  WRITE  (6,11901  J.PCAPI J I ,PK( J) 

1C00  CONTINUE 

IFINC.NE.7t  GO  TO  9030 
WR I TE I 6, 9020) PAR  I VE 
9C20  FORHATIIX.'PARIVE  » ' * E 12. 5 ) 


are  used  to  write  the  mode  sequence  probabilities,  and  read  and  write  offensive  and 
defensive  capabilities  as  appropriate  to  the  event  type  and  MISDEM  mode.  If  the  current 
event  is  not  a vehicle  simulation,  an  additional  set  of  statements  are  executed.  For  offensive 
events,  the  output  effectiveness  equations  are: 

JCAP 

pK.m=  jJj  pk,J  x pJ,n 

N .1 

Rt(N)=  Z PK,n 
n = 1 
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The  statements 


C 

C?9  tF  EVENT  IS  OFFENSIVE .COMPUTE  ET  - MISSION  EFFECTIVENESS 

C 

903C  II-O 

IF  INC.NE.O)  GO  TO  1030 
00  l 010  J-l » I NCE I 
IVI J»«J 
1)10  CONTINUE 
SUM-0.0 

00  1020  J-l.JCAP 
SUN- SUM^PCAPI  J ) *PK(  Jl 
1C2C  CONTINUE 

WRITE  (6.1200)  (EVENT, SUM 
ET-EUSUN 
WRITE  (6.1210)  ET 
GO  TO  80 


are  used  to  set  the  prior  event  type  variable  11  to  zero  in  preparation  of  the  next  event,  and 
remains  so  unless  the  current  event  is  not.  If  offensive,  the  vulnerability  indices  are 
initialized  equal  to  subsystem  ordinal  number,  and  cumulative  expected  number  of  targets 
killed  is  updated  and  printed  out.  If  the  event  is  not  offensive,  control  is  transferred  to  the 
next  block  of  code. 

The  statements 

C 

c 

CO  IF  EVENT  IS  DEFENSIVE  DETERMINE  SYSTEM  SURVIVABILITY  PARAMETERS 

C 

C 

1030  IF ( NC. EQ. 7)  GO  TO  2500 
II-l 

00  1040  10-1,1 NCE I 
IV(IQt-miQ) 

1040  CONTINUE 

IF  (NC.NE.6)  GO  TO  1070 
1 1 -6 


are  used  to  set  prior  event  type  equal  to  1 and  vulnerability  indices  equal  to  the  rank  order 
dictated  by  this  current  nuclear/conventional  defensive-electronics  event. 


The  statements 

C 

C 3 1 QUICK  CONVENTIONAL  KILLS 
C 

RE AO  (5.1290)  LMAX, KMA X.(R(K),K«l, KM AX ) 
WRITE  (6.1300)  (R(K).K-l.KMAX) 

00  1050  M-l.INCEl 
00  1050  L-l.LMAX 

READ  (5,1310)  (PCKILLIM.L.KI.K-l.KMAX) 
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1050  WRITE  (6,1320)  M,L  , (PC  KILL  ( H,  L ,K) , Ml  ,KM  AX  ) 
KMAXi-KMAX-1 
00  1060  I JK*i • JCOUNT 
PM  I SSI KMAX, l JK)*l.O 
4* I JN( IJK) 

SIGMA*CEP( J)/l.l78 
OE  NUM*2.0*SIGMA**2 
0EN1 75*1 75. 0*DENUM 
EXl-1.0 

00  1060  K*l,KMAXl 

EX2-0.0 

RM2*R( Kl P*2 

IF  (RM2.LT. OEN175I  EX2*EXP(-RM2/0€NUNI 
PMISSI K,I JK) -EXI-EX2 

PMISS(KMAX,IJKI*PMISS( KMAX , I JK ) -PM I SSI K, UK) 
EXl*EX2 
1 C 60  CONTINUE 
GO  TO  80 


are  used  first  to  load  all  subsystem  kill  probabilities  associated  with  KMAX  miss  distances 
and  LMAX  offset  angles,  for  use  in  computing  the  conventional  weapon  transition 
probability  in  the  next  event.  The  statements  after  1050  are  used  to  compute  the 
probability  of  a threat  warhead  miss  in  a circular  zone  whose  inner  radius  is  R(K).  The  three 
equations  used  in  this  block  of  code  are  given  below.  The  standard  deviation  (sigma)  is 
derived  from  the  threat  weapon  CEP  as  follows: 

o = CEP/1.178 

The  result  of  the  integration  of  the  bivariate  normal  distribution  over  limits  and 
Rk+1  > Rk> is: 


pmiss<K)  = exP  <Rk2  I2  <72)-  exp(R(*  + ,/2  a2) 


The  transition  probability  associated  with  such  an  event  is  then  computed  as: 


A'* 

Pj/j  = Pj/i  (given  burst  point  k)  x P (burst  point  k) 

k = 1 


The  statements 


C 

C '2  SL 
C 

2500  11*7 


SLOW  KILLS 


DO  2122  I-l.MLTH 

RE  A0  I 5, 2 120)  FTFNNII),  (FTONNKtJ),  J*1,MLTH) 
2120  FORMAT! 17F3.2) 

IFF*  I- 1 

IFII.E0.1)  GO  TO  2125 
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PTFNtll  - IFTFNNI  II-FTFNNI  IFF)  I 
GO  TO  2129 

2125  PTFNt l l-FTFNNt 1 > 

2129  CONTINUE 

00  2130  J* 1 * NLTH 
JFF-J-l 

IFIJ.EQ.1IG0  TO  2140 
P TONI  I « J I -F  TONNI I, Jl-FTDNNI I.JFFI 
GO  TO  2130 

2140  PTONU.JI  * FTONNI  I * J I 
213G  CONTINUE 
2122  CONTINUE 
GO  TO  80 


are  used  first  to  read  in  the  flight  time  and  abort  deteetion  time  distribution  functions.  The 
statements  after  21  20  are  used  to  compute  the  discrete  probabilities  of  flight  time  and  abort 
detection  time  for  all  values  up  to  and  including  the  duration  of  the  mission  MLTH.  The 
flight  time  algorithm  ends  at  2129;  the  abort  detection  time  algorithm  ends  at  2130. 


The  statements 


C 

C3 3 QUICK  NUCLEAR  KILLS 

C 

1C7C  CONTINUE 

00  1090  I JK-l » JCOUNT 
Q0( I JK  1*0.0 
J* l JN( I JK I 
SIGMA-CEPt Jl/1.178 
DE  NUM-2. 0*S I GMA**2 
DENI  75  = 1 75. 0*0ENUM 
EXl-0.0 
RM2«RM(l)**2 

IF  I RM2.LT. 0EN175)  EX1 * EXPt -RM2/DENUMI 
QOI l JK)-EX1**FA< J) 

DO  1080  M*1 • I NCE I 
QPRMIM.I JK) -0.0 
RMM2-RMC  M ) **2 
RMM12-RM(M*1)**2 

IF  (RMM2.LE.RMM12.0R.RMM12.GE.DEN175)  GO  TO  1080 

E XMl »E XP ( -R MMl 2 /OE NUMI 

EXM-0.0 

IF  (RMH2.LT.OEN175)  EX M-EX P ( -RMM2/ DENUM ) 

QPRM I M ( I JK ) -EX  Ml **F  A I J ) -EXM**F A f J | 

1080  CONTINUE 
1090  CONTINUE 
GO  TO  80 
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are  used  to  compute  the  probability  of  a nuclear  weapon  miss  in  the  circular  zone  whose 
outer  radius  is  RM  (M).  Since  the  miss  is  assumed  normally  distributed,  the  miss  probability 
is  negligibly  small  beyond  a few  standard  deviations.  This  fact  is  used  in  setting  up  a test  for 
RM  to  avoid  underflow  in  the  subtraction  of  two  vanishing  numbers  in  the  fourth  statement 
from  the  bottom.  The  equation  for  miss  probability  is: 

q’k.i  = exP  (-SPtlrE  + 1 1~  °f ) - exP  (-SPdr£  /2  o^) 

Tire  following  statements  are  used  to  rewind  the  tapes  and  stop  the  execution  of  the 
program. 

1 ICO  REWIND  4 
REWIND  3 
STOP 

The  statements 


C 

C34  DEFINE  FORMATS 

C 

1110  FORMAT  I 101 3) 

1120  FORMAT  I 26H1  SYSTEM  CONFIGURATION! 

1130  FORMAT  I 1H0 , 8X .9HEQUIP MENT , 25X , 4HMT8F , 7X, 7HTH0M ( 11 , 6X, THTHDMt 2 » , 7 X 
It 7HTH0MI 3! , 8X, 7MTHDMI 4 I /56 X ,3HG  D* 11X.1HN* 12X. IHB.1SX, 1HT//1XI 
1140  FORMAT  1 1 3 , 1 X, 8 A4, F10. 2 /4E 1 2. 5 . 2X, 2F6.2 , 2X, A2 I 
1150  FORMAT  I IH  ,lX,I3f2X,8A4,F10.2»2X,4(E12.5.2X)) 

1160  FORMAT  I 50H0  TOO  MANY  INPUTS  OR  M IS  OUTSIDE  ALLOWABLE  RANGE! 

1170  FORMAT  4 I3t lX,F6.2flXt I 3 , IX , 7A4, 2X , L l ) 

1130  FORMAT  4 1H0.5X.9HEVENT  NO. , 1 2 , l X , 3H I S , 3A4/ 5X, 1 5HE VENT  OCCURRED  .F 
16. 2f  20H  HOURS  AFTER  TAKEOFF /5X  » 21HEVENT  DESCRIPTION  IS  . 7A4  /5X.L1 
2! 

1190  FORMAT  41HO,lOX,3HJ  I 3 ,2X ,9H PCAP 4 J ! * , El  2. 5, 2X, 7HPK 4 J I -.E12.5I 

120C  FORMAT  4 lHOt lOXf 40HEFFECTIVENESS  FOR  OFFENSIVE  EVENT  NUMBER, 13, IX, 
12HIS.E12.5I 

1210  FORMAT  4 1H0, 10X , 36HCUMULAT I VE  MISSION  EFFECTIVENESS  tS  ,E12.5I 
1220  FORMAT  4 1 3, 2X, 6E12 . 574 4 5X ,6E 12- 5 > ! 

1230  FORMAT  (23H1  EVENT  DESCRIPTION! 

1240  FORMAT  I 13! 

1250  FORMAT  I IH0  //5X.8HRM  ARRAY / 4 1 4X ,4 4 E 12 . 5 ,2X ) » ! 

1260  FORMAT  43E12.5I 

1270  FORMAT  41H0,10X,3HJ  * » I 3 , 2 X , 9H PC AP ( J ) «, E12.5.2X, 8HCEP4 J I «,E12.5, 
12X,7HFA«J»  » ,E 12.51 

127*  FORMAT 4 1H0,10X,3HJ  *, I 3 ,2X ,9HPCAP 4 J ! *,E12.5! 

1280  FORMAT  C 251 3! 

1290  FORMAT  42I3,10F7.0! 

1300  FORMAT  4 IHO. 2X .31HCOMPONENT  PROBABILITIES  OF  KILL/ 

1 17H  COMPi  ELEV#  R», 10IF8.0, 1X1! 

1310  FORMAT  18E10.4! 

1320  FORMAT  4 216,24 4X,E 10.4 ! I 
ENO 


define  the  input  and  output  formats. 
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Subroutine  TWQDNZ 


The  purpose  of  this  routine  is  to  obtain  the  value  of  an  unknown  which  lies  along  a 
given  curve  by  the  linear  interpolation  method. 


Tire  statements 

SUBROUTINE  TWOONZ  < XG, NI Z, Zl ,NX, NZ, ANSI  I 
C 

C*  TWO  OIMENSIONAL  LINEAR  INTERPOLATION  ROUTINE 

C 

DIMENSION  Zitll*  ANSI III 


are  used  to  pass  information  in  and  out  of  the  subroutine  and  allow  storage  for  the  arrays  by 
use  of  the  DIMENSION  statement.  The  arguments  in  the  calling  sequence  are  identified  as 
follows: 


XB  = Given  value  of  independent  variable 

NIZ  = Column  in  which  independent  variable  is  stored 

/.I  = Name  of  array  in  which  the  variables  are  stored.  Although  the  array  is  actually 
a onc-diinensional  array,  it  can  be  considered  to  consist  of  two  columns,  one 
for  the  independent  variable  and  one  column  for  each  dependent  variable. 

NX  = Number  of  values  in  any  column  (same  for  all  dependent  variables) 

NZ  3 Number  of  different  variables  (including  the  independent  variable) 

ANSI  = Answer  vector 

NOTE:  Zl  is  stored  in  the  following  fashion: 

Zl(l).  Zl(2).  . . .ZI(NX).  Z2(l).  . . ,Z2(NXV--ZNZ(1).  . . .ZNZ(NX). 

Zl  is  stored  in  increasing  order  as  a single  vector  with  a one-dimensional 
subscript. 


The  statements 


IXL-(NIZ-lt*NX»l 

C?  IF  GIVEN  VALUE  IS  LESS  THAN  LONER  LIMITS  OF  TABLE.  SET  RESULTS 

C EQUAL  TO  LOWER  LIMIT 

IF  t XG-ZlttXL) ) 70,70,10 
10  LL-IXLM 
LU-I XL ♦NX-l 

C 3 SEARCH  FOR  INTERVAL  IN  WHICH  GIVEN  VALUE  LIES 

00  20  J-LL.LU 
IF  IXG-ZltJI)  30,80,20 
20  CONTINUE 

C<*  IF  GIVEN  VALUE  IS  GREATER  THAN  URRER  LIMIT  OF  TABLE,  SET 

C RESULTS  EQUAL  TO  UPPER  LIMIT 

<!■ 1 XL*NX-l 
GO  TO  BO 
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C 5 GIVEN  VALUE  WITHIN  TABLE*  CALCULATE  INTERPOLATION  FACTOR 

30  RAT-IZK  JI-XGI/IZH  J)-ZltJ-lH 
AO  JP- J-l N! Z-I I *NX 

Cfi  DO  LINEAR  INTERPOLATION  FOR  RESULTS 

00  60  K-l.NZ 

C 7 CHECK  FOR  ZERO  SUBSCRIPT 

IF  UP.EQ.il  GO  TO  SO 

ANS1IKI-Z1I  JP|-RAT*IZ1 1 JPI-Z1I JP-lll 
GO  TO  60 
SO  ANSI  111*21(11 
60  JP*JP*NX 
RETURN 
70  J-IXL 
8C  RAT-0.0 
GO  TO  AO 
END 


perform  the  interpolation  using  a ratio  factor.  Assuming  that  the  values  of  the  independent 
variable  are  called  Vj,  i = 1,.  . .n,  the  'alue  of  Vr  corresponding  to  a given  value,  Tg  is  found 
by  the  following  method: 


where 


I 


Tj<Tg<Tj+1 

Vr  = R [ Vj  + i - Vj  ] + Vj 


There  is  no  extrapolation  beyond  the  values  given  in  the  table.  When  the  given  value  is  less 
than  the  smallest  value  of  the  independent  variable  in  the  table,  this  smallest  value  is  used; 
similarly  when  the  given  value  is  greater  than  the  largest  value  of  the  independent  variable  in 
the  table,  the  largest  value  is  used. 


USER  INFORMATION 
Machine  Requirements 

MISDEM  was  written  for  use  with  the  IBM  370/168  computer.  This  FORTRAN 
program  contains  three  external  references  (EXP,  MAXO,  and  AM  AX  1 ),  all  of  which  are 
basic  ASCII  rbutines.  Since  MISDEM  is  almost  totally  self-contained,  it  could  easily  be 
converted  for  use  on  other  machines  with  some  minor  programming  changes.  Program  2 is 
the  larger,  and  occupies  approximately  168K  bytes  of  storage  to  execute  as  currently 
dimensioned.  Running  time  for  the  test  (verification)  cases  required  less  than  one  minute  for 
both  the  electronics  and  vehicle  modes  combined. 
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Those  specifications  are  for  running  the  M1SDEM  model  at  an  IBM  370/168  facility. 
However,  the  exact  requirements  are  both  machine  and  facility  dependent  and  should  be 
verified  before  running  the  MISDEM  model  at  a specific  facility. 

Conversion  to  CPC  Machines 

Conversion  to  a CDC  machine  is  exceedingly  simple.  In  program  one.  line  4. 

replace  INTEGER*  2 IJN(256) 
with  DIMENSION  1JN(256) 

Between  lines  70  and  71  the  read  statement  must  be  changed  to  read  as  follows: 

REAP  (5.  1010)  I I VI  NT.  T2.  NC.  (WEAPN(l).  1=1.  7).  MOT  II  (101(5))  840.  71  I 
CONTINUE 

The  changes  necessary  in  program  2 are  very  similar,  starting  with  the  first  dimension 
statement: 

replace  INTEGER*  2 IJN(256) 
with  DIMENSION  1JN(256) 

Between  lines  85  and  90  the  read  statement  must  be  changed  to: 

READ  (5.  I 170)  I EVENT.  TC,  NC.  (WEAPNf  1).  1=1.  7).  MDT  II  (EOE(5))  1 100.  87 
CONTINUE 

Uncontrolled  Errors 

Machine  limitations  are  probably  not  significant  in  most  applications.  As  currently 
dimensioned,  six  decimal  places  have  been  allowed  for  the  output  probabilities.  Input  data 
are.  at  present,  probably  resulting  in  inaccurate  output  beyond  the  third  decimal  place.  As 
long  as  the  machine  precision  exceeds  the  output  parameter  precision,  there  is  no  ambiguity 
or  error  in  the  output  caused  by  the  machine.  Uncontrolled  mathematical  errors  occur  only 
in  the: 

1.  Approximation  in  the  reliability  algorithms 

2.  Approximation  in  the  multiple  nuclear  warhead  zone  miss  probability  algorithm 

3.  Assumed  unity  kill  probability  of  a nuclear  warhead  within  a lethal  radius 

4.  Assumption  of  subsystem  failure  independence,  as  discussed  previously 

5.  Assumption  of  failure  independence  between  the  vehicle  and  the  electronics. 
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Controlled  Errors  Versus  Running  Time 

USE  OF  NZT.  In  the  conventional  electronics  mode,  controllable  errors  will  occur  as  a 
result  of  suppressing  states  having  (acceptably)  low  probability.  This  is  accomplished  by 
setting  NZT  <M  where  NZT  is  the  maximum  number  of  allowable  zeros  (failures)  in  the 
state  (binary  -M)  vector,  and  M is  the  number  of  subsystems.  Because  the  states  having  a 
large  number  of  failures  will  normally  be  assigned  low  probability,  their  deletion  will  not 
significantly  affect  the  output. 

This  is  not  true  in  the  nuclear  electronics  case.  A significant  amount  of  probability 
mass  will  be  allocated  to  states  having  several  vulnerable  subsystems  failed.  This  is  due  to  the 
tendency  of  a user  to  make  the  damage  thresholds  common  for  several  subsystems,  due  to 
lack  of  contrary  information.  When  the  nuclear  threat  is  significant,  it  is  because  failures 
occur.  Then  the  probabilities  assigned  to  group  failures  are  high  due  to  the  unity  correlation 
of  failures  in  the  group  having  the  same  lethal  radius.  St)  the  larger  the  size  of  such  groups, 
the  larger  must  be  the  value  of  NZT  to  ensure  accounting  for  significant  probability  mass. 

rhe  use  of  NZT  < M is  a way  to  reduce  machine  use  time  in  electronics  modes; 
therefore,  its  choice  is  important.  The  error  caused  by  NZT  < M is  difficult  to  predict,  but  it 
can  be  (in  some  cases)  observed  in  the  output  by  adding  all  values  of  PCAP(J)  (mode 
sequence  probability)  at  a given  event.  If  the  sum  falls  short  of  1.0.  the  defect  is  due  to 
NZT.  The  cases  where  this  technique  works  is  where  states  have  not  been  suppressed  by  the 
subsystem  criticality  ( MCR,  MPR,  MAV)  tests  (which  would  otherwise  result  in  an 
additional  loss  of  probability  mass,  that  would  thereby  hide  the  loss  caused  by  NZT  alone). 

The  variable  NZT  cannot  be  used  to  reduce  the  running  time  in  the  vehicle  mode 
because  it  would  introduce  large  errors  as  follows.  If  NZT  is  less  than  INCFI  (the  size  of  the 
state  vector),  the  state  having  all  zeros  is  suppressed.  This  is  the  state  that  causes  flight 
failure  in  a single  epoch,  and  could  carry  a significant  probability  mass  if  the  damage 
mechanisms  are  relatively  quick. 

QUANTIZATION  ERRORS.  The  use  of  a relatively  small  number  of  events  to 
represent  a larger  number  of  events  in  the  real  world  results  in  time-quantization  errors.  For 
the  electronics  case,  the  simulated  system  is  not  allowed  to  change  its  response  to 
encounters  (i.e„  employ  a different  mode  sequence)  except  at  the  specified  event  time, 
whereas  in  the  real  world,  the  response  could  have  changed  several  times  in  the  time  interval. 
The  state  probability  distribution  could  be  in  error  at  the  end  of  the  interval  for  defensive 
events,  as  a result  of  the  instability  caused  by  the  feedback  of  countermeasures  effectiveness 
to  survivability  of  countermeasures.  When  precise  results  are  required,  the  quantization 
interval  can  be  decreased  to  suit.  When  the  mission  results  are  within  desired  values  of  the 
apparent  asymptotic  values,  the  interval  need  be  decreased  no  further.  The  running  time 
increase  is  directly  proportional  to  the  number  of  events  when  it  is  large. 

In  the  vehicle  mode,  there  is  a quantization  error  due  to  the  fact  that  the  scenario 
event  times  may  not  coincide  with  the  regular  intervals  required  in  this  mode.  One  cure  for 
this  is  the  use  of  a larger  number  of  events,  although  the  impact  on  running  time  is  much 
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greater  than  in  the  electronics  case,  The  running  time  in  the  vehicle  ease  is  proportional  to 
2 (4  logT  N)  for  large  N.  where  N is  the  number  of  events.  Another  possible  cure  for  this 
problem  is  manipulation  of  the  input  data  to  provide  flight  time  and  abort-detection  time 
distribution  functions  at  regular  intervals,  which  are  then  compatible  with  M1SDEM. 

USE  OF  MCR,  MPR  AND  MAV.  The  use  of  critical  subsystems  (MCR.  MPR  and  MAV ) 
in  the  electronics  case,  to  reduce  running  time,  suppresses  those  states  having  zeros  in  those 
subsystems,  and  results  in  a loss  of  probability  mass  in  those  modes  not  requiring  such 
subsystems.  If  these  latter  modes  are  not  considered  significant  to  the  output,  the  error  is 
acceptable.  The  three  variables  (MCR.  MPR  and  MAV)  are  identical  in  function. 

Application  Notes 

DELETION  OF  SUBS1  STEMS.  In  the  electronics  mode,  subsystems  may  be  deleted 
on  a mission  basis  only,  by  means  of  the  variable  ONOFF  in  the  system  configuration  data. 
This  is  not  applicable  to  the  vehicle  mode,  which  does  not  use  subsystems  in  its  state  vector. 

CONTROL  OF  SUBSYSTEM  ON  AND  OFF  TIMES.  For  purposes  of  computing 
reliability  in  the  electronics  mode,  the  program  requires  a single  turn-on  time  and  a single 
turn-off  time  for  each  subsystem.  The  variables  TMN  and  TMF  (for  time-on  and  time-off. 
respectively)  cannot,  however,  be  used  to  delete  subsystems  (a  function  reserved  for 
ONOFF).  In  fact,  it  is  necessary  that  the  user  select  TMN  and  TMF  in  such  a way  that  the 
subsystems  are  “on”  for  every  event  whose  mode  sequence  logic  requires  them  “up”  for  any 
mode  (unless  the  user  is  willing  to  accept  the  error  resulting  from  the  inconsistent  inputs). 

SOURCES  OF  DATA.  There  is  no  known  documentary  source  of  data  for  the 
abort-detection  time  distribution  function  at  present.  Some  examples  of  these  functions  are 
given  to  aid  the  user  in  deciding  how  to  select  input  values  of  this  function.  Consider  first 
the  M1SDEM  application  to  involuntary  aborts.  Die  latter  aborts  are  defined  to  be 
independent  of  the  time  the  crew  detects  a need  to  abort.  To  model  the  involuntary  aborts, 
the  input  abort-detection  time  distribution  functions  could  be  set  to  zero  for  all  time.  The 
abort  then  occurs  whenever  the  flight  time  runs  out.  not  before.  Alternatively,  it  may  be 
assumed  that  the  crew  has  perfect  knowledge  of  the  amount  of  flight  time  remaining.  In 
order  to  model  this  situation,  the  abort-detection  time  distribution  functions  would  be  set 
to  1.0  for  all  time.  The  abort  then  occurs  in  accordance  with  the  mode  sequence  logic;  e.g.. 
when  the  flight  time  remaining  becomes  less  than  the  normal  mission  duration. 
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List  of  Abbreviations  and  Symbols 
(Simulation  Model  — Program  One) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

AMAXI 

IBM  utility  routine  to  choose  the  maximum  between 
two  numbers 

BLK 

A blank  which  is  used  to  clear  the  subsystem  names  in 
core 

CN(N) 

Ultimately  contains  the  event-type  descriptions  ‘'offen- 
sive" or  “defensive” 

CNAME(J.I) 

Permanent  storage  location  for  the  1th  subsystem  name 

DEF(N) 

A word  containing  the  Holerith  character 
“DEFENSIVE" 

DNAME(J) 

Temporary  storage  location  for  the  subsystem  names 

FNAME 

Name  of  the  subfunction  and/or  mode 

HDMT(MPRM) 

Um'.R) 

Temporary  location  for  the  damage  threshold  for 
designated  MPRM  damage  mechanism  (gamma  “dot", 
neutron,  blast,  and  thermal) 

1 EVENT 

Event  number 

UK 

Subscript  of  UN  representing  prior  stale  number 

IJN(IJK) 

Mode  sequen  e number  for  the  indicated  system  state 

INCEI 

Number  of  arbsystems  read-in 

IU 

k 

Subscript  < f JU  and  Kl'  arrays  representing  subsystem 
ordinal  nui  iber 

JCAP 

Ultimately  the  number  of  mode  sequences 

JCOUNT 

K 

Ultimately  the  number  of  system  states  generated 

JES(L.K,J) 

Final  binary  array  defining  the  subfunction  and  mode 
of  the  Jth  mode  sequence  "■ 

JU(IU,J) 

Subsystem  ordinal  numbers  used  in  the  Jth  mode 
sequence 

4 

KS(L.K) 


Temporary  binary  array  defining  the  subfunction/mode 
number  in  the  current  mode  sequence 
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List  of  Abbreviations  and  Symbols  (contd) 
(Simulation  Model  — Program  One) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

KU(1U) 

Temporary  binary  array  defining  the  ordinal  numbers 
of  subsystems  used  for  the  current  mode  sequence 

KW(I) 

An  array  consisting  of  zeros  (0)  and  ones  ( 1 ) describing 
the  state  of  the  system 

LCEI 

Allowable  number  of  subsystems  in  equipment 
configuration 

LCNb 

Dimensional  constant.  Allowable  number  of  equipment 
name  segments 

LE 

Number  of  subfunctions  in  the  current  event 

LEE(l) 

Number  of  modes  in  subfunction  1 

LEQ 

Synonym  for  LEE,  used  as  do-loop  limit 

LF 

Subfunction  ordinal  number 

LI 

Maximum  number  of  modes  (all  subfunctions) 

L40 

Dimensional  constant.  Maximum  number  of  sub- 
functions 

L72 

Dimensional  constant.  Maximum  number  of  mode 

L 

sequences 

LLF(K) 

Mode  ordinal  number  for  Kth  subfunction 

LLCKLK) 

Next  mode  number,  if  current  subfunction  and  mode 
can  be  accomplished 

LLQS 

Designated  value  for  the  next  mode  in  the  mode 
sequence  definition 

LL1 

Location  of  current  state  zero  no.  1 

LL2 

Location  of  current  state  zero  no.  2 

LL3 

Location  of  current  state  zero  no.  3 

LL4 

Location  of  current  state  zero  no.  4 

LL5 

Location  of  current  state  zero  to.  5 
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List  of  Abbreviations  and  Symbols  (contd) 
(Simulation  Model  Program  One) 


Abbreviation 
or  symbol 

Fquivalent  in 
mathematical 
model 

Definition 

LLb 

Location  of  current  state  zero  no.  b 

LL7 

Locution  of  current  state  zero  no.  7 

LL8 

Location  of  current  state  zero  no.  K 

LMA(l.K.L) 

The  binary  urray  defining  the  subsystem  numbers  (L) 
required  for  subfunction  1 and  mode  K 

LMAT 

Ordinal  numbers  of  subsystems  needed  for  the  I.J 
subfunction/mode  input  to  show  location  of  a (1)  in 
the  IMA  array 

LMATT(L) 

Temporary  storage  for  IMA 

LQ 

Next  subfunction  number  if  current  subfunction  and 
mode  equipment  requirements  are  met 

LQSVL 

Designated  value  for  the  next  subfunction  in  the  mode 
sequence  definition 

LTF.ST 

Flow  control  Hag  for  blocking  the  selection  of  the  next 
subfunction  if  the  equipment  requirements  arc  not  met 

MAV 

Ordinal  number  of  critical  subsystem  (one  of  three) 

M('R 

Ordinal  number  of  critical  subsystem 

MD(I,J) 

Required  mission  condition  for  subfunction  and  mode 

MDT 

Mission  descriptor  (T  or  F)  specifying  mission  condi- 
tions (cloud  cover,  visibility,  etc.)  for  the  event 

MI(MM) 

Flag  for  each  subsystem,  MM,  set  to  one  or  zero 
depending  on  ONOFF 

MM 

Subscript  of  Ml  representing  subsystem  ordinal  number 

MMQ(I.J) 

Next  mode  number  if  current  subfunction  and  mode 
do  not  meet  requirement  of  MDT 

MPR 

Ordinal  number  of  critical  subsystem 

MPKM 

m 

Subscript  for  HDMT  representing  the  number  of  a 
damage  mechanism 
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List  of  Abbreviations  and  Symbols  (contd.) 


(Simulation  Model  - Program  One) 


Abbreviation 
or  symbol 

equivalent  in 
mathematical 
model 

Definition 

Units 

MQll.K) 

Next  subfunction  number  if  current  subfunction  and 
mode  does  not  meet  requirements  of  the  “mission 
descriptor"  MDT 

NC 

A flag  which  identifies  birth  whether  the  current  event 
type  is  offensive,  defensive-nuclear  (also  identifies  the 
number  of  different  damage  mechanisms),  defensive- 
conventional,  or  vehicle  mode 

NZ 

Current  number  of  zeros  (0)  in  the  KW  array 

NZT 

Maximum  number  of  zeros  (0)  allowed  in  the  KW  array 
(the  maximum  allowable  by  the  program  is  8) 

OFF 

A word  containing  the  Holerith  characters 
“OFF!  NSIVF" 

ONOFF 

A (lag  used  to  include  or  exclude  a subsystem  from  the 
system 

OUT 

Mask  word  against  which  the  input  variable  ONOFF  is 
tested 

TBFM 

MTBF 

Mean  time  between  failures 

T2 

Time  of  current  event 

WFAPN(7) 

Description  of  weapon  used,  whether  event  is  offensive 
or  defensive 
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List  of  Abbreviations  and  Symbols 


(Simulation  Model  Program  Two) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

Units 

CEP(J) 

CEP 

Threat  circular  error  probable  for  Jth  mode 
sequence 

feet 

CN(3) 

Contains  the  holerith  characters  for  the 
words  OFFENSIVE  or  DEFENSIVE  used  in 
event  output  description 

CURVH(l.J) 

UM'.r) 

Damage  level  at  radius  J for  l-typc  damage 
mechanism 

DEF(N) 

Contains  the  holerith  characters  for  the 
word  DEFENSIVE 

DF.LTAT 

At 

Elapsed  subsystem  time  during  transition 
used  in  reliability  calculation 

hours 

;j 

DFNAME(8) 

Alphanumeric  subsystem  name  characters 

DENUM 

Denominator  of  exponent  in  calculation  of 
survivability  calculation 

feet- 

DEN  175 

Test  value  for  miss  distance  used  to  prevent 
underflow  in  zone  miss  probability  calcu- 
lation 

ET 

E-r(N) 

Expected  number  of  targets  killed,  cumu- 
lative 

EXM 

Temporary  value  for  second  term  in  nuclear 
zone  miss  probability 

EXI 

Temporary  value  for  first  term  in  q'  calcula- 
tion and  in  conventional  miss  probability 
calculation 

Ji 

• 

f| 

EX2 

Second  term  in  conventional  miss  proba- 
bility calculation 

EXM  1 

First  term  in  nuclear  zone  miss  probability 
calculation 

1 

I 

FA(J) 

SPd 

Expected  number  of  weapons  arriving  in 
target  vicinity  lor  Jth  mode  sequence 

J 

FTDNN(I.J) 

Probability  of  abort-detection  before  time  J, 
given  flight  time  equal  to  J 

I j 

I 

* j 
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List  of  Abbreviations  and  Symbols  (contd) 
(Simulation  Model  - Program  Two) 


Abbreviation 
or  symbol 

Equivalent  in 
inathemat  ical 
model 

Definition 

Units 

FTFNN(I) 

Probability  of  flight  time  less  than  1 

GRAPH(IJ) 

Damage  threshold  vs.  nuclear  weapon  mis 
data  being  read  into  interpolation  subroutine 

HDMT(I) 

Gamma  dot  damage  threshold 

rads/  sec 

HDMT(2) 

Neutron  damage  threshold 

ncut  rons/cm  * 

HDMT(3) 

Blast  damage  threshold 

lb/in~ 

HDMT(4) 

Thermal  damage  threshold 

calorics/in’ 

ICOUNT 

i 

Prior  event  state  number 

IEVLNT 

Prior  event  number  (normal  mission) 

II 

Prior  event  type: 

II  = 0 offensive 

I -5  defensive-nuclear 

6 defensive  quick  conventional  kills 

7 defensive  slow  kill 

III 

Location  of  prior  state  zero  no.  1 

112 

Location  of  prior  state  zero  no.  2 

113 

Location  of  prior  state  zero  no.  3 

114 

Location  of  prior  state  zero  no.  4 

IIS 

Location  of  prior  state  zero  no.  5 

116 

Location  of  prior  state  zero  no.  6 

117 

Location  of  prior  state  zero  no.  7 

118 

Location  of  prior  state  zero  no.  8 

IINCKI 

One-half  the  number  of  elements  in  the  state 
vector  (for  the  vehicle  mode  only) 

UK 

i 

Prior  state  number 

IJN(IJK) 

•1(0 

Mode  sequence  number  for  prior  state 
number  UK 
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List  of  Abbreviations  and  Symbols  (contd) 
(Simulation  Model  - Program  Two) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

INCE1 

Ultimately,  the  number  of  subsystems  read 
from  input  data 

INE 

One  more  than  INCEI 

1N1T 

'‘Next"  subsystem  number,  to  be  tested  for 
lethal  radius  against  "this"  subsystem 
number 

1 >S 

Temporary  value  of  current  damage  mecha- 
nism type 

ISUB(I) 

Temporary  damage  mechanism  type  number 
for  1th  damage  mechanism 

IT(M) 

Prior  state  binary  bit  (before  reordering)  for 
Mth  subsystem 

ITDN 

r\ 

Prior  event  abort-detection  time,  referenced 
to  current  event  time 

1TFN 

‘i 

Prior  event  flight  time,  referenced  to  current 
event  time 

1TORE 

Temporary  value  of  IZ(L) 

ITT(M) 

Reordered  prior  state  binary  bit  for  sub- 
system M 

1TTD 

Prior  event  abort-detection  time,  referenced 
to  prior  event  time 

IV(L) 

Prior  event  vulnerability  index  number  for 
Lth  subsystem 

IZ(L) 

IZ(k) 

Current  event  vulnerability  index  number 
for  subsystem  L 

JCAP 

Number  of  mode  sequence  read  from  IJN 
tape 

JCOUNT 

Current  event  state  number 

JEVENT 

Event  number  in  normal  mission  used  to 
initialize  state  probability  for  an  abort 
mission 
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Abbreviation 
or  symbol 

J1NCE1 

JV 

KM  AX 
KM  AX  I 

KW(M) 

KWW(I) 

L72 

LCEI 

LCN8 

LL1 

LL2 

LU 

LL4 

LL5 

LL6 

LL7 

LU 

LMAT(L) 


List  of  Abbreviations  and  Symbols  (contd) 
(Simulation  Mcxlel  - Program  Two) 


Equivalent  in 
mathematical 
model 


Definition 


One  plus  one-half  the  number  of  elements  in 
the  state  vector  (for  the  vehicle  mode  only) 

Temporary  value  for  IV(M) 

Number  of  different  offset  zones 

One  less  than  the  number  of  offset  weapon 
zone  radii  (conventional  weapons) 

Current  state  binary  bit  for  subsystem  M 

Reordered  current  state  binary  bit  for 
subsystem  I 

Maximum  number  of  mode  sequences  gener- 
ated by  program  one  for  any  event  (40) 

Allowable  number  of  subsystems  in  the 
equipment  configuration 

Allowable  number  of  equipment  name 
segments 

Location  of  current  state  zero  no.  I 

Location  of  current  state  zero  no.  2 

Location  of  current  state  zero  no.  .f 

Location  of  current  state  zero  no.  4 

Location  of  current  state  zero  no.  5 

Location  of  current  state  zero  no.  6 

Location  of  current  state  zero  no.  7 

Location  of  current  stale  zero  no.  8 

Ordinal  number  of  subsystem  required  (elec- 
tronics mode)  or  location  of  “one"  in  state 
vector  (vehicle  mode) 


Units 


I 


LMAX 


Number  of  different  elevation  angles  associ- 
ated with  warhead  offset  trajectory 
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List  of  Abbreviations  and  Symbols  (contd) 


(Simulation  Model  - Program  Two) 


Abbreviation 
or  symbol 

Hquivalcnt  in 
mathematical 
model 

Definition 

Units 

LOO 

“This”  subsystem  number  used  in  defining 
vulnerability  index  numbers 

M 

K 

Subsystem  ordinal  number 

MAV 

Ordinal  number  of  critical  subsystem 

MCR 

Ordinal  number  of  critical  subsystem 

MDT 

Mission  descriptor: 

T = scenario  constraint  imposed 

F = scenario  constraint  not  imposed 

MI(M) 

On/otT  Hag  for  subsystem  M 

MLTH 

Mission  length  measured  from  initiation  of 
first  threat  exposure 

MPR 

Ordinal  number  of  critical  subsystem 

MPRM 

9 

m 

Damage  mechanism  type  number 

MTBF(M) 

MTBF 

Mean  time  between  failures 

hours 

NABORT 

Number  of  the  event  at  which  the  abort  path 
is  initialed 

N(' 

A flag: 

NC  * 0 offensive  event 

NC  = 1,2,  3,  4,  or  5 nuclear  defensive 
event  with  NC  damage  mechanisms 

NC  * 0 conventional  defensive  event 

NC  « 7 vehicle  simulation  event 

NF 

A flag 

NF  * 0 nuclear  damage  is  possible 

NF  * 1 nuclear  damage  is  not  possible 

NPOINT 

Number  of  points  in  the  miss  distance  vs. 
damage  '*vel  table  lookup 

NT  ON 

*) 

Current  event  >bort-detection  time,  refer- 
enced to  current  e*c"'  time 

150 


JTCG/A.S-76-S-004 


I. 


I 

List  ot'  Abbreviations  and  Symbols  tcontd) 


(Simulation  Model  Program  Two) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

Units 

NTFN 

s 

Current  event  flight  time,  referenced  to 
current  event  time 

NTTD 

Current  event  abort-detection  time,  refer- 
enced to  prior  event  time 

NTTF 

Current  event  flight  time,  referenced  to  prior 
event  time 

NZ 

A counter  to  keep  track  of  current  number 
of  zeros  in  prior  state  vector 

NZI 

A counter  to  keep  track  of  current  number 
of  zeros  in  current  state  vector 

NZT 

Maximum  number  of  zeros  allowed  in  state 
vector  for  the  mission 

N2 

Two  times  the  number  of  lethal  radii 

feet 

OFF(N) 

Contains  the  Holerith  characters  for  the 
word  OFFENSIVE 

ONOFF 

Flag  used  to  include  or  exclude  a subsystem 
from  the  system  (T  = include.  F = exclude) 

OUT 

Variable  against  which  input  variable 
ONOFF  is  tested 

PARIVE 

Probability  of  arrival 

PCAP(J) 

Probability  of  Jth  mode  sequence 

PCKILL(M.L.K) 

k 

Probability  of  kill  of  subsystem  M given 
offset  zone  K and  elevation  L (conventional) 

PCSURV(M.LK) 

vk 

Probability  of  survival  of  subsystem  M 
(conventional) 

PI(IJK) 

Pi 

Probability  for  prior  state  no.  UK 

PJ(L) 

Pi.n 

Probability  for  current  state  no.  L 

PK(J) 

Weapon  delivery  kill  probability  for  Jth 
mode  sequence 
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List  of  Abbreviations  and  Symbols  (contd) 


(Simulation  Model  - Program  Two) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

Units 

PM(M) 

dk 

Probability  of  reliable  operation  in  transition 
(before  reordering) 

PMISS(K,IJK) 

PmissM 

Probability  of  threat  miss  distance  in 
zone  K,  given  state  IJK 

PMM(M) 

dk 

Reliability  of  each  subsystem  M in  transition 
(after  reordering) 

PTDN(1,J) 

Probability  that  abort-detection  time  is  J, 
given  that  flight  time  is  I 

PTFN(I) 

Probability  that  flight  time  is  I 

QO(IJK) 

q' 

The  probability  of  nuclear  weapon  miss 
distance  exceeding  the  maximum  lethal 
radius 

QPRM(M.IJK) 

r 

qk 

Probability  of  nuclear  weapon  miss  in  zone 
associated  with  Mth  subsystem,  given  the 
IJKth  state 

R(K) 

Rk 

Inner  radius  of  the  Kth  offset  zone 

feet 

RA(M) 

Temporary  value  of  subsystem  lethal  radius 
obtained  from  linear  interpolation  routine 

feet 

RESUL(2) 

rk 

Lethal  radius  output  of  interpolation 
subroutine 

feet 

RM2 

rk2 

Square  of  the  radius  RM(1)  or  R(k) 

feet^ 

RMM12 

Square  radius  associated  with  second  term  in 
nuclear  zone  miss  probability  calculation 

feet  2 

RMM2 

Square  of  the  radius  RM(M) 

feet  2 

SIGMA 

0 

Standard  deviation  of  miss  distance  whether 
nuclear  or  convention  weapon 

feet 

STORE 

rk 

Temporary  value  of  RM(k) 

feet 

SUM 

**k,n 

Average  (over  all  modes)  target  kill  proba- 
bility in  nth  event 

JTCG/AS-7  6-S-004 


w 


1 

r 

1. 


i. 


p 

1. 


i. 


List  of  Abbreviations  and  Symbols  (contd) 


(Simulation  Model  - Program  Two) 


Abbreviation 
or  symbol 

Equivalent  in 
mathematical 
model 

Definition 

Units 

SUMPTD 

Ultimately,  the  probability  of  abort- 
detection  time  (due  to  an  encounter)  greater 
than  or  equal  to  the  current  abort  detection 
time 

SUMTRA 

Ultimately,  the  sum  of  all  transition  proba- 
bilities resulting  from  all  LMAX  offset 
trajectories  for  one  zone 

SUMTRK 

Ultimately,  the  sum  of  all  transition  proba- 
bilities resulting  from  all  KMAX  offset 
trajectories 

Tl 

Time  of  prior  event 

hours 

T2 

Time  of  current  event 

hours 

TBFM 

MTBF 

Mean  time  between  failures 

hours 

THDM(M.IS) 

Temporary  value  of  damage  threshold  for 
subsystem  M.  damage  type  IS 

(See  IIDMT) 

TIMEF(M) 

Time  that  subsystem  M is  turned  off 

hours 

TIMEN(M) 

Time  that  subsystem  M is  turned  on 

hours 

TMF 

Time  that  subsystem  is  turned  off 

hours 

TMN 

Time  that  subsystem  is  turned  on 

hours 

TRANS 

Pj/i,n 

Transition  probability  for  state  pair  under 
consideration 

TRANS1 

Temporary  value  of  TRANS  in  the  vehicle 
transition  algorithm  (slow  threat  damage) 

TWODNZ 

Name  of  subroutine  for  two-dimensional 
interpolation 

WEAPN(7) 

Description  of  weapon  whether  offensive  or 
threat 
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Subscript  used  to  fill  answer  vector 


Current  subscript  of  upper  value  of  result 
interval 

Subscript  of  lower  value  in  result  interval 

Subscript  of  upper  value  in  result  interval 

Number  of  values  in  any  column  (same  for 
all  parameters) 

Number  of  different  parameters 

Given  value  of  independent  variable 

Name  of  table  to  perform  lookup  on 

Subscript  of  first  value  in  independent 
variable  column 

Column  of  independent  variable 
Ratio  factor  for  linear  interpolation 
Answer  vector 

Subroutine  name  and  acronym  for  two- 
dimensional  linear  interpolation 
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